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Abstract 



The non-classical, non-local quantum correlations so-called entanglement, which arc possibly 
the greatest mystery of quantum mechanics, have been proved to be a very useful resource. 
Entanglement may be used to make exponentially faster computations with a quantum 
computer than the equivalent ones with a classical computer. Also, it has been shown 
to introduce advantages in the quantum communication protocols over the classical ones. 
Moreover, nowadays it is possible to codify provably secure messages through the quantum 
cryptography protocols (which have a deep relationship with entanglement), in opposition 
to the classical protocols, which are just conditionally secure. 

This Thesis is devoted to the analysis of entanglement in relevant physical systems. 
Entanglement is the conducting theme of this research, though I do not dedicate to a single 
topic, but consider a wide scope of physical situations. 

I have followed mainly three lines of research for this Thesis, with a series of different 
works each, which are: 

• Entanglement and Relativistic Quantum Theory. 

I show the unbounded entanglement growth appearing in S-matrix theory of scattering 
for incident fermions with sharp momentum distributions, in the context of quantum 
electrodynamics. I study the possibility of spin entanglement generation through spin- 
independent scattering of identical particles. I also analyze the properties of Lorentz 
invariance and relativity of entanglement distillability and separability. I propose the 
complete simulation of Dirac equation and its remarkable relativistic quantum effects 
like Zitterbewegung or Klein's paradox in a single trapped ion. 

• Continuous-variable entanglement. 

I demonstrate that an arbitrary, unbounded degree of entanglement may be achieved 

between two atoms by measurements on the light they emit, when taking into account 
additional ancillary photons. I detect the spin entanglement loss due to transfer of cor- 
relations to the momentum degree of freedom of s = ^ fermions or photons, through 



local interactions entangling spin and momentnm. I also develop a mathematical 
method for computing analytic approximations of the Schmidt modes of a bipartite 
amplitude with continuous variables. I study the momentum entanglement gener- 
ated in the decay of unstable systems and verify that, surprisingly, the asymptotic 
entanglement is smaller for wider decay widths, related to stronger interactions. 

• Multipartite entanglement. 

After a careful analysis of the 1 ^ M approximate quantum cloning for qubits se- 
quentially implemented, I show that it can be done with just linear resources of the 
ancilla: the dimension of the ancilla Hilbcrt space grows linearly with the number 
of clones, while for arbitrary multiqubit states sequentially generated it would grow 
exponentially in the number of qubits. This has remarkable experimental interest 
as it provides a procedure for reducing approximate quantum cloning to sequential 
two-body interactions, which are the ones experimentally feasible in laboratory. I 
also propose an inductive scheme for the classification of A?^-partite entanglement, for 
arbitrary N, under stochastic local operations and classical communication, based on 
the analysis of the coefficient matrix of the pure state in an arbitrary product basis. 
I give the complete classification of genuine 4-qubit entanglement. 
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Chapter 1 

Introduction 



/ would not call that one but rather the characteristic trait of quantum me- 
chanics 

-Erwin Schrodinger jSch35| 

1.1 Motivation 

Quantum computation and quantum information is the field that deals with the information 
processing and transmission making use of quantum systems. This is a discipline that has 
begun its development in recent years, and to which a lot of effort is being devoted (for 
a review of the field, see for example |GMD02] INCOOj ). Its interest stems from the fact 
that there are physical tasks regarding information processing and transmission that can be 
possibly done more efficiently with quantum systems than with the ordinary classical ones. 
Some examples are the following: 

1. Regarding information processing we highlight 

• The factorization of large integer numbers into primes [Sho97] . 

This takes a time that grows exponentially in the number of digits N of the 
integer being factorized with the most powerful classical algorithms known. That 
is the property in which many classical cryptography systems are based, i.e. 
public key distribution. With Shor's algorithm a quantum computer would only 
use a polynomial time in A^. 

• Search algorithms |Gro97] . 

"Quantum mechanics helps in searching for a needle in a haystack" (L. Grover). 
With a classical computer it takes a time of order N to find a given item in a 
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disorganized list of size N. In a quantum computer it would take a time that 
grows as ^/N, with Grover's algorithm. 

2. Regarding information transmission, we mention 

• Quantum cryptography [Wie83|, |BBB"'"92] . 

The private distribution of secure keys, for use in cryptography, is a reality. 
Nowadays there are some companies which offer this service. This feature relies 
in the impossibility of distinguishing with certainty two non-orthogonal quantum 
states. 

• Dense coding |BW92] . 

With quantum mechanics it is possible to transmit two bits of information with 
a single quantum bit, in a secure way. 

• Quantum teleportation |BBC+93L IbBM+981 lBPM+97j . 

This is a feature very useful in quantum information. It is a procedure to move 
quantum states around without making use of a communication channel. With 
this protocol a copy of a certain quantum state is obtained in another specified 
place. This implies the destruction of the initial state, because the cloning of 
quantum states is forbidden by quantum linearity (no-cloning theorem). 

Some of the previous applications make use of a genuine quantum resource, namely En- 
tanglement. In the early days of quantum mechanics, Erwin Schrodinger already realized 
the importance of this quantum property. He said, "Entanglement is the characteristic trait 
of Quantum Mechanics". Roughly speaking, entanglement are genuine quantum correla- 
tions between spatially separated physical systems. To show the difference between these 
non-local, quantum correlations and purely classical ones, we consider an example conceived 
by Asher Peres |Per78t lPer95j : Suppose a bomb, initially at rest, which explodes into two 
fragments carrying opposite momenta Ji, J2 = — Ji (see FIG. II. ip . An observer measures 
the magnitude sign(Q; • Ji), where d is a unit vector with a fixed arbitrary direction. The 
result of the measurement, called a, is +1 or —1. Additionally, a second observer measures 
sign(/3 • J2), where /3 is another unit vector with another fixed arbitrary direction. The 
result can only be b = ±1. 

The experiment is repeated times. We call aj and bj the results measured by these 
observers for the jth bomb. If the observers compare their results, they find a correlation 

(ab) =^ajbj/N. 
3 
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Figure 1.1: Gedanken experiment in which a bomb, initially at rest, explodes into two 
fragments with opposite momenta, J2 = — Ji- 

For example, if a = (3, they obtain (ab) = — 1. 

To compute {ab) for arbitrary a and /?, we consider the sphere shown in FIG. II. 2[ The 
plane orthogonal to a divides the sphere in two hemispheres. We have a = 1 if Ji points 
through one of these hemispheres and a = —1 if it points through the other hemisphere. 
Similarly, the regions where 6 = ±1 are limited by the intersection of the sphere with the 
plane orthogonal to /3. This way, the sphere is divided into four sections, as shown in FIG. 
II. 2[ The shaded sections have (ab) = 1, while the unshaded ones have (ab) = —1. The 
classical correlation for uniformly distributed Ji results 

{ab) = [e-{Tr-9)]/7r = -1 + 29/71. (1.1) 

Regarding now the quantum mechanical case, we consider two spatially separated spin 1/2 
particles 1 and 2 in the singlet state 

l^) = ^(l Tii2)-UiT2)), 

where the arrows denote the third component of spin along an arbitrary direction. An 
observer measures the observable a-ai, while another one measures [3 ■ a2, being ai and (T2 
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Figure 1.2: Auxiliary figure for computing tlie classical correlation (jl.ip . The shaded areas 
correspond to (ab) = 1, while the unshaded ones correspond to (ab) = — 1. 

the Pauli spin matrices associated to particles 1 and 2 respectively, a and /3 are arbitrary 
unit vectors. We denote as before a and b the results of these measures, taking values ±1. 
It can be shown that in this case, the correlation is 

(ab) = -a- (3 = -cose. (1.2) 

We plot (jl.ip and ()1.2p in FIG. 11.31 This figure shows how the quantum correlation is 
always stronger than the classical one, except where both are or ±1. This qualitatively 
different behavior of the quantum and classical correlations has very profound implications, 
as John Bell shown |Bel64] . upon the arguments of Einstein, Podolsky and Rosen |EPR35j : 
realism+locality is incompatible with quantum mechanics. Up to now, experiments favor 
the latter. There is, indeed, a general consensus about the validity of quantum mechanics, 
and, in particular, entanglement, in opposition to local realistic theories. Although it is 
difficult to prove with total security the correctness of quantum mechanics, there are no 
relevant experiments that contradict i10. 



^For a recent reference showing the allowed correlations volume predicted by local realistic theories, 
quantum mechanics, and in general no-signalling theories, see |Cab05| 



1.1 Motivation 




Figure 1.3: Quantum correlation (jl.2p (solid) and classical one (jl.ip (dashed), (ab), as 
functions of 6. 

1.1.1 Preliminaries: Basic notions of entanglement 

In this section we give some relevant definitions about entanglement, that will be used along 
the Thesis. 

We consider a composite system S described by a Hilbert space 7i, either finite- or 
infinite-dimensional. This space is built upon the tensor products of the Hilbert spaces 
associated to the subsystems Sa of 5. We consider for the time being bipartite systems S, 
for simplicity. Thus, a = 1, 2, and TC = TCi 'H2- 

• Pure states 

Definition 1. (product state) A vector state \^) of system S is a product state if it 
can be expressed as 

1^) = |${i))|^(2)), (1.3) 
where l^'(^)) G Hi and j^f^^)) £ ■H2. 

Definition 2. (entangled state) A vector state \^) of system S is entangled if it is 
not a product state. 

A remarkable example, for dim('Hi) = dim('H2) = 2, are the so-called Bell states. 































-I^«)l^?)), 


(1.4) 
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that possess the maximum achievable entanglement for these dimensions (1 ebit, or 
entangled bit). 

A very useful tool for analyzing the entanglement of pure bipartite states is the 
Schmidt decomposition [EKQSj iLWEOOj . It basically consists in expressing the pure 
bipartite state as sum of biorthonormal products, with positive coefficients, ^/X^, 
according to 

n=0 

where {|\l'n ^)}, are orthonormal bases associated to 7ii and 7i2, respectively 

(see Appendix E]) . In Eq. (jl.Sp . d = min{dim('Hi), dim('H2)}; and it can be infinite, 
like for systems described with continuous variables as momentum, energy, position, 
frecuency, or the like. In those cases, the states would be wave functions, 

{p\^t^) = i^^^Hp), « = 1,2, (1.6) 
where p denotes the corresponding continuous variable. 

For pure bipartite states the relevant entanglement measure is the entropy of entan- 
glement, S. Given a certain state |^), it is defined as the von Neumann entropy of 
the reduced density matrix with respect to 5i o ^2, 

d-l 

S ■.= -Y,Klog2Xn, (1.7) 

n=0 

where the A„ coefficients are the eigenvalues of the reduced density matrix of |^) with 
respect to either of the two subsystems, and are the ones appearing in Eq. p.5|) . In 
general, 5 > 0. S = to. a product statfl and, the more entangled is a state, the 



larger is S. For a maximally entangled statqj, S = log2 d, and if d = oo, then S 
diverges. 

Another interesting entanglement measure for pure states is the Schmidt number, K. 
It is defined 

K:=-J—^. (1.8) 



K gives the effective number of terms appearing in the Schmidt decomposition (jl.Sp 
of a pure bipartite amplitude. K = 1 iov product states, and, the larger the 
larger the entanglement. Along this Thesis we will be using the notation K both for 



In this case, Aq = 1 and A,i = 0, n > 0. 



2 

^In this case, \„ — n — d ~ 1. 
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the Schmidt or the Slater number. The latter is defined analogously, although it is 
related to pure bipartite amplitudes of identical fermions, that can always be written in 
terms of superpositions of biorthonormal Slater determinants (it is the generalization 
of the Schmidt decomposition for identical fermions, called Slater decomposition). 
The Slater number gives the effective number of Slater determinants appearing in the 
Slater decomposition. This is just half the Schmidt number. On the other hand, the 
degree of genuine entanglement is the same in both cases, given that for identical 
fermions the correlations due to antisymmetrization must be substracted: they are 
not genuine entanglement. 

• Mixed states 

Definition 3. (separable state) A separable state can always be expressed as a convex 
sum of product density operators fWerS^. In particular, a separable bipartite state 
can be written as 

p = Y,Qpt^^pf'\ (1.9) 

i 

where C-i > 0, Cj = 1, and p^f* and pf''* are density operators associated to sub- 
systems A and B. 

Definition 4. (entangled state) An entangled mixed state is a quantum state that is 
not separable. 

A remarkable example are the so-called Werner states, canonical examples of mixed 
states obtained from a Bell state that suffers decoherence. They are defined in the 
way 

p^B F\^'){^-\ + (^|^+)(^+| + |^")($"| + |«>+)($+|^ . (1.10) 

Where < F < 1, and F gived the degree of mixture of the state. It is well-known 
|Wer89j that the Bell state 1^") is distillablfQ from Eq. (frTO|) iff F > 1/2. 

For mixed bipartite or multipartite states there are not known universal entanglement 
measures. In fact, neither there are criteria that may allow to determine whether a 
certain state is or not entangled (separability criteria), and they are just known in 
some particular cases. 

For dim(Wi) = dim(7^2) = 2 or dim(7^i) = 2 and dim(W2) = 3 it exists the PPT 
separability criterion [Per 961 IHHH96j . It establishes that a certain bipartite mixed 

^it is, may be obtained from a certain number of copies of pg^ by Local Operations and Classical 
Communication. 
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state is separable iff its partial transposed matrix (PT, transposed with respect to 
one of the two subsystems) is positive (with positive eigenvalues). Otherwise it is 
entangled. 

This criterion gives rise to defining an entanglement measure for 2x2 and 2x3 
dimensions, so-called negativity N, according to 



where Amin is the smallest of the eigenvalues of the PT matrix. A separable state has 
= 0, an entangled state, N > 0, and, the more entangled is a state, the larger A^. 

For the sake of completeness we mention that there exist many more entanglement 
measures for low dimensions (Concurrence, Entanglement of Formation, Tangle, etc.) 
and other separability criteria for arbitrary dimensions like the ones based on entan- 
glement witnesses [LBC+OOllPVOT] . 

At the same time that the field of quantum information and computation evolves, more 
and more evidence appears that shows the crucial role of entanglement in this field. Bi- 
partite and multipartite entanglement is one of the features that give rise to many of the 
developments of quantum computation and information, like quantum cryptography |Wie83t 
lEkeOlj . dense coding |BW92j . quantum teleportation [BBC+931 [BBM+98[ [BPM+97| and 
aspects of quantum computation [NCOOj . among others. A lot of effort is being devoted 
to obtain separability criteria (to decide whether a given mixed or pure state is entangled 
or not), and to measure and characterize entanglement. For a review, see Ref. |LBC+00] 
or Ref. |PV07j . For a compilation of bibliographic references on entanglement and other 
topics on foundations of quantum mechanics and quantum information, see Ref. [CabOO] . 
The evaluation of the entanglement of a composite state is thus a main task to be done. 

This Thesis is mainly a series of theoretical results about entanglement. The impor- 
tance of entanglement, both from the theoretical, fundamental point of view, and for the 
experimental applications in information processing and communication, is a well-founded 
motivation for having carried out the lines of research I have developed here. With respect 
to the focusing, being this a theoretical Thesis, it is mainly related to entanglement prop- 
erties from the physical point of view. This is indeed a Thesis about physics, more than 
about mathematics or computer science. However, Part IIIII (Multipartite entanglement) in 
this Thesis is more related to mathematics and information theory, with two more abstract 
chapters. 

In the following I briefly expose the main results of this Thesis. 



A^ := max{0, — 2Amm} 
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1.2 Contributions 

The research hnes I fonowed for carrying out this work are mainly three: 

• Entanglement and Special Relativity. 

• Entanglement of pure states described by continuous variables. 

• Multipartite entanglement. 



1.2.1 Entanglement and Relativistic Quantum Theory 

Peres, Scudo y Terno introduced [PST02j the relative character under Lorentz transforma- 
tions of the entropy of entanglement of a s = ^ particle. This result may have profound 
implications given that the quantum information theory is being mainly developed in the 
laboratory frame, without considering different observers in relative uniform motion, so it 
is not covariant. It would be interesting to investigate what happens when processing or 
transmitting quantum information in relativistic regimes. In this respect several papers 
have appeared, which explore the relationship among quantum information theory and spe- 
cial relativity |Cza97[ IPST02[ IAM021 IGA021 IGBA031 [PS03l ITTMI |AjLMH03i iPTOl IMY04[ 



IFSMn5[ IHarn5[ lAFSMTOGl IHarOBbi IHiHIHil IHWOGI I.TSSOBbl I.TSSnOa] 

Following this line we have contributed with 



• Dynamics of momentum entanglement in lowest-order QED 

This is a work [LLS06b] somewhere in between the entanglement of continuous vari- 
ables and the relativistic aspects of entanglement, so it could be placed in either part 
of the Thesis. Here we study the dynamics of momentum entanglement generated 
in the lowest order QED interaction between two massive spin-1/2 charged parti- 
cles, which grows in time as the two fermions exchange virtual photons. We observe 
that the degree of generated entanglement between interacting particles with initial 
well-defined momentum can be infinite. We explain this divergence in the context of 
entanglement theory for continuous variables, and show how to circumvent this appar- 
ent paradox. Finally, we discuss two different possibilities of transforming momentum 
into spin entanglement, through dynamical operations or through Lorentz boosts. 

• Generation of spin entanglement via spin-independent scattering 

Here we consider |LL06| the bipartite spin entanglement between two identical fermions 
generated in spin-independent scattering. We show how the spatial degrees of freedom 
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act as ancillas for the creation of entanglement to a degree that depends on the scat- 
tering angle, 9. The number of Slater determinants generated in the process is greater 
than 1, corresponding to genuine quantum correlations between the identical fermions. 
The maximal entanglement attainable of 1 ebit is reached at = tx 12. We also analyze 
a simple Q dependent Bell's inequality, which is violated for 7r/4 < < tt/2. This 
phenomenon is unrelated to the symmetrization postulate but does not appear for 
unequal particles. 

• Relativity of distillability 

In this work we study |LMDS06] entanglement distillability of bipartite mixed spin 
states under Wigner rotations induced by Lorentz transformations. We define weak 
and strong criteria for relativistic isoentangled and isodistillable states to characterize 
relative and invariant behavior of entanglement and distillability. We exemplify these 
criteria in the context of Werner states, where fully analytical methods can be achieved 
and all relevant cases presented. 

• Dirac equation and relativistic effects in a single trapped ion 

We present |LLSS07b] a method of simulating the Dirac equation, a quantum-relativistic 
wave equation for spin-1/2 massive particles, in a single trapped ion. The four- 
component Dirac bispinor is represented by four metastable ionic internal states, 
which, together with the motional degrees of freedom, could be controlled and mea- 
sured. We show that the proposed scheme would allow for a smooth transition from 
massless to massive particles, as well as for access to parameter ranges and physical 
regimes not provided by nature. Furthermore, we demonstrate that paradigmatic 
quantum relativistic effects unaccesible to experimental verification in real fermions, 
like Zitterbewegung, Klein's paradox, Wigner rotations, and spontaneous symmetry 
breaking produced by a Higgs boson, could be studied. 

1.2.2 Continuous variable entanglement 

The entanglement of continuous variables has raised a lot of interest in the past years 
[Vai94l lFSB+981 [LB991 IGieOll IGECP031 EEOSl IBvL05| . For a thorough review of the field, 
see |BvL05| . 

We will concentrate in the continuous variable entanglement of pure bipartite states, 
which is very relevant for the applications and corresponds to the ideal case with no deco- 
herence. 

The results we have obtained in this line are 
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• How much entanglement can be generated between two atoms by detecting 
photons? 

We prove [LGRCOT] that in experiments with two atoms an arbitrary degree of entan- 
glement between them may be reached, by only using linear optics and postselection 
on the light they emit, when taking into account additional photons as ancillas. This 
is in contrast to all current experimental proposals for entangling two atoms, that 
were only able to obtain one ebit. 

• Spin entanglement loss by local correlation transfer to the momentum 

We show [LLS06a] the decrease of spin-spin entanglement between two s = fermions 
or two photons due to local transfer of correlations from the spin to the momentum 
degree of freedom of one of the two particles. We explicitly show how this phenomenon 
operates in the case where one of the two fermions (photons) passes through a local 
homogeneous magnetic field (optically-active medium), losing its spin correlations 
with the other particle. 

• Schmidt decomposition with complete sets of orthonormal functions 

We develop [LL05a] a mathematical method for computing analytic approximations of 
the Schmidt modes of a bipartite amplitude with continuous variables. In the existing 
literature, various authors compute the Schmidt decomposition in the continuous 
case by discretizing the corresponding integral equations. We maintain the analytical 
character of the amplitude by using complete sets of orthonormal functions. We give 
criteria for the convergence control and analyze the efficiency of the method comparing 
it with previous results in the literature related to entanglement of biphotons via 
parametric down-conversion. 

• Momentum entanglement in unstable systems 

We analyze the dynamical generation of momentum entanglement in the decay of 
unstable non-elementary systems described by a decay width T |LL05bj . We study 
the degree of entanglement as a function of time and as a function of T. We verify 
that, as expected, the entanglement grows with time until reaching an asymptotic 
maximum, while, the wider the decay width T, the lesser the asymptotic attainable 
entanglement. This is a surprising result, because a wider width is associated to a 
stronger interaction that would presumably create more entanglement. However, we 
explain this result as a consequence of the fact that for wider width the mean life 
is shorter, so that the system evolves faster (during a shorter period) and can reach 
lesser entanglement than with longer mean lives. 
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1.2.3 Multipartite entanglement 

The multipartite entangled states stand up as the most versatile and powerful tool for realiz- 
ing information processing protocols in quantum information science [BDOOj . The controlled 
generation of these states becomes a central issue when implementing the applications. In 
this respect, the sequential generation proposed by Schon et al. [SSV"'"05l ISHW^OG) is a 
very promising scheme to create these multipartite entangled states. 
We have contributed to the field with 



Sequential quantum cloning 

Not every unitary operation upon a set of qubits may be implemented sequencially 
through successive interactions between each qubit and an ancilla. Here we analyze 
|DLL"'"06 the operations associated to the quantum cloning sequentially implemented. 



We show that surprisingly the resources (Hilbert space dimension D) of the ancilla 
grow just linearly with the number of clones M to obtain. Specifically, for universal 
symmetric quantum cloning we obtain D = 2M and for symmetric phase covariant 
quantum cloning, D = M + 1. Moreover, we obtain for both cases the isometrics for 
the qubit-ancilla interaction in each step of the sequential procedure. This proposal is 
easily generalizable to every quantum cloning protocol, and is very relevant from the 
experimental point of view: three-body interactions are very difficult to implement in 
the laboratory, so it is fundamental to reduce the protocols to sequential operations, 
which are mainly two-body interactions. 



• Inductive classification of multipartite entanglement under SLOCC 

Here we propose [LLSSOGj ILLSS07a] an inductive procedure to classify A^-partite 
entanglement under stochastic local operations and classical communication (SLOCC) 
when the classification for — 1 qubits is supposed to be known. The method relies 
in the analysis of the coefficients matrix of the state in an arbitrary product basis. We 
illustrate this method in detail with the well-known bi- and tripartite cases, obtaining 
as a by-product a systematic criterion to establish the entanglement class of a pure 
state without using entanglement measures, in opposition to what has been done up 
to now. The general case is proved by induction, allowing us to obtain un upper 
bound for the number of entanglement classes of A^-partite entanglement in terms of 
the number of classes for — 1 qubits. Finally, we give our explicit calculation for 
the highly nontrivial case A^ = 4 [LLSS07aj . 



1.3 Description of the Thesis 
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1.3 Description of the Thesis 

• In Chapter [2] I analyze the momentum entanglement generation among two elec- 
trons which interact in QED by exchanging virtual photons. I show that surprisingly, 
S matrix theory produces pathological results in this case: the entanglement in M0ller 
scattering would be divergent for incident particles with well-defined momentum. In 
order to manage with these divergences, that would be physical (entanglement is a 
measurable magnitude, with a physical meaning), I made the calculation for electrons 
with Gaussian momentum distributions which interact for a finite time. The diver- 
gences disappear, but, remarkably, the attainable entanglement would not be bounded 
from above. 

• In Chapter[3]I consider the spin entanglement among two or more identical particles, 
generated in spin-independent scattering. I show how the spatial degrees of freedom 
act as ancillas creating entanglement between the spins to a degree that will depend in 
general on the specific scattering geometry considered. This is genuine entanglement 
among identical particles as the correlations are larger than merely those related to 
antisymmetrization. I analize specifically the bipartite and tripartite case, showing 
also the degree of violation of Bell's inequality as a function of the scattering angle. 
This phenomenon is unrelated to the symmetrization postulate but does not appear 
for unlike particles. 

• In Chapter |4] I analyze the Lorentz invariance of usual magnitudes in quantum infor- 
mation, like the degree of entanglement or the entanglement distillability. I introduce 
the concepts of relativistic weak and strong isoentangled and isodistillable states that 
will help to clarify the role of Special Relativity in the quantum information theory. 
One of the most astonishing results in this work is the fact that the very separability 
or distillability concepts do not have a Lorentz-invariant meaning. This means that a 
state which is entangled (distillable) for one observer may be separable (nondistillable) 
for another one that propagates with a finite v < c speed with respect the first one. 
This is an all- versus- nothing result, in opposition to previous results on relativistic 
quantum information, which showed that a certain entanglement measure was not 
relativistically- invariant (but always remained larger than zero). 

• In Chapter [5]I present a method for simulating Dirac equation, a quantum-relativistic 
wave equation for massive, spin-i particles, in a single trapped ion. The four- 
component Dirac bispinor is represented by four metastable, internal, ionic states, 
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which, together with the motional degrees of freedom, could be controlled and mea- 
sured. I show that paradigmatic effects of relativistic quantum mechanics unaccesible 
to experimental verification in real fermions, like Zitterbewegung, Klein's paradox, 
Wigner rotations, and spontaneous symmetry breaking produced by a Higgs boson, 
could be studied. 

• In Chapter [6] I prove that in experiments with two atoms an arbitrary degree of 
entanglement between them may be reached, by only using linear optics and postse- 
lection on the light they emit, when taking into account additional photons as ancillas. 
This is in contrast to all current experimental proposals for entangling two atoms, that 
were only able to obtain one ebit. 

• In Chapter \7\ I show the decrease of the initial spin-spin entanglement among two 
s = ^ fermions or two photons, due to local correlation transfer from the spin to the 
momentum degree of freedom of one of the two particles. I explicitly show how this 
phenomenon works in the case where one of the two fermions (photons) traverses a lo- 
cal homogeneous magnetic field (optically active medium), losing its spin correlations 
with the other particle. 

• In Chapter [8] I develop a mathematical method for computing analytic approxi- 
mations of the Schmidt modes of a bipartite amplitude with continuous variables. I 
maintain the analytical character of the amplitude by using complete sets of orthonor- 
mal functions. I give criteria for the convergence control and analyze the efficiency of 
the method comparing it with previous results in the literature related to entangle- 
ment of biphotons via parametric down-conversion. 

• In Chapter [9] I apply our method to a relevant case: the entanglement of two pho- 
tons created by parametric down-conversion. I compare our results (leading to well 
known, continuous functions) with those computed by standard numerical methods 
that produce sets of points: discrete functions. Both procedures agree remarkably 
well. 

• In Chapter 1101 1 consider the final products of an unstable system like an excited 
atom that emits a photon and decays to the ground state, or a nucleus that radiates 
a particle entangled with it. I analyze the momentum entanglement of these final 
particles. I study its dependence on the evolution time t and on the decay width 
r. I observe that the entanglement grows with time, until it reaches an asymptotic 
maximum, while the wider the F, the lesser the entanglement. I also compute the 
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power-law corrections in t to the exponential decay, and obtain the entangled energy 
dependence of these corrections. 

• In Chapter 1111 1 show that the decomposition of the unity in L?'{R) is in fact the 
Schmidt decomposition of the Dirac delta. It has maximum (infinite) entanglement, 
well-known result that is very easily verified from this point of view. 

• In Chapter 1121 1 analyze the operations associated to the quantum cloning sequen- 
tially implemented. I show that surprisingly the resources (Hilbert space dimension 
D) of the ancilla grow just linearly with the number of clones M to obtain. Specifi- 
cally, for universal symmetric quantum cloning I obtain D = 2M and for symmetric 
phase covariant quantum cloning, D = M + 1. Moreover, I obtain for both cases the 
isometrics for the qubit-ancilla interaction in each step of the sequential procedure. 
This proposal is easily generalizable to every quantum cloning protocol, and is very 
relevant from the experimental point of view: three-body interactions are very diffi- 
cult to implement in the laboratory, so it is fundamental to reduce the protocols to 
sequential operations, which are mainly two-body interactions. 

• In Chapter 1131 1 propose an inductive procedure to classify A^-partite entanglement 
under stochastic local operations and classical communication (SLOCC) when the 
classification for — 1 qubits is supposed to be known. The method relies in the 
analysis of the coefhcients matrix of the state in an arbitrary product basis. I illustrate 
this method in detail with the well-known bi- and tripartite cases, obtaining as a by- 
product a systematic criterion to establish the entanglement class of a pure state 
without using entanglement measures, in opposition to what has been done up to 
now. The general case is proved by induction, allowing us to obtain un upper bound 
for the number of entanglement classes of A^-partite entanglement in terms of the 
number of classes for — 1 qubits. I also include the complete classification for the 
A^ = 4 case. 

• In Appendix [A] I review the Schmidt procedure for expressing a general bipartite 
pure state as 'diagonal sum of biorthogonal products'. I describe the finite dimensional 
case and the continuous case. 

• In Appendix |B] I review the no-cloning theorem of quantum mechanics, and also 
some examples of optimal approximate quantum cloning (to a certain fidelity): sym- 
metric universal quantum cloning and symmetric economical phase-covariant quantum 
cloning. 
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In Appendix [O I review the protocol |Vid03| for expressing a multiqubit pure state 
in its matrix-product form (c/. lEckOSllPnVWCOeQ . 



Part I 



Entanglement and Relativistic 
Quantum Theory 



Chapter 2 



Dynamics of momentum 
entanglement in lowest-order QED 



In the last few years two apparently different fields, entanglement and relativity, have experi- 
enced intense research in an effort for treating them in a common framework |Cza971 IPST02t 
IAM021 IGIM IGBA031 lPS03l ITTOI |AjLMH03t iPTOil IMY041 HfeFOSl iFSMOSl IAFSMT061 
IHarOBbl IHarnOal IHWnOl I.TSSOGbl I.TSSOGaj . Most of those works investigated the Lorentz co- 
variance of entanglement through purely kinematic considerations, and only a few of them 
studied ab initio the entanglement dynamics. For example, in the context of Quantum 
Electrodynamics (QED), Pachos and Solano |PS03| considered the generation and degree 
of entanglement of spin correlations in the scattering process of a pair of massive spin-i 
charged particles, for an initially pure product state, in the low-energy limit and to the 
lowest order in QED. Manoukian and Yongram [MY04| computed the effect of spin polar- 
ization on correlations in a similar model, but also for the case of two photons created after 
e^e~ annihilation, analyzing the violation of Bell's inequality |Bel64] . In an earlier work, 
Grobe et al. |GRE94j studied, in the nonrelativistic limit, the dynamics of entanglement 
in position/momentum of two electrons which interact with each other and with a nucleus 
via a smoothed Coulomb potential. They found that the associated quantum correlations 
manifest a tendency to increase as a function of the interaction time. 

In this chapter, we study to the lowest order in QED the interaction of a pair of identi- 
cal, charged, massive spin-^ particles, and how this interaction increases the entanglement 
in the particle momenta as a function of time |LLS06b] . We chose to work at lowest order, 
where entanglement already appears full-fledged, precisely for its simplicity. In particular 
this allows to set-aside neatly other intricacies of QED, whose influence on entanglement 
should be subject of separate analysis. In this case, the generation of entanglement is a 
consequence of a conservation law: the total relativistic four-momentum is preserved in the 
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system evolution. This kind of entanglement generation will occur in any interaction verify- 
ing this conservation law, like is the case for closed multipartite systems, while allowing the 
change in the individual momentum of each component. The infinite spacetime intervals 
involved in the S-matrix result in the generation of an infinite amount of entanglement for 
interacting particles with well-defined momentum. This apparent paradox is surpassed by 
considering finite-width momentum distributions. However, it is remarkable that the at- 
tainable entanglement is not bounded from above, as we will show here. We will also discuss 
two different possibilities, with dynamical operations or with Lorentz boosts, of establishing 
transfer of entanglement between the momentum and spin degrees of freedom in the collec- 
tive two-particle system. In Section 12.11 we analyze at lowest order and at finite time the 
generation of momentum entanglement between two electrons. In Section [2.21 we apply the 
method developed in Chapter [8] (see also Refs. |LL05al ILamOSj ) to calculate the Schmidt 
decomposition of the amplitude of a pair of spin-i particles, showing the growth of momen- 
tum entanglement as they interact via QED. We obtain also analytic approximations of the 
Schmidt modes (jS.Sp and (j8.9p both in momentum and configuration spaces. In Section [2. 3^ 
we address the possibilities of transferring entanglement between momenta and spins via 
dynamical action, with Local Operations and Classical Communication (LOCC), using the 
majorization criterion |Nie99] . or via kinematical action, with Lorentz transformations. 



2.1 Two-electron Green function in perturbation theory 

To address the properties of entanglement of a two electron system one needs the amplitude 
(wave function) ^p{xl,X2) of the system, an object with 16 spinor components dependent on 
the configuration space variables xi, X2 of both particles. The wave functions were studied 
perturbatively by Bethe and Salpeter |SB51j and their evolution equation was also given 
by Gell-Mann and Low |GML5lj . The wave function development is closely related to the 
two particle Green function, 

ir(l,2;3,4) = {^o,T[il^{xi)^{x2)i^{x^)i^{xi)]-^o) (2.1) 

which describes (in the Heisenberg picture) the symmetrized probability amplitude for one 
electron to proceed from the event X3 to the event xi while the other proceeds from X4 to 
X2- If Ups(3) describes the electron at 3 and tip's' (4) that at 4, then 

^{xi,X2) = y"da^(3)da,(4)i^(l,2;3,4)7g)7(^4)nps(3)np's'(4), (2.2) 

will be their correlated amplitude at 1, 2, where 7(a) denotes the Dirac matrix associated 
to vertex a, and d(T^{a) is the differential element lying in the hypersurface orthonormal 
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Figure 2.1: Feynman diagrams for the QED interaction between two electrons (second 
order) . The minus sign denotes the antisymmetry of the amplitude associated to the fermion 
statistics. 

to the /i coordinate. In the free case this is just Ups{l) Upis'{2), but the interaction will 
produce a reshuffling of momenta and spins that may lead to entanglement. The two body 
Green function K is precisely what we need for analysing the dynamical generation of 
entanglement between both electrons. 
Perturbatively [SB51] . 

i^(l,2;3,4) = 5f(1,3)Sf(2,4) - j d^x^d^xaSp (1,5) 

X Sf{2, 6) Tf^) DpiB, 6) 7(6)^ Sf{5, 3) Sf{6, 4) 

+ ••• -{1^2} (2.3) 

where Spia, b) is the free propagator of an electron that evolves from b to a, and Dp{a, b) 
is the free photon propagator for evolution between b and a. We may call to the 

successive terms on the right hand side of this expression. They will describe the transfer of 
properties between both particles due to the interaction. This reshuffling vanishes at lowest 
order, which gives just free propagation forward in time: 

I d^xid'x2d^X3d'x^ul^,^ (1) nt,^,^(2)i^W(i, 2; 3, A^'^^^Up^sM 7(4) ^p^^, (4) 

= e{tl-t3) d{t2 - U)5s,s^ <5,,,,5(3) (pi - p,) 5(3) (p2 - p5) (2.4) 

where Ups{x) = (27r)~3/^(m/£')^/^ exp(— ipa;)ns(p). The first effects of the interaction 
appear when putting K^'^^ instead of K^^^ in the left hand side of the above equation. The 
corresponding process is shown in Fig. 12.11 To deal with this case we choose ii = t2 = 
t, ts = t4 = —t and introduce the new variables and t_ given by 
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-{t5 + tQ), t+e{-t,t), 



1 



ih-h), t.£i-it-\t+\),t-\t+\), 



in Eq. (|2.3p . which gives 



(2.5) 
(2.6) 



k^^\l,2-a,b;t) 



JlaJ/^2b 



(27r)4 ^2Ei2E22Ea2Eb 



X 



S^^HpI + P2 - Pa - Pb) 



-oo (fe°)^-(Pa-Pl)2 + ie J-i 

X / ^t_^-i{E,-E,+E,-E.+2kO)t. _ {1 ^ 2} 

it~\t+\) 



(2.7) 



whereE:(2)(l,2;a,6;t) is a shorthand notation for what corresponds to (j2.4p at second order, 
and = Ukj^ui. After some straightforward calculations, we obtain 

i?(^)(l,2;a,M) = ;i'^^^Hpi + P2-Pa-p.) 
^ ^ 47r3 ^2Ei2E22Ea2Eb 



with. 



St{t) 
Tt(t) 



sin [(^1 +E2-Ea- Ei,)t] 



E1—E2+E1, — Eg 

2 



1 



Pa - Pll 
1 



- (Pa - Pl)2 
1 1 



/i(S2-^2) i,(S2-zy2 



El + E2 — Ea — Eh 
S sin(Sf) 



S2 



+ 



E2 -1/2 



cos(St) + 



1 



S2 



,e-^^* + 



1 



-iut 



S2 



(2.8) 

(2.9) 
(2.10) 



fl = A + 2|pa-pi|, 

= -A + 2|pa-pi|, 
A = El — E2 + Eb — Ea, 



(2.11) 
(2.12) 
(2.13) 
(2.14) 
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and 



Su{t)^Stit),Tuit)^Tt{t) 



1 ^ 2 



(2.15) 



St^u are the only contributions that remain asymptotically {t — > oo) leading to the standard 
scattering amplitude, while Tt^u vanish in this limit. We recall that these are weak limits: 
no matter how large its modulus, the expression in Eq. (|2.1U|) will vanish weakly due to its 
fast oscillatory behavior. On the other hand, the sine function in Eq. (j2.9p enforces energy 
conservation 



This limit shows also that the entanglement in energies increases with time |LL05aj , reaching 
its maximum (infinite) value when t — > oo, for particles with initial well-defined momentum 
and energy. This result is independent of the chosen scattering configuration. Exact con- 
servation of energy at large times, united to a sharp momentum distribution of the initial 
states, would naturally result into a very high degree of entanglement. The better defined 
the initial momentum of each electron, the larger the asymptotic entanglement. The phys- 
ical explanation to this unbounded growth is the following: The particles with well defined 
momentum (unphysical states) are spread over all space, and thus their interaction is ubiq- 
uitous, with the consequent unbounded degree of generated entanglement. This is valid for 
every experimental setup, except in those pathological cases where the amplitude cancels 
out, due to some symmetry. In the following section, and for illustrative purposes, we will 
single out these two possibilities. 

1. The case of an unbounded degree of attainable entanglement, due to an incident 
electron with well defined momentum. We consider, with no loss of generality, a fuzzy 
distribution in momentum of the second initial electron, for simplicity purposes. 

2. Basically the same setup as [T]) but with a specific spin configuration, which leads to 
cancellation of the amplitude at large times due to the symmetry, and thus to no 
asymptotic entanglement generation. 

On the other hand, for finite times, nothing prevents a sizeable contribution from Eq. (|2.10p . 
In fact, in the limiting case where is large compared to the energies relevant in the 
problem, it may give the dominant contribution to entanglement. Whether the contribution 
from Tf(t) and T„(i) is relevant, or not, depends on the particular case considered. 




(2.16) 
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2.2 Two-electron entanglement generation at lowest order 

The electrons at x^, X4 will be generically described by an amplitude F 

V'f(x3,X4)=^ (fpa (fpbF{pa,Sa-,Pb,Sb)Up^^Sa{x3)Upb,Stix4:) (2.17) 

that should be normalizable to allow for a physical interpretation, i.e., 

[ d'^Pa j <i^Ph\F{Pa, Sa, Pb, Sb)\^ = I- (2.18) 

For separable states where F{a;b) = fa{pa, Sa)fb{Pb, Sb), fa and could be Gaussian 
amplitudes g centered around a certain fixed momentum p*^ and a certain spin component s^, 



which in the limit of vanishing widths give the standard -well defined- momentum state 
<5..o5(3)(p_pO). 

In the absence of interactions, a separable initial state will continue to be separable 
forever. However, interactions destroy this simple picture due to the effect of the correlations 
they induce. Clearly, the final state 

= E/ d^PaJd'pbK^^Hl,2;a,b;t)Fipa,Sa;Pb,Sb) (2.19) 

can not be factorized. 

In the rest of this section we analyze the final state F^^\pi, si;p2, S2',t) in Eq. (j2.19p 
to show how the variables pi and p2 get entangled by the interaction. We consider the 
nonrelativistic regime in which all intervening momenta and widths p,a <^ m, so the 
characteristic times t under consideration are appreciable. We single out the particular case 
of a projectile fermion a scattered off a fuzzy target fermion b centered around Pf| = 0. As 
a further simplification, we consider the projectile momentum sharply distributed around 
Pa i^a ^ Pa) that the initial state can be approximated by 

F(a;6) «5,^,o<5(3)(p.-p°)— ;4^e-(P-P°)V<.g. (2.20) 

Our kinematical configuration would acquire complete generality should we introduce a 
finite momenta p^ for the initial electron b. The reference system would be in this case 
midway between the lab. system and the c.o.m. system. In short, the choice p^ = will 
not affect the qualitative properties of entanglement generation. 
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Figure 2.2: Experimental setup considered in the calculations. 

We will work in the lab frame, where particle b shows a fuzzy momentum distribution 
around = 0, and focus in the kinematical situation in which the final state momenta 
satisfy pi •p2 = and also Pa-Pa = ^/V^PaPa^ a = 1,2 (see Fig. 12. 2p . This choice not only 
avoids forward scattering divergencies but also simplifies the expression of the amplitude in 
Eq. (I2.19|] . due to the chosen angles. For sure, the qualitative conclusions would also hold 
in other frames, like the center-of-mass one. We obtain 



El + E2 — Ea 



Pa— Pa 







Pb=Pl+P2 Pa 



-^{pi+p, - V2pl) + Oiip'jmjy,), 
\/2m 



(Pl+P2-Pa) _ {Pl-P\ 



2? , {p2-pVV2f 



(Pl - Pa) 
(P2 - Pa) 



(P2 



pI 
pI 



2)' + (P°)V2, 
2f + (pO)V2. 



(2.21) 



Here, boldface characters represent trivectors, otherwise they represent their associated 
norms. We perform now the following change of variables. 



turning the amplitude in Eq. (j2.19p into 
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, n 




sm[(p + q)t\ 
t 

1 



p2 ^ I eV'^ 



q2 ^ I El 



cos[(p + q)i] - e pS"' 



(2.23) 



where = ^{p + q), jl = \/2^ p'^ + ; and i = ^^t. In the following, we analyze 

different specific spin configurations in the non-relativistic limit with the help of Eq. (|2.23p . 
We consider an incident particle energy of around 1 eV<C m {p^ = 1 KeV) , and a momentum 
spreading a one order of magnitude less than p^ . We make this choice of p^ and a to obtain 
longer interaction times, of femtoseconds {t = ^r^i). Thus the parameter values we consider 



in the subsequent analysis are Pa/i^^ 
spin state for particles a and b as 



0.002 and a/m = 0.0002. We consider the initial 



ITi), 



(2.24) 



along an arbitrary direction that will serve to measure spin components in all the calculation. 
The physical results we are interested in do not depend on this choice of direction. The 
QED interaction, in the non-relativistic regime considered, at lowest order, is a Coulomb 
interaction that does not change the spins of the fermions. In fact, {jig:j^i2b) — ^n^'^^s^s^ 



S2' 



{ji^jfi2a) — 4m-^5g0gj(5s052- Given the initial spin states of Eq. ()2.24p . depending on whether 
the channel is t or n, the possible final spin states are 



\siS2)t = ITi), 
\siS2)u = lit)- 



(2.25) 
(2.26) 



Due to the fact that the considered fermions are identical, the resulting amplitude after ap- 
plying the Schmidt procedure is a superposition of Slater determinants |SCK"'"0l| IESBL02| 
IGM04| . Whenever this decomposition contains just one Slater determinant (Slater number 
equal to 1) the state is not entangled: its correlations are just due to the statistics and are 
not useful for the applications because they do not contain any additional physical informa- 
tion. If the amplitude contains more than one determinant, the state is entangled. Splitting 
the amplitude in the corresponding ones for the t and u channels, we have 
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F(^)(,,T;,,i;0. oc !M^^^e-^/^e-V2 (2.27) 

+ ^ I tsm[{p + q)i] - ^ ( cos[(p + g)t] -e"'^^*He-P'/V«'/2 

with 

p^q. (2.28) 

In the infinite time hmit the sine function converges to S{p + q), which is a distribution 
with infinite entanglement |LL05aj . The presence of the sine function is due to the finite 
time interval of integration in Eq. (12. 7p . This kind of behavior can be interpreted as a time 
difi'raction phenomenon [Mos52j . It has direct analogy with the diff'raction of electromag- 
netic waves that go through a single slit of width 2L comparable to the wavelength A. The 
analogy is complete if one identifies t with L and p + q with 2ti/\. 

In Fig. 12.31 we plot the modulus of Eq. (j2.27p versus p, q, at times i = 1,2,3,4. This 
graphic shows the progressive clustering of the amplitude around the curve q = —p, due 
to the function This is a clear signal of the growth in time of the momentum 

entanglement. Fig. 12.31 puts also in evidence the previously mentioned time diffraction 
pattern. 

We have applied the method for obtaining the Schmidt decomposition given in Ref. [LLOSaj 
(see Chapter [8l or a more complete description in Ref. [LamOSj ) to Eq. (j2.27p , considering 
for the orthonormal functions {0^^\p)}, {0^'^\q)} Hermite polynomials with their weights, 
to take advantage of the two Gaussian functions. We obtain the Schmidt decomposition for 
t = 1, 2, 3,4, where the error with matrices Cmn 12 x 12 or smaller is dmo.no — ' 10~^ in all 
considered cases. We plot in Fig. 12.41 the coefficients A„ of the Schmidt decomposition of 
Eq. (|2.27p as a function of n, for times t = 1,2, 3, 4. The number of A„ different from zero 
increases as time is elapsed, and thus the entanglement grows. 

The complete Schmidt decomposition, including channels t and u, is given in terms of 
Slater determinants jSCK"'"Ol] . and is usually called Slater decomposition. It is obtained 
antisymmetrizing the amplitude for channel t 

>{p,si;q,S2;t) oc 2^ y A„(t) y= , 

n V 2 

(2.29) 
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Figure 2.3: \F^'^\p,];q, l;t)t\ versus p, q at i = 1,2,3,4. 

where tiie modes Tpn\k,i) and Tpn\k,i) are the Schmidt modes of the channel t obtained 
for particles 1 and 2 respectively, and they correspond to the modes of the channel u for 
particles 2 and 1 respectively. 

A measure of the entanglement of a pure bipartite state of the form of Eq. (j2.29p . 
equivalent to the entropy of entanglement S, is given by the Slater number |GRE94j 




(2.30) 



K gives the number of effective Slater determinants which appear in a certain pure bipartite 
state in the form of Eq. (12.290 . The larger the value of K, the larger the entanglement. 
For K = 1 (one Slater determinant) there is no entanglement. This measure is obtained 
considering the average probability, which is given by Yl'^=o^n iYl'^=o = 1, and thus 
{A„} can be seen as a probability distribution). The inverse of the average probability is 
the Slater number. Its attractive properties are that it is independent of the representation 
of the wavefunction, it is gauge invariant, and it reaches its minimum value of 1 for the 
separable state (single Slater determinant). In Fig. 12. 5^ we show the Slater number K as a 
function of elapsed time t, verifying that the entanglement increases as the system evolves. 
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n n 



Figure 2.4: Eigenvalues A„ versus n at times t = 1,2,3,4. 



It can be appreciated in this figure the monotonic growth of entanglement, due to the 
fact that we have considered an incident electron with well defined momentum. In realistic 
physical situations with wave packets, this growth would stop, due to the momentum spread 
of the initial electrons. The general trend is that the higher the precision in the incident 
electron momentum, the larger the resulting asymptotic entanglement. The fact that the 
entanglement in momenta between the two fermions increases with time is a consequence 
of the interaction between them. We remark that the entanglement cannot grow unless the 
two particles "feel" each other. The correlations in momenta are not specific of QED: the 
effect of any interaction producing momentum exchange while conserving total momentum 
will translate into momentum correlations. 

The Schmidt modes in momenta space for the amplitude of Eq. (j2.27p are given by 

no 

^P^^\k,i) ^ e-^'/2^(V^2"n!)-V2^W(£)^„(fc) 

n=0 

a = 1,2, (2.31) 

where no is the corresponding cut-off and the values of the coefficients Amn{i) are obtained 
through the method given in Ref. |LL05aj . The modes in momenta space depend on time 
because they are not stationary states: the QED dynamics between the two fermions and 
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Figure 2.5: Slater number K as a function of the elapsed time t. 



the indeterminacy on the energy at early stages of the interaction give this dependence. By 
construction, the coefficients Amn{i) do not depend on p, q. 

We plot in Fig. 12.61 the Schmidt modes iljn^ {p, i) at times i = 1, 2, 3, 4 for n = 0,1,2,3 
(we are plotting specifically the real part of each mode only, which approximates well the 
whole mode, because Eq. ()2.27p is almost real for the cases considered). The sharper modes 
for each n correspond to the later times. Each Schmidt mode is well approximated at early 
times by the corresponding Hermite orthonormal function, and afterwards it sharpens and 
deviates from that function: it gets corrections from higher order polynomials. The fact 
that the modes get thinner with time is related to the behavior of Eq. (j2.27p at large times. 
In particular the sine function goes to d{p + q) and thus the amplitude gets sharper. 

Now we consider the Schmidt modes in configuration space. To obtain them, we just 
Fourier transform the modes of Eq. ()2.3ip with respect to the momenta pi, p2 



where a = 1,2. The dependence of p on pi and q on p2 is given through Eq. (j2.22p . The 
factor e~^2^^ in Eq. ()2.32p is the one which commutes the states between the interaction 
picture (considered in Eq. (|2.3ip and in the previous calculations in Sees. 12.11 and 12. 2p and 
the Schrodinger picture. 

The Hermite polynomials obey the following expression [GR48j 




(2.32) 




(2.33) 
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Figure 2.6: Schmidt modes ipn\p,i) at times i = 1,2,3,4 for n = 0,1,2,3. The sharper 
modes for each n correspond to the later times. 



With the help of Eq. ()2.33p and by linearity of the Fourier transforms, we are able to obtain 
analytic expressions for the Schmidt modes in configuration space (to a certain accuracy, 
which depends on the cut-offs considered). This is possible because the dispersion relation 
of the massive fermions in the considered non-relativistic limit is Ea = and thus the 



integral of Eq. (j2.32|) can be obtained analytically using Eq. (j2.33|) . 

The corresponding Schmidt modes in configuration space are then given by 



no 



n=0 



where the orthonormal functions in configuration space are 

-in arctan(o't)+j5-~-'- f/2) 



-e 2(i+»5i)iJ„ 



(2.34) 



,(2.35) 



a 



1,2, a 



7Pi 



^t. The modes in Eqs. (12:341) and (fOKj) are 



In Eqs. ()2.34p and (j2.35p . we are using dimensionless variables, x 

and the dimensionless time defined before, t 
normalized in the variables Xa- The orthonormal functions of Eq. ()2.35p propagate in space 
at a speed and they spread in their evolution. Additionally, the modes of Eq. (j2.34p 

have also the time dependence of Amn{i). The Slater decomposition in configuration space. 
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obtained Fourier transforming the modes of Eq. (j2.3ip is 



n ' 

(2.36) 

The coefficients A„(t) are unaffected by the Fourier transformation, and thus the degree of 
entanglement in configuration space is the same as in momenta space. 
We consider now the initial spin configuration 

Is'ast) = ITT), (2.37) 
where, the only possible final state in the non-relativistic limit is 

\s1s2) = ITT)- (2.38) 

In this case, the sine term goes to zero, because the momentum part of this term is antisym- 
metric in p^, and the sine function goes io 5{p + q), which has support (as a distribution) 
on g = —p. We point out that the sine contribution to this amplitude is negligible because of 
the particular setup chosen. In other experiment configurations the amplitude in Eq. (j2.19p 
associated to the spin states of Eqs. ()2.37p and (|2.38p would have appreciable sine term 
and thus increasing momenta entanglement with time. On the other hand, in this case the 
contribution from Tt{t) in Eq. (|2.10p and T„(t) is even smaller than the sine term, and 
converges weakly to zero. 

We plot in Fig. 12.71 the real and imaginary parts of the term associated to T((t) and 
Tu{t) in Eq. (j2.23p . which we denote by g{p,q,t), for spin states of Eqs. (j2.37p and (j2.38p 
as a function of time t S (1, 1.001) and having p = 1, q = 1.2. We want to show with it 
the strong oscillatory character of the amplitude with time, and how all the contributions 
interfere destructively with each other giving a zero final value. This is similar to the 
stationary phase procedure, in which only the contributions in proximity to the stationary 
value of the phase do interfere constructively and are appreciable. What we display here is 
the weak convergence to zero for the functions Tt{t) and Tu{t). 

In this section, we investigated the generation of entanglement in momenta between 
two identical spin-^ particles which interact via QED. We showed how the correlations 
grow as the energy conservation is increasingly fulfilled with time. The previous calculation 
had, however, the approximation of considering a projectile particle with perfectly well 
defined momentum, something not achievable in practice. This is a first step towards a 
real experiment, where both fermions will have a dispersion in momenta and thus infinite 
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Figure 2.7: Real and imaginary parts of the amplitude g{p,q,i) for p = 1, q = 1.2, and 
t € (1, 1.001) (arbitrary units). 

entanglement will never be reached, due to the additional integrals of the Dirac delta 6{AE) 
over the momentum spread. We believe these results have much interest both from the 
fundamental point of view but also from the experimental one, for example in fermion gases 
ISALW04] . 

2.3 Entanglement transfer between momentum and spin 
2.3.1 Dynamical transfer and distillation 

In the previous section we computed the entanglement in momenta for a pair of electrons 
which interact through exchange of a virtual photon. The results could be summarized by: 
the more collimated in momentum is the incident electron, and the more time is elapsed, the 
larger the entanglement in momenta obtained. Heisenberg's principle, on the other hand, 
establishes a limit to the precision with which the momentum may be collimated and hence 
to the achievable degree of entanglement. 

It is possible in principle to transform the entanglement in momenta into entanglement 
in spins. This is easily seen in terms of the majorization criterion |Nie99t iNCOOj . This is of 
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practical interest because the experimenter usually manipulates spins, thus spin entangle- 
ment seems to be more useful. Here we analyze this entanglement transfer. 

Majorization is an area of mathematics which predates quantum mechanics. Quoting 
Nielsen and Chuang, "Majorization is an ordering on d-dimensional real vectors intended to 
capture the notion that one vector is more or less disordered than another" . We consider a 
pair of d-dimensional vectors, x = {xi, ...,Xd) and y = (yi, ...,yd)- We say x is majorized by 
y, written x ^ y, if Yl^=i — X]^=i Vj k = 1, d, with equality instead of inequality 
for k = d. We denote by the components of z in decreasing order (z| > Z2 > ... > zj^). 
The interest of this work in the majorization concept comes from a theorem which states 
that a bipartite pure state may be transformed to another pure state \(p) by Local 
Operations and Classical Communication (LOCC) if and only if -< A,^, where A^, A^ are 
the vectors of (square) coefficients of the Schmidt decomposition (jA.lip (finite dimension 
case) or ()A.14p (continuous case) of the states {(p), respectively. LOCC adds to those 
quantum operations effected only locally the possibility of classical communication between 
spatially separated parts of the system. According to this criterion, it would be possible in 
principle to obtain a singlet spin state beginning with a momentum entangled state \ip) 
whenever A^ ^ A,;!,. 

The possibility of obtaining a singlet spin state from a momentum-entangled state can 
be extended to a more efficient situation: the possibility of distillation of entanglement. 
This idea consists on obtaining multiple singlet states beginning with several copies of a 
given pure state The distillable entanglement of {ip) consists in the ratio n/m, where 
m is the number of copies of {tp) we have initially, and n the number of singlet states we 
are able to obtain via LOCC acting on these copies. It can be shown [NCO O] that for pure 
states the distillable entanglement equals the entropy of entanglement, S (jll.2p . Thus, in 
the continuous case (infinite-dimensional Hilbert space), the distillable entanglement is not 
bounded from above, because neither is S. According to this, the larger the entanglement 
in momenta the more singlet states could be obtained with LOCC. 

To illustrate the possibility of entanglement transfer with a specific example, we consider 
a momentum-entangled state for two distinguishable fermions 



This state has associated a vector A^ = (1/2, 1/2,0,0, ...). On the other hand, the singlet 
state 




(2.39) 




(2.40) 



has associated a vector A 



(1/2,1/2,0,0...). 
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These vectors obey -< and thus the state entangled in momenta may be trans- 
formed into the state entangled in spins via LOCC. The operations needed to achieve this 
are in general unitary operations acting on the degrees of freedom of momentum and spin 
of each individual electron, local projective measurements, and classical communication 
between the separated parts of the system (the experimental setups located where the 
electrons reach). The unitary local operations may be implemented by, for example (in 
the non-relativistic case), electric fields (which modify the momentum) and magnetic fields 
(which modify the spin and the direction of the momentum) combined to give the desired 
effect. The specific setup needed to do this might be rather complicated. In addition, 
decoherence effects should be avoided. What we point out is that majorization states this 
transformation would be possible in principle. 

2.3.2 Kinematical transfer and Lorentz boosts 

Another approach to the study of entanglement transfer between momentum and spin 
degrees of freedom is the kinematical one. In fact, the Lorentz transformations may entangle 
the spin and momentum degrees of freedom. To be more explicit, and following the notation 
of Ref. |GA02j . we consider a certain bipartite pure wave function gxai'P^'^) for two spin-i 
fermions, where A and a denote respectively the spin degrees of freedom of each of the two 
fermions, and p and q the corresponding momenta. This would appear to an observer in a 
Lorentz transformed frame as 

5A<x(p,q) ^ J]C/itr'PV^;r'^^5AV'(A-'p,A-'q), (2.41) 

^ A' a-' 

where 

^Iv ^=<''^(^(A,P)) (2.42) 
is the spin ^ representation of the Wigner rotation i?(A, p). The Wigner rotations of 
Eq. (I2.42|] can be seen as conditional logical operators, which rotate the spin a certain angle 
depending on the value of the momentum. Thus, a Lorentz transformation will modify 
in general the entanglement between momentum and spin of each individual electron. We 
distinguish the following particular cases. 
Product state in all variables. In this case, 

5Aa(p,q) =5i(p)52(q)|A)|a), (2.43) 

and the entanglement at the rest reference frame is zero. Under a boost, the Wigner 
rotations of Eq. (j2.42p entangle the momentum of each fermion with its spin, and thus the 
entanglement momentum-spin grows |PST02| . 
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Entangled state spin-spin and/or momentum-momentum. We consider now a state 

<?A.(p,q) =/(p,q)|</') (2.44) 

witli \<j)) an arbitrary state of the spins, and /(p, q) an arbitrary state of tlie momenta. In 
this case, a Lorentz boost will entangle in general each spin with its corresponding momen- 
tum, and a careful analysis shows that the spin-spin entanglement never grows [GA02j . Of 
course, by applying the reversed boost the entanglement momentum-spin would be trans- 
ferred back to the spin-spin one, and thus the latter would grow. This particular case shows 
that, for the state we considered in Sec. 12.21 given by Eqs. (|2.24p . (j2.25p and (j2.27p . the en- 
tanglement could not be transferred from momenta into spins via Lorentz transformations. 
Thus, the dynamical approach would be here more suitable. 

Entangled state momentum-spin. According to the previous theorem, the momentum- 
spin entanglement may be lowered, transferring part of the correlations to the spins, or 
increased, taking some part of the correlations from them. To my knowledge, there is not 
a similar result for the momentum, that is, whether the momentum entanglement can be 
preserved under boosts, or it suffers decoherence similarly to the spins, and part of it is 
transferred to the momentum-spin part. 
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Generation of spin entanglement 
via spin-independent scattering 

3.1 Spin-spin entanglement via spin-independent scattering 

A compound system is entangled when it is impossible to attribute a complete set of prop- 
erties to any of its parts. In this case, and for pure states, it is impossible to factor the 
state in a product of independent factors belonging to its parts. In this chapter we will 
consider bipartite systems composed of two s = ^ fermions. Our aim is to uncover |LL06j 
some specific features that apply when both particles are identical. They appear itemized 
three pages below. 

States of two identical fermions have to obey the symmetrization postulate. This implies 
that they decompose into linear combinations of Slater determinants (SLs) of individual 
states. Naively, as these SLs cannot be factorized further, indistinguishability seems to 
imply entanglement. This is reinforced by the observation that the entropy of entanglement 
(EoE) is bounded from below by 5 > 1, well above the lower limit S = for a pair of 
non-entangled distinguishable particles. So, it looks like there is an inescapable amount of 
uncertainty, and hence of entanglement, in any state of two identical fermions. 

The above issue has been extensively examined in the literature [SCK"'"OH IESBL02"| 
IGM04| with the following result: Part of the uncertainty (giving S = 1) corresponds to the 
impossibility to individuate which one is the first or the second particle of the system. This 
explains why the lower limit for the EoE is 1. Consider for instance two identical s = 
fermions in a singlet state 

Xs :=4s[x(l)^x(2)^-x(l)^x(2)^]. 
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The antisymmetrization does not preclude the assignment of properties to the particles, but 
only assigning them precisely to particle 1 or particle 2. The reduced density matrix of any 
of the particles is p = with an EoE S{p) = 1. 

The portion of S above 1 (if any) is genuine entanglement as it corresponds to the 
impossibility of attributing precise properties to the particles of the system [GM04] . Assume 
for instance that we endow the previous fermions with the capability of being outside 
{x = "0) or inside (x = 92) the laboratory ((^/;*, 'i/''') = '5*"') = ^^'^ 1 ii^^-,^^) = Oj 

=T;i)- We now have two different possibilities: either the fermion outside has spin up 
or spin down {il^^). Hence, there are two different SLs for a system built by a pair of 
particles with opposite spins, one outside, the other inside the laboratory 

SL(l,2)i = -^[V^(l)T^(2)i-^(l)i^(2)T], 

SL(1,2)2 = -^[V^(1)V(2)^-V^(1)^^(2)^]. (3.1) 

They form two different biorthogonal states, the combination [SL(l,2)i — SL(1, 2)2]/\/2 
corresponding to the singlet and [SL(1, 2)i + SL(1, 2)2]/ \/2 to the triplet state (with respect 
to the total spin s = si + S2). An arbitrary state <I>(1, 2) would then be a linear combination 
of these two SLs: 

$(1, 2) = ciSL(l, 2)1 + C2SL(1, 2)2, led' = 1> (3-2) 

i 

giving an EoE 

S = log2|Qp > 1. (3.3) 

i 

Clearly, when ci or C2 vanish, we come back to S" = 1, as the only uncertainty left is 
the very identity of the particles. Summarizing, while indistinguishability is an issue to be 
solved by antisymmetrization within each SL, entanglement is an issue pertaining to the 
superposition of different SLs [SCK+Oll lESBlM IGM04] . At the end, we could even decide 
to call 1 to the variables of the outside particle, and forget about symmetrization 

«>(1, 2) ^ ci^(l)T (/?(2)i + C2V'(1)^ (3.4) 

as both particles are far away from each other. In this case, the EoE, equal to S" = 
— |cjp log2 |cjp > is lesser than the one corresponding to antisymmetrized states by 
a quantity of 1, which is just the uncertainty associated to antisymmetrization. From 
now on we will consider the latter definition of 5, which gives the genuine amount of 
entanglement between the two particles. Notice that for half-odd s, the number t^SL of 
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Slater determinants is bounded by #SL < (2s + l)(i/2, where d is the dimension of each 
Hilbert space of the configuration or momentum degrees of freedom for each of the two 
fermions. 

Much in the same way as above, we could consider one of the particles as right moving 
ix = V'o) the other as left moving {x = V'tt)) giving rise to two SLs in parallel with the 
above discussion. This is the first step towards the inclusion of the full set of commuting 
operators for the system. In addition to the spin components (si, ,52) or helicities, there are 
the total P and relative p momenta. In the center of mass (CoM) frame we could consider 
the system described by the continuum of SLs 

SL(l,2;p), = i=[V^(l)gV'(2);^-V'(l)rV'(2)S], 

SL(l,2;p)_, = -L[V;(l)-^(2)^-V(l)^V(2)o^], (3.5) 

where V'(l)o ~ (^IP"*) ^^'^ V'(l)^ — ~ P^)- The labels and vr are the azimuthal angles 
when we laid the axes along p. Finally, there is a pair of SLs for each p, so that a general 
state made with two opposite spin particles with relative momentum p could be written in 
the form: 

$(1,2)0 ^ ^ c,(p)SL(l,2;p)„ (3.6) 

s=±l/2 

with Yls=±i/2 ks(p)P ~ 1- Again, we run into the impossibility to tell which is 1 and 
which is 2. In addition there may be some uncertainty about the total spin state, whether 
a singlet or a triplet, or conversely, about the spin component of any of the particles, ■00 or 

After this discussion it should be clear to what extent entanglement and distinguishabil- 



ity belong to different realms SCK"'"0l[lESBL02llGM04] . A requirement to include identical 



particles is to symmetrize the expressions used for unlike particles. Until now, we have only 
considered the free case. We have to examine the case of two interacting particles, as in- 
teraction is expected to be the source of subsequent entanglement |LRM761 ITor851 lALPOli 
IP5031 IMY041 IAAMS041 ISALW041 lHa?05l ITK051 iLLSOGbl IWanOSl IHar06bl IHarOGal IH W06j . 
Obviously, the answer may depend on a tricky way on the detailed form of the interaction, 
of its spin dependence in particular . It also seems that the role of particles identity, if any, 
will be played through symmetrization. 

In the following we will show that spin entanglement is generated for the case of two 
interacting spin-^ identical particles, with the following features: 

• Spin-spin entanglement is generated even by spin independent interactions. 
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• In this case, it is independent of any symmetrization procedure. 

• This phenomenon does not appear for unlike particles. 

We first tackle the scattering of two unequal s = particles A and B which run into 
each other with relative CoM momentum p. We set the frame axes by the initial momentum 
p of particle A, and let the spin components he Sa = s and Sb = —s along an arbitrary but 
fixed axis. We will consider a spin independent Hamiltonian H, so the evolution conserves 
Sa and Sfc. We denote by Aq (Bg) the state of particle A (B) that propagates along direction 
6 with spin s. In these conditions the scattering proceeds as: 

a>in = A'o = m B-%, (3.7) 

where 6 is the scattering angle and fp{6) the scattering amplitude. We will consider 9 
different from or vr to avoid forward and backward directions. While the increase of 
uncertainty due to the interaction is clear, because a continuous manifold of final directions 
with probabilities opened up from just one initial direction, spin remains untouched. 

The information about Sa is the same before and after the scattering; as much as we knew 
the initial spin of A, we know its final spin whatever the final direction is. In other words, 
spin was not entangled by the interaction. We will now translate these well known facts to 
the case of identical particles, where they do not hold true. 

Let particle B be identical to A. Consider the same initial state as before: A particle A 
with momentum p and spin s runs into another A with momentum — p and spin —s. Notice 
there is maximal information on the state. We could write <I>jj^ = AqA~'^, and eventually 
symmetrize. We now focus on the final state. It is no longer true that particle A will come 
out with momentum p' and spin s with amplitude fp{9) while the amplitude for coming 
out with momentum p' and spin —s vanishes. Recalling that B above did become A, the 
two cases fp{9)AgB~^Q and fp{-K — 9)A^_gBg'^ fuse into a unique state 

^out(^) = fpiO) At A-% + fpiTT - 9) AU A^, (3.8) 

as shown in Figure 13.11 Notice the uncertainty acquired by the spin: Now particle A 
comes out from the interaction along 9 either with spin s or with spin —s, with relative 
amplitudes fp{9) and fp{-K — 9) respectively. In other words, spin was entangled during the 
spin independent evolution. Here, it is not the spin dependence of the interaction, but the 
existence of additional degrees of freedom which generate spin-spin entanglement. These, 
act as ancillas creating an effective spin-spin interaction that entangles the two fermions. 
The ancilla and the degree of entanglement depend on the scattering angle 9. Notice that 
for 9 = 7r/2 both amplitudes fp{9) and /^(vr — 9) become equal, so that the degree of 
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Figure 3.1: Schematic picture of the two channels that contribute to the spin-independent 
scattering of two identical fermions. The shaded regions denote an arbitrary spin- 
independent interaction between the two fermions. The vertical arrows |, | indicate the 
corresponding third component of spin. 

generated entanglement is maximal, 1 ebit. On the other hand, for ^ ~ 0, it generally holds 
fp{6) ^» fp^TT — 0), so that in the forward and backward scattering almost no entanglement 
would be generated. However, this depends on the specific interaction. In Section 13.21 we 
will clarify this point with Coulomb interaction. 

Symmetrization does not change this, it only expresses that we can not tell which one 
is 1 and which one is 2. The properly symmetrized initial state is 

= SL(l,2;p), = ±[A{ir^A{2)~^ - A(l);M(2)g]. (3.9) 

The scattering process could be written in terms of SLs as 

SL(1, 2; p), fp{e) SL(1, 2; p')s - - 9) SL(1, 2; p')-., (3.10) 

where p' is the final momentum and the Slater determinants are given in (j3.5p . Both, 
this expression and Eq. (j3.8p . describe the same physical situation and lead to the same 
entanglement generation. 

The bosonic case may be analyzed in an analogous way. The modification for two- 
dimensional spin Hilbert spaces (i.e. photons) would be a sign change in Eqs. (|3.5p . ()3.9p 
and (|3.10p . as bosonic statistics has associated symmetric states. The equivalent of Eq. 
()3.10p for bosons is a genuine entangled state for 7^ 0, vr, much as in the fermionic case. 
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3.2 A specific example: Coulomb interaction at lowest order 

We now consider Coulomb interaction at lowest order to illustrate the reasonings presented 
above. In this case 



fpiO) 



m 



U^-O) = (3.11) 

where A^(e) is a numerical factor depending on the charge e. t{6) and u{9) are two of the 
Mandelstam variables, associated to t and u channels respectively, and depending on the 
scattering angle 9. for initial p and final p' relative 4- momenta of the scattering fermions, 
they are given hy t = {p — p')'^, u = {p + p')^- the CoM frame, 

t{0) := 2(m2 -^2)(1 - cos 61), 

u{e) := 2{m^ - E^){l + cose), (3.12) 

where m is the mass of each fermion and 2E is the available energy. 

According to this, the spin part of the state (j3.8p for this case, properly normalized, is 

\xe) = f+{0)\n) - f-mm, (3.13) 

being 

-^^24=. (3.14) 

y 2(1 + cos^ 0) 

The two amplitudes /+ and /_ vary monotonously as 9 grows, becoming equal for 9 = Tr/2. 
The physical meaning for this is that for — > 0, the knowledge about the system is maximal 
and the entanglement minimal (zero), and for increasing 9 the knowledge of the system 
decreases continuously until reaching its minimum value at 9 = it/2. Accordingly, the 
entanglement grows with 9 until reaching its maximum value for 9 = tt/2. 

We plot in Figure the EoE |TK05| S{9) = - f+{9f\og2 f+{9f - f^{9f\og2f-{6? 
of state p.l3p as a function of 9, for < < 7r/2. The entanglement grows monotonically 
until 9 = '7r/2, where it becomes maximal (1 ebit). 



3.3 Violation of Bell's inequality as a function of the scatter- 
ing angle 6 

In order to analyze the role the 9 scattering angle plays in the generation of these genuine 
quantum correlations, we consider now the degree of violation of Bell's inequality as a 
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Figure 3.2: EoE S{9) as a function of 6. 



function of 6. To this purpose, we define |Bel64t ICS78j the observable 

^(a,b) := -a^d^^) -6)1$) 

= -[azbz + 2fp{9)fp{7r - 9){axhx + ayhy)], 



(3.15) 



where |<I>) := \^ovit{Q)) is the (normalized) state ()3.8p and a, b are arbitrary unit vectors. 
In Eq. (I3.15P we consider the amplitudes fp{9) and /p('7r — 9) normalized for each 9, in 
the form |/p(^)P + |/p('?r — ^)P = 1- We consider three coplanar unit vectors, a, b and c. 
(a, b) = 7r/3, (a, c) = 27r/3 and (b, c) = 7r/3. We have 



|^(a,b)-E(a,c)| = 1, 

F{9) := 1 + E{h, c) = \- \a9)U{n - 9). 



The Bell's inequality, given by |Bel64l [CS78] 

|£'(a, b) - E{k, c)| < 1 + E{h, c) 

will then be 

F{9) > 1. 



(3.16) 



(3.17) 



(3.18) 



For the particular case of Coulomb interaction at lowest order here considered, 2fp{9)fp{7r- 
9) = 2/+(0)/_(0) = (1 — cos^ + cos^ 9) and thus the critical angle for which the in- 
equality becomes violated is 9c = '7r/4 for F{9c) = 1- For 9c < 9 < tt/2 the Bell's inequality 
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Figure 3.3: F{0) as a function of 9. The classical-quantum border corresponds to F{9c) = 1, 
with 9c = vr/4. 

does not hold. We show in Figure [3^ the 9 dependence of F{9) together with the classical- 
quantum border, F = 1, at 6c = vr/4. Thus, for experiments with 6 > 7r/4 one could be 
able in principle to discriminate between local realism and quantum mechanics. This is in 
contrast with recent analysis of Bell's inequality violation in elementary particle systems 
|Go04[ IBEG051 ITor86| . where the emphasis was placed on flavor entanglement, K^K^, 
B^B^, and the like. These analysis presented j BBGH04] two kinds of drawbacks coming 
from the lack of experimenter's free will and from the unitary evolution with decaying states. 
These issues reduce the significance of the experiments up to the point of preventing their 
use as tests of quantum mechanics versus local realistic theories. The spin-spin entangle- 
ment analyzed in this chapter does not have this kind of problems and could be used in 
principle for that purpose. 



3.4 The tripartite case 

Our protocol can be straightforwardly generalized to more complicated settings. To show 
this, we consider now three identical, massive, s = 1 particles which interact through some 
kind of two- and three-body interaction, which for our practical purposes is completely 
general, but spin independent. We suppose the particles evolve in a plane, and the pro- 
longations of the directions of propagation intersect at one point, having relative angles 
of 27r/3 degrees among each two neighbor directions, see Figure [331 We will label these 
in-directions with angles 0, 27r/3 and 47r/3, which will also be the final angles for illustrative 
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purposes. Our initial state is now 

'^'in = ^o^L/3^4"J/3' (3-19) 

where Ag denotes, Uke previously, the state of particle of type A with propagation direction 
9 and spin s. Our final state will now be a superposition of all the possible final states 
associated to the different channels. There are as much as permutations of the three states 
^4", A^ and A~^, which can evolve each through directions 0, 2ir/3 and An/ 3. 

For example, in case there is no interaction, the three states will evolve according to 
A|] ^ A[5, ^^^/3 ^27r/3' ^4^/3 ~^ ^4^/3- ^^^^ ^ive a Weight /o,o,o to this channel, 
where the labels mean the in- and out-angle for each particle remains the same. 

In the case where there is a two-body interaction among the particles, the possible 
channels are those which permutate two of the final legs of the Feynman diagram. There 
are three channels associated to two-body scattering: {0 — > 2tt/3, 2tt/3 — > 0, 47r/3 — > Att/S}, 
{0 ^ 47r/3,27r/3 ^ 27r/3,47r/3 ^ 0}, and {0 ^ 0,27r/3 ^ 47r/3,47r/3 ^ 27r/3}. One 
example would be the following evolution, A^ ^27r/3' ^2tt/3 ~^ ^O' ^^'^ ^4^/3 ~^ ^4^/3' 
As in previous sections, the spins remain unchanged. We will give a weight /27r/3,27r/3,o to 
the channels associated to two-body interactions, which means that two of the particles 
experience a trajectory change of 27r/3. 

Finally, in case there is a three-body interaction among the three particles, the final 
channels will be those that permutate the three legs of the Feynman diagram. The two 
available channels for three-body interactions are {0 — > 27r/3, 27r/3 — > 47r/3, 47r/3 0} 
and {0 47r/3,27r/3 0, 47r/3 27r/3}. One example would be the evolution 

^27r/3' ^27r/3 ^ ^L/3' ^^'^ ^4^/3 ^ ^O^' We wih give a weight /27r/3,27r/3,27r/3 tO the 

channels associated to three-body scattering, which means that the three particles change 
their propagation direction by 27r/3. Of course, three-body interactions are unusual, so 
typically we wih have /o,o,o > /27r/3,27r/3,o > /27r/3,27r/3,27r/3- Accordingly, the final state 
obtained this way will be 

^OUt = /0,0,0^0^27r/3^4V3 

+ /27r/3,27r/3,o[-42,r/3"^0"^4j/3 + "^47r/3^27r/3^0 ^ + ^0^47r/3^2V3] 
+/27r/3,27r/3,27r/3[^27r/3^47r/3^0^ + ^47r/3^0^2"v3] ' (^-^O) 

We plot in Figure 13.41 the six different diagrams associated to the six channels of the 
spin- independent scattering of three identical s = 1 particles. 

In the asymptotic limit, once the particles are far apart, their wave packets will not 
overlap anymore and we may consider the particles as distinguishable. We will denote by 



46 



Generation of spin entanglement via spin-independent scattering 




Figure 3.4: Schematic diagrams for the tripartite spin-independent scattering of three iden- 
tical s = 1 particles. The spins associated to the incident directions are s = (Oin = 0), 
s = 1 (^in = 27r/3), and s = —1 (^in = 47r/3) in all the six cases. 

52,53) the spin state of the three particles, where si is associated to the particle that 
propagates at the end along 0, S2 to the one that propagates along 27r/3 and S3 to the one 
that moves along 47r/3. The asymptotic spin state associated to Eq. (I3.20p will be, in this 
notation, 

Ix) = /o,o,o|0,l,-l> 

+ /2V3,2V3,0[|1, 0,-1) +1 -1,1,0) +10,-1,1)] 

+/2./3,2./3,2./3[|l, -1,0) + 1-1,0, 1)]. (3.21) 

Notice that in case we considered only two spins, 1 = and —1 (either three spin-^ 
particles or three spin-1 particles with just two spin states of the 3-dimensional s = 1 
Hilbert space), the part associated to two-body interaction would be a state, a state 
with genuine tripartite entanglement. Thus, in this way it could be possible to generate 
genuine multipartite entanglement via spin-independent scattering. A generalization to the 
A^-partite case of our procedure is straightforward from the tripartite and bipartite cases 
and we will omit it here. 

In summary, we analyzed the relation between entanglement and antisymmetrization 
for identical particles, in the context of spin-independent particle scattering. We showed 
that, in order to create genuine spin-spin quantum correlations between two s = ^ fermions, 
spin-dependent interactions are not compulsory. The identity of the particles along with 
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an interaction between degrees of freedom different from the spin, suffice for this purpose. 
The entanglement generated tliis way is not a fictitious one due to antisymmetrization, but 
a real one, and violates a certain Bell's inequality foi 6 > 9c = tt / 4:. 



Chapter 4 

Relativity of distillability 



4.1 Lorentz transformed spin density matrix 

As we have pointed out along this Thesis, nowadays entanglement is considered a basic 
resource in present and future applications of quantum information, communication, and 
technology [NCOOl IGMD02] . However, entangled states are fragile, and interactions with 
the environment destroy their coherence, thus degrading this precious resource. Fortu- 
nately, entanglement can still be recovered from a certain class of states which share the 
property of being distillable. This means that even in a decoherence scenario, entangle- 
ment can be extracted through purification processes that restore their quantum correla- 
tions [BBP+96[ IBMDOSj . An entangled state can be defined quantum state that is 
not separable, and a separable state can always be expressed as a convex sum of prod- 
uct density operators |Wer89j . In particular, a bipartite separable state can be written as 
p = Ylii^iP^f' ® Pi'^^ where Q > 0, Yli^i — 1' ™d p[^^ and pf*^ are density operators 
associated to subsystems A and B. 

In the context of quantum field theory, special relativity (SR) [RinOl] and quantum me- 
chanics are described in a unified manner. From a fundamental point of view, in addition, 
it is relevant to study the implications of SR on the modern quantum information theory 
(QIT) |PT04j . Recently, Peres et al. jPST02j have observed that the entropy of entangle- 
ment of a single spin-^ density matrix is not a relativistic invariant, given that Wigner 
rotations |Wig39| entangle the spin with its momentum distribution when observed in a 
moving referential. This astonishing result, intrinsic and unavoidable, shows that entan- 
glement theory must be reconsidered from a relativistic point of view. On the other hand, 
the fundamental implications of relativity on quantum mechanics could be stronger than 
what is commonly believed. For example, Wigner rotations induce also decoherence on two 
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entangled spins |AM02l lP503l IC;A021 lUBAOSj . However, they have not been studied yet in 
the context of mixed states and distillable entanglement [Per 961 IHHH98j . 

A typical situation in SR pertains to a couple of observers: one is stationary in an inertial 
frame S and the other is also stationary in an inertial frame S' that moves with velocity v 
with respect to S. The problems addressed in SR consider the relation between different 
measurements of physical properties, like velocities, time intervals, and space intervals, 
of objects as seen by observers in S and S' . However, in QIT, it is assumed that the 
measurements always take place in a proper reference frame, either S or S' . To see the 
effects that SR can bring about in QIT problems |PT04] . we need to enlarge the typical 
situations where quantum descriptions and measurements take place. 

In order to analyze the new possibilities that SR offers, we introduce |LMDS06] the 
following concepts 

• Weak isoentangled state p^^^: A state that is entangled in all considered reference 
frames. This property is independent of the chosen entanglement measure £. 

• Strong isoentangled state /o|^^: A state that is entangled in all considered reference 
frames, while having a constant value associated with a given entanglement measure 
S. This concept depends on the £ chosen. 

• Weak isodistillable state p^^^: A state that is distillable in all considered reference 
frames. This implies that the state is entangled for these observers. 

• Strong isodistillable state A state that is distillable in all considered reference 
frames, while having a constant value associated with a given entanglement measure 
£. This concept depends on the £ chosen. 

In general, the following hierarchy of sets holds (see Fig. 14.11 for a pictorial representa- 
tion) 

To illustrate the relative character of distillability, let us consider |LMDS06] the specific 
situation in which Alice (A) and Bob (B) share a bipartite mixed state of Werner type with 
respect to an inertial frame S. Moreover, in order to complete the SR+QIT scenario, we 
also consider another inertial frame S' , where relatives (A') and (B') of Alice (A) and Bob 
(B) are moving with relative velocity v with respect to S. Using the picture of Einstein's 
trains, we may think that (A) and (B) are at the platform of the railway station sharing 
a set of mixed states, while their relatives (A') and (B') are travelling in a train sharing 
another couple of entangled particles of the same characteristics. The mixed state is made 
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Figure 4.1: Hierarchy for tiie sets of states WIE, SIE, WID, and SID. 

up of two particles, say electrons witli mass m, having two types of degrees of freedom: 
momentum p and spin s = ^. The former is a continuous variable while the latter is a 
discrete one. By definition, we consider our logical or computational qubit to be the spin 
degree of freedom. Each particle is assumed to be locahzed, Uke in a box, and its momentum 
p wiU be described by the same Gaussian distribution. We assume that the spin degrees of 
freedom of particles A and B are decoupled from their respective momentum distributions 
and form the state 

pF ■■= F\^p){%\ + + + (4-2) 

Here, is a parameter such that < F < 1, 
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(4.3) 



where Pa and Pb are the corresponding momentum vectors of particles A and B, as seen in 
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S, and 




<^(Pb) : = 

with Gaussian momentum distribution ^(p) := vr" 
It) and 11) represent spin vectors pointing up and down along the z-axis, respectively. If 
we trace momentum degrees of freedom in Eq. (j4.3p . we obtain the usual spin Bell states, 
ll^*^), 1^^), 1*1* ), |*1*^)}- If we do the same in Eq. (|4.2p . we remain with the usual spin 
Werner state |Wer89j that can be written in its matrix form as 



3/4^-3/2 



(4.4) 



exp(— p^/2?i;^), being p := |p| 
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(4.5) 



It is known that Bell state ) is distillable out of Eq. (j4.5|) if, and only if, F > 1/2. 

We consider also another pair of similar particles, A' and B', with the same state as 
A and B, Pgi^' = Pg^, but seen in another reference frame 5'. The frame 5' moves with 
velocity v along the x-axis with respect to the frame S. 

When we want to describe the state of A' and B' as observed from frame S, rotations 
on the spin variables, conditioned to the value of the momentum of each particle, have to 
be introduced. These conditional spin rotations, considered first by Wigner |Wig39 , are 
a natural consequence of Lorentz transformations. In general, Wigner rotations entangle 
spin and momentum degrees of freedom for each particle. We want to encode quantum 
information in the two qubits determined by the spin degrees of freedom of our two spin- 
\ particles. As a consequence, the reduced spin system, after a Lorentz transformation, 
increases its entropy and reduces its initial degree of entanglement. If we consider the 
velocities of the particles as having only non-zero components in the z-axis, each state 
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vector of A' and B' in Eq. (j4.4p transforms as 



*.(p)^(,;,)^Ai*.(p)i^(-::M5(p). 



(4.6) 



where cos 6p and sin 9p express Wigner rotations conditioned to tlie value of tiie momentum 
vector. 

The most general bipartite density matrix in the rest frame for arbitrary spin-^ states 

(a) (a) 

and Gaussian product states in momentum, is spanned by the tensor products of , ^2 > 
and ^-^^^ and can be expressed as 

p= E Q,kM^\p^) ^ ^f\p^)[^t\p'J ^ '^f'HM^ ■ (4.7) 

ijkl=l,2 

Under a boost, Eq. (j4.7p will transform into 
ApAt = Q.-fc;A(")[*S")(pa)] ® AW[^&f (pb)]{A(-)[M/(")(p'J] AW[*f)(p'b)]}t. 

ijkl=l,2 

(4.8) 

Tracing out the momentum degrees of freedom, we obtain 

TVp,,p,(AM^) = E C^.,HTrp,(A(^)[^S'^(Pa)]{A(^)[^l^Va)]}t) 

ijkl=l,2 

®TrpJA(^)[M/f)(pb)]{A(^)[^;^)(pb)]}t). (4.9) 

Following Peres et al. [PST02j . we compute first the Lorentz transformed density matrix 
of state ^I'l, after tracing out the momentum. The expression is given, to leading order of 
w/m <C 1, by 

IVp[AM/i(AvI/i)t] = \(^^ + ""'^ ^ j , (4.10) 

where 

and cosh a := 7 = (1 — /3^)~^/^. Demanding w/m to be small does not restrict the generality 
of our results, on the contrary, it is a requirement of physical consistency in our model. 
First, the Newton- Wigner localization problem prevents us from considering momentum 
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distributions with w ^ m. In that case, particle creation would manifest and our model, 
relying on a bipartite state of the Fock space, would break down. Second, w ~ m would 
produce a very fast wave-packet spreading, having as a consequence an undesired particle 
delocalization. Finally, thinner packets favor the possibility of finding isoentangled and 
isodistillable states. 

This can be straightforwardly generalized to the other three tensor products involving 
"^i and *2, 



TVp[A^'2(A^'2)^] 

TVp[A^'i(A^'2)^] 
TVp[A^'2(A^'i)t] 



1 / 1 - \ 
2\ 1 + n'J' 

if 1 + 

2(^-(l-n'J 

l( -(!-<) 

2 I 1 + 



(4.12) 
(4.13) 
(4.14) 



With the help of Eqs. (j4.9ti4.14l) . it is possible to compute the effects of the Lorentz trans- 
formation on any density matrix of two spin-i particles, after tracing out the momentum, 
for a boost in the x-direction. 



4.2 Application to Werner states 



With these tools, we compute now the Lorentz transformation, under a boost, of the state 
in Eq. (j4.2p and obtain the reduced spin state 



/ 1 I /2 
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/ 2 
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/ 2 
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(4.15) 



where cp := We can apply now the positive partial transpose (PPT) criterion |Per96t 

IHHH96j to know whether this state is entangled. Due to the box-inside-box structure of 
Eq. (j4.15p . it is possible to diagonalize its partial transpose in a simple manner, finding the 
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following eigenvalues 

2F + 1 

rci = 

6 

1-F l-AF ,2 

xo = \ 

3 6^ 

1-F l-AF ,2 
= ^ —""^ 

X, = (4.16) 

Given that F > 0, xi and X4 are always positive, and X3 can be shown to be positive 
for < n'jj < 1. The eigenvalue X2 is negative if, and only if, F > N'^, where A'^^ := 
(2 + n'^)/{2 + 4nl,^). The latter implies that in the interval 

\<F<N', (4.17) 

distillability of state |^'~) is possible for the spin state in A and B, but impossible for the 
spin state in A' and B', both described in frame S. We plot in Fig. 14.21 the behavior of 
N'^ as a function of the rapidity a. The region below the curve (ND) corresponds to the F 
values for which distillation is not possible in the Lorentz transformed frame. On the other 
hand, the region above the curve (D), corresponds to states which are distillable for the 
corresponding values of n'^. Notice that there are values of F for which the Werner states 
are weak isodistillable and weak isoentangled, corresponding to the states in the region D 
above the curve for the considered range of n'^. On the other hand, there are states that 
will change from distillable (entangled) into separable for a certain value of n'^ , showing the 
relativity of distillability and separability. 



4.3 Analysis of existence of SIE and SID states 

The study of strongly isoentangled and strongly isodistillable two-spin states is a much 
harder task that will depend on the entanglement measure we choose. We believe that 
these cases impose demanding conditions and, probably, this kind of states does not exist. 
However we would like to give a plausibility argument to justify this conjecture. Our 
argument is based on two mathematical points: (i) analytic continuation is a mathematical 
tool that allows to extend the analytic behavior of a function to a region where it was 
not initially defined, and (ii) an analytic function is either constant or it changes along all 
its interval of definition. Point (i) will allow us to extend analytically our calculation to 
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Cf 1 2 3 



ci 

Figure 4.2: A^^ of Eq. (|4.17p as function of the rapidity a, for w/2m = 0.1. 

= 0, an unphysical but mathematically convenient limit. Point (ii) will be applied to 
any well-behaved entanglement measure. We consider then a general spin density matrix 



/ ai hi 62 &3\ 

b\ 02 ci 



V^3 



as 
d* 



C2 

d 

04/ 



(4.18) 



where ai, a2, 03, and 04 are real, and = 1. The analytic continuation of the Lorentz 

transformed state, according to our procedure in Eqs. (I4.9M.141) . in the limit n'^ — > 0, is 
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-i(Qb2+Qc2) 
2 

2 



2 

1/4 
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(4.19) 
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where 3ft and 9 denote the real and imaginary parts. This state is separable because its 
eigenvalues, given by 



coincide with the corresponding ones for the partial transpose matrix. In this case, Ai <-> A4, 
and A2 ^ A3. So, according to the PPT criterion, the analytic continuation of the Lorentz 



Our analytic calculation holds for n'^ < 1, leaving out of reach the case n'^ = 0. However, 
any analytic measure of entanglement, due to this behavior of the analytic continuation at 
n'-, = 0, is forced to change with re^ for n'^ < 1, except for states separable in all frames. 
In this way, our approximation techniques give evidence of the non-existence of strong 
isoentangled and strong isodistillable states, for variations of the parameter n'^. 

From a broader perspective, our analysis corresponds to the invariance of entanglement 
and distillability of a two spin-i system under Lorentz- Wigner completely positive (CP) 
maps |HHH96j . This is an important problem that, to our knowledge, has not received much 
attention in quantum information theory, and that will require a separate and more abstract 
analysis. Moreover, for higher dimensional spaces, like a two spin-1 system (qutrits), the 
notion of relativity of bound entanglement will also arise |HHH98j . 

In summary, the concepts of weak and strong isoentantangled and isodistillable states 
were introduced, which should help to understand the relationship between special relativity 
and quantum information theory. The study of Werner states allowed us to show that 
distillability is a relative concept depending on the frame in which it is considered. We have 
proven the existence of weak isoentangled and weak isodistillable states in our range of 
validity of the parameter n'^. We also conjectured the non-existence of strong isoentangled 
and strong isodistillable states in special relativity. We give evidence for this surprising 
result relying on the analytic continuation of the Lorentz transformed spin density matrix 
for a general two spin-^ particle state. This analytic continuation converges to a separable 
state in the limit 0. Thus, any entanglement measure analytic in the parameter is 

forced to change with n'^ in every interval, unless the state is separable in every frame. 




(4.20) 



transformed density matrix of all the two spin- 2 particle states with Gaussian momentum 
part factorizable considered here, converges to a separable state in the limit of — > 0. 



Chapter 5 



Dirac equation and relativistic 
effects in a single trapped ion 

5.1 Dirac equation simulation in a trapped ion 

The search for a fully relativistic Schrodinger equation gave rise to the Klein-Gordon and 
Dirac equations [Sak67j . P. A. M. Dirac looked for a Lorentz-covariant wave equation that 
is linear in spatial and time derivatives, expecting that the interpretation of the square 
wave function as a probability density holds. As a result, he obtained a fully covariant wave 
equation for a spin-1/2 massive particle (fermion), which incorporated ab initio the spin 
degree of freedom. It is known that the Dirac formalism plays a central role in the context 
of quantum field theory, where creation and annihilation of particles are allowed. However, 
the one-particle solutions of the Dirac equation in relativistic quantum mechanics predict 
some astonishing effects, like Zitterbewegung, Klein's paradox, Wigner rotations, and mass 
acquisition through spontaneous symmetry breaking produced by a Higgs boson. 

In recent years, a growing interest has appeared regarding simulations of relativis- 
tic effects in controllable physical systems. Some examples are the simulation of Un- 
ruh effect in trapped ions [ADM05] . Zitterbewegung for massless fermions in solid state 
physics |SLW05] . and the simulation of black-hole properties in the realm of Bose-Einstein 
condensates [GACZOO] . In this respect, simulation of the paradigmatic Dirac equation in 
a completely flexible and controllable physical system would be desirable, given its many 
counterintuitive and unverified predictions. This problem arises mainly due to the fact 
that the required range of energies (or frequencies) is of the order of the fermion mass m, 
unaccessible to present experiments. 



60 



Dirac equation and relativistic effects in a single trapped ion 



On the other hand, the fresh dialog between quantum information and special relativ- 
ity |PST021 IAM021 IPS03t IPT04j has raised important issues concerning the content and 
transfer of quantum information under Lorentz transformations. In this sense, the so-called 
Wigner rotations entangle the spin with the momentum degrees of freedom, yielding a 
non-covariant spin entropy for a single particle and reducing the degree of entanglement in 
multiparticle systems |LMDS06] . 

In this chapter, I study the simulation of the Dirac equation in a single trapped ion. We 
show |LLSS07b] that it is possible to implement realistic interactions on four metastable 
ionic internal levels, coupled to the motional degrees of freedom, so as to reproduce such a 
fundamental quantum relativistic wave equation. In this manner, the Dirac dynamics could 
be fully simulated and many relevant quantum relativistic effects reproduced and measured. 

We consider a single ion inside a Paul trap, where four metastable ionic internal states, 
\a),\b),\c), and \d), may be coupled pairwise to the center-of-mass (CM) motional degrees of 
freedom in directions x, y, and z. We will make use of three standard interactions in trapped 
ion technology, allowing for the coherent control of the vibronic dynamics |LBMW03] . First, 
a carrier interaction consisting of a coherent driving acting on any pair of internal levels, 
while leaving untouched the motional degrees of freedom. It can be described by the Hamil- 
tonian H^^^ = hflxif^^ + ) = h^xO'x, where and a~ are the raising and lowering 
spin-1/2 operators, respectively, and is the coupling strength. The driving phases and 
frequencies could be adjusted so as to produce, alternatively, the Hamiltonians H^j^ = MlyGy 
and H(j^ = hQz^^z, where cr^, ay, and az, are Pauli operators in the conventional directions 
X, y, and z. Second, a Jaynes-Cummings (JC) interaction, usually called red-sideband 
excitation, consisting of a coherent driving acting on two internal levels and one of the 
CM modes. Typically, a resonant JC coupling induces an excitation in the internal levels 
while producing a deexcitation of the motional harmonic oscillator, and viceversa. The 
JC Hamiltonian reads iJjc = hg{a^a + a^a)), where a and a) are the annihilation and 
creation operators associated with a harmonic oscillator, and g is the effective coupling 
strength. Third, an anti-JC (AJC) interaction, consisting of a JC-like coupling tuned to 
the blue motional sideband with Hamiltonian Hajc = fig{(T^a) + (T~a). Here, an internal 
level excitation accompanies an excitation in the considered motional degree of freedom, 
and viceversa. 

Later, we will give details on the experimental aspects of our proposal. For the moment, 
we want to stress that all these interactions could be applied simultaneously and addressed to 
different pairs of internal levels coupled to different CM modes. For example, it is possible 
to combine a JC and an anti-JC dynamics to form the Hamiltonian H^^ = hgx(Jx{ax + 
al) = hgxCTxX, whose physics is far from being described in terms of Rabi oscillations 
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associated with independent JC or anti-JC interactions |Zhe981 ISdMPZOTj IS AW03j . In 
turn, it yields a conditional displacement in the motional degrees of freedom depending 
on the internal state, producing the so-called Schrodinger cat states. By manipulating 
directions and phases, we could also produce H^y = hgyayy, H^^ = hgzCTzZ, and i/aiJi = 
hgxCTxiax - oD/i = figx(^xPx-, H^l = ^QyCTyPy, H^l = Hg^a^pz, as shown in an impressive 
recent experiment |HBD"'"05| . Although conventional, the directions of the Pauli matrices 
could be combined so as to produce, for instance, simultaneously or independently, H§y = 
hgyx(^yPx and Hll = hg^ya^Py 

We define the wave vector associated with the four ionic internal levels as 1^*) = ^a\o) + 
+ ^c|c) + ^'dM), that is. 



We build now the following Hamiltonian acting on j^*), 

= figf{\a){b\ + \b){a\){a., - al)/{V2i) - hgf{\c){d\ 
+h^T{\a){c\ + \c){a\) + h^,^{\h){d\ + \d){h\), 



(5.1) 



\d){c\){a.-al)/{^i) 



(5.2) 



that can be written as 



bd 
'x ' 



(5.3) 



where we have assumed that g'^ = g^'^ = g and Q'^^ = Q'^ = 17, and we have used the 
relationship 



^x -Px 



V2i h 

We may rewrite Eq. ()5.2p in a more compact manner. 



(5.4) 



o'xPxC mc 
^2 



mc 



CTxPxC 



(5.5) 



where each entry represents a 2 x 2 matrix, c := hgAx and mc^ := are the simulations 
of the speed of light and the electron's rest energy in our implementation. The Schrodinger 
equation associated with this evolution, -ffol^) = ifu)\^) / dt^ yields the same dynamics 
as the one- dimensional Dirac equation for a free spin-1/2 particle, where |^') represents 
the four-component Dirac bispinor. Extensions to the cases of 2D or 3D Dirac equations 
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are naturally expected. For example, the standard Dirac equation could be obtained by 
replacing and implementing axPx — cr.p = axPx + <^yPy + (^zPz- Given the flexibility of 
coherent control in trapped ions, we could consider the use of different couplings for the 
terms UxPx-, (^yPy, and UzPz-, producing an anisotropic Dirac equation. The properties of 
these new families of Dirac equations, not given by nature but implemented artificially in 
our proposed scheme, are still to be studied. At this point, we will turn our attention to 
describe the diverse quantum-relativistic effects that could simulated and tested under the 
dynamics described above. 

We will be hereafter using the 4-spinor notation and we shall work in the chiral repre- 
sentation, where 

3 

Hp = cy^^piaj + /3mc^, (5.6) 

with 

a = diag((T, -(t),/3 = off-diag(/2, -^2)- (5.7) 

In terms of the internal states |a), |c), and Id), we have 

ax = \a){b\ + \h){a\-\c){d\-\d){cl (5.8) 

ay = i{-\a){b\ + \b){a\ + \c){d\-\d){c\) , (5.9) 

az = \a){a\-\b){b\-\c){c\ + \d){d\, (5.10) 

(3 = |a)(c| + |6)(d| + |c)(a| + |d)(6|. (5.11) 

In the chiral representation ^ = ((/?, x)^, — '^^j (0;X)^ = ■^^j where u {v) corre- 

sponds to a positive (negative) energy state in the standard representation. 

It is interesting to note that the Dirac equation holds even in two spatial dimensions. 
This will be of relevance in our simulations given that the control in trapped ions is, cur- 
rently, mainly constrained to two dimensional]. In fact, Eq. ()5.6|) remains as it is in the 
2 + \ case, but with the index i running from 1 to 2 only (x and y components). In 3-|-l 
dimensions the Lorentz group is generated by 3 rotations M = (M^^, M^^, M^'^) and 3 
boosts N = (M°\M°2,M°3). The linearly independent combinations J = ^(M - zN) and 
K = |(M -|- iN) generate two independent SU2 algebras whose eigenvalues = j{j + 1) 
and = k{k + l) serve to label (j, k) the representations of the group. As parity commutes 
with M and anticommutes with N, it trades J by K. Hence, it is necessary to introduce di- 
rect sums of representations {j, k) © {k,j) (doubling the number of components) to account 

^What we mean with this assertion is that the interaction term a^p^ could not be easily implemented 
with current technology, in opposition to the (JxPx and (JyPy terms. The latter are just conditional rotations, 
while the former is a conditional energy shift which is harder to implement, like with an AC-Stark shift. 
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for parity. In particular, Dirac fermions belong to the (1/2,0) © (0, 1/2) representation of 
the group. It is straightforward to show that, by using the Pauli matrices to represent SU2 
for j = 1/2, the group generators are given by M^*^ = ^[7^,7'^]. Explicitly 

W = U'' '.],K = U^' r). (5.12) 

Notice that N* = and that from these expressions we get 7" = /3 and 7* = as given 
in (|5.7p : incidentally, /? acts as the parity operator in this representation. 

In 2+1 dimensions there are only two boosts and N"^ and the rotation M^^ in the 
X, y plane. There is only one SU2 algebra |Bar471 IBF65j generated by the combinations 
F' = {K^ -J\K^- J2, + J3). By parity F^'^ ^ -F^'^ and F^ F^, which is another 
representation of the same algebra with the same eigenvalues. Hence, the 2+1 Lorentz group 
only requires two component spinors for the fundamental representation. In particular this 
means that Dirac fermions will be in a sum of irreducible representations of this group. It 
is important to recognize that the behavior of 2+1 parity will continue to be represented 
by /? and the boosts by |a*, i = 1, 2. From now on we will analyze the 2+1 Dirac equation 
physics. 



5.2 Relativistic effects to be simulated 
5.2.1 Zitterbewegung 

The first phenomenon we analyze is the electron's Zitterbewegung (ZB). It is the helicoidal 
motion that arises due to the non-commutativity of the electron's velocity operator com- 
ponents, CQ!j. Explicitly, the electron's position operator, in the Heisenberg picture, evolves 
with time as |Sak67| 



Xi{t) = xi{0) + ^t + 
no 



ihc 



ai 



CPi \ 

Hd J 2Hr) 



exp 



-2i- 



Hot 



h 



(5.13) 



where we restrict to i = x, y as explained above. 

The resulting expression consists of an initial position, a motion proportional to time, 
and an unexpected oscillation term with an amplitude equal to the Compton wavelength. 
That oscillation term is the ZB. 
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We consider an initial electron wavepacket having associated positive and negative en- 
ergy components |Sak67j {h = 1), in the Schrodinger picture, 

V I I r=l,2 

+g-{p) ^^^'■^(p)e(^P-^+^l^l*)], (5.14) 

r=3,4 

where u^^\p) are the usual fermionic spinors verifying H^u^^^ = \E\u^^\ r = 1,2, Hdu^^^ = 
— \E\u^'^\ r = 3,4. g±{p) are the corresponding distributions associated to positive (+) 
or negative (-) energy components, which may be Gaussians, for example. We will take 
S+iP) — S-iP) —■ Qip)- The value of the position operator in the 2-dimensional x — y 
plane, averaged with the state (|5.14p is, in the limi1§| \p\ <C mc |Sak67| . 

x|a(p)P{x(^-)V,x(^+)[e-2*l^l* - 1] - H.c}, (5.15) 

where i = x,y and x^'**'' ^tre the two-component Pauli spinors associated to u^^ °^ (-I-) 
and u^^ " (-). Whenever the matrix element x^^'^"^ ^iX^^^^ is different from zero, which 
happens for properly chosen initial states (j5.14|) (for example, when x^^~'^ = X^'^"'''' is an 
eigenvector of cTj), then (x(t)) will oscillate harmonically in the x — y plane. In our ionic 
simulation, we have that {x{t)) will be the average position of the ion, which will make the 
ZB oscillation (j5.15p for initial ionic state given by Eq. (|5.14p and the dynamics of Eq. 
(j5.5p (in the 2D-case), with the substitutions c — > hgA^ and mc^ Ml. 

The only condition for the validity of Eq. (j5.15p is \p\ <^ mc which translates, in the 
ionic case, into (7A^.|p| <C fi. This can always be attained by properly tuning the laser 
intensities and g. We point out that the limit |p| <^ mc is not necessary in order that 
the ZB takes place, and we have only considered it for the ease of the analytical expression 
(I513D . 

Notice that, in our proposed simulation, the ZB frequency could be tunable to feasible 
experimental requirements, like wzb ■= 2|i?| ~ hQ ~ 1 MHz. This is remarkable given that 
the ZB frequency for the real free electron, wzb ~ ""t-c^ ~ 10^^ Hz is far from experimental 
reach. Additionally, the ZB amplitude in our implementation could be around 10 position 
phonons, clearly measurable, while the ZB amplitude for the real free electron is around 
10~^ A, totally unattainable. 

^We point out that we take this limit in order to have a fully-analytical, easy-to-compute expression. 
However, the ZB oscillation would also take place in other energy regimes 
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5.2.2 Klein's paradox 

Our following proposal is to simulate Klein's paradox. In 1929, Klein noticed |Kle29j the 
anomalous behavior of Dirac particles in regions where a high potential energy V exists 
{Hd = cJ2i oz-iPi + I3m(? + VI/i): When V > 2mc^ , negative energy electrons may swallow 
V , acquiring positive energy and behaving as ordinary electrons, while leaving a hole in the 
Dirac sea. This stems from the fact that H£) implies (pc)^ = {E — V + m(?){E — V — mc^) 
which is positive when either both factors on the r.h.s. are positive, or both are negative. In 
the first case, E > mc^ + V which is the usual condition; in the second one E < —m(? + V , 
which may be larger than mc^ as noticed by Klein. In this case an electron - positron pair 
could be created from y. Interactions as simple as -ff/ = V{\a){a\+\b){h\ + \c){c\+\d){d\)Q{t— 
to) can simulate this phenomenon. Notice that Hj does not induce any transition by itself. 
It is the tuning of F to a value larger than 2m(P' at t = to that produces the effect. This 
phenomenon could be simulated in our implementation with a global Stark shift. 



5.2.3 Wigner rotations 

The next physical phenomenon we consider is the Wigner rotation induced by Lorentz 
boosts. We will be using mainly the notation of Ref. |AM02j . We begin with a particle of 
s = \ with momentum pz = mc sinhr/ and energy Ez = mc? cosh. r]. We make a Lorentz 
boost along direction x, with rapidity u;, being tanh(— w) = Vx/c, to analyze the explicit 
dependence of the Wigner rotation on the momentum pz of the particle and on the boost 
velocity Vx- The Wigner angle 0p that the spin rotates conditional to the momentum is 



^ sinhrysinhu; ;^sinhti; /^-^x 

tan e„ = — — = . (5.16) 

cosh rj + cosh u + cosh uo 



On the other hand, we have sinhw/ coshw =: tanhw and cosh uj — sinh a; = 1, which gives 

. , 2 tanh^ u) vi . ^ 

smh^w = = ^ ^ , (5.17) 

1 - tanh^ UJ -vl' 

(5.18) 

According to this, and considering pz ^ mc, Vx ^ c, we have 

0p--^-. (5.19) 
zmc c 

This interaction may be simulated with a Hamiltonian of the kind H^l = hgzAzCzPz as 
introduced above, where here hgzAz = hvx/i^mc^t) = Vx/{2Qt), being t the elapsed time. 
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5.2.4 Spontaneous symmetry breaking induced by a Higgs boson 

The last phenomenon we consider is the spontaneous symmetry breaking (SSB) that gives 
the electron a finite mass through the acquisition of a finite vacuum expectation value (v.e.v) 
of a Higgs field |PS95j . In this case, the electron's mass m is given by the expression 

mc2 = r(0), (5.20) 

where F is a Yukawa coupling and (0) is the v.e.v. of the Higgs field. As this v.e.v. grows 
from a zero value into a finite one, the electron's mass becomes nonzero. It is possible 
with our experimental proposal to simulate general behaviors of the Higgs field that would 
produce many different dynamics of mass acquisition. Given that, according to Eq. (jS.Sp . 
mc^ := HQ, and here we have mcP = the Higgs dynamics may be simulated by tuning 

hil.{t) = T{(j){t)). An appropriate choice of i^{t), could simulate the phase transitions that 
lead to the true vacuum and massive fermions. 

5.3 Experimental procedure 

The basic ingredients to implement the proposed experiments are four independent elec- 
tronical (internal) states in one ion, harmonically confined in a rf-Paul trap. We will have to 
achieve the initialization of the states and four pairwise and independent couplings, as out- 
lined in Eq. (|5.2p . to be driven simultaneously. Finally, one will have to read out the result 
of the calculation by measuring the population of each electronic state and their distribu- 
tion over the harmonic oscillator states of motion individually. For a thorough reference 
containing detailed experimental procedures in trapped ions, see |LBMW03] . 

The required states could be composed by four ground-state hyperfine levels of an earth 
alkaline atomic ion, e.g. of ^^Mg+ by |-F = 3;mf = —3), \F = 3;mf = —2), \F = 2;mf = 
—2) and \F = 2;mf = —1) % ( \a), \b), \c) and Id) respectively), as depicted in Fig. 15.31 
A constant external magnetic field will define the quantization axis and lift the degeneracy 
of the levels. We will exploit the formal equivalence between a two- level system and a 
spin- 1/2 magnetic moment in a magnetic field for each individual coupling of two out of the 
four states. To achieve the required pairwise and independent couplings, we would suggest 
two photon stimulated Raman transitions. Choosing the beam orientation relative to each 
other and relative to the quantization axis appropriately will allow to provide beams of the 

^The nuclear spin of / = 5/2 of ^^Mg+ would allow for 3 ancilla levels needed for detection purposes 
described in the following. 
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required polarization to drive all necessary transitions and to avoid disturbing AC-Stark- 
shifts. 




Figure 5.1: Schematic of relevant energy levels of one ^^Mg"*" ion (not to scale). Shown 
are the ground-state hyperfine levels supplying the four ququat states (|a), |c) and \d)) 
and the harmonic oscillator levels due to the harmonic confinement in a trap similar to 
that described in |LBMW03] . We subsumed excited levels of the P1/2 and P3/2 states. 
Typically, the energy splitting of the motional levels and the Zeeman shift induced by an 
external magnetic field are of the same order of magnitude within 1-10 MHz, therefore 
much smaller than the hyperfine splitting of 1.8 GHz, the fine structure splitting of 2750 
GHz and the optical transition frequency of the order of 1015 THz. We depict the resonant 
transition state sensitive detection and the two relevant types of two photon stimulated 
Raman transitions between states \a){c\ and |c)(d|. The required transitions between states 
\b){d\ and \a){b\ would be provided by additional beams of slightly different frequencies or 
detunings. 

At the start of each experiment, the ions will be laser cooled close to the motional 
ground state and optically pumped into state \a). Two carrier transitions and two red/blue 
side band transitions will be driven simultaneously to realize the calculation, one of each 
type depicted in Fig. 15.31 The coupling of states \a){c\ would be equivalent to that of 
the not shown \b){d\ via an additional carrier. The red/blue sideband transitions between 
\a){b\ would be driven by an additional pair of beams comparable but at slightly different 
frequencies or detunings to those shown driving |c)((i|. 
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At the end of the calculation, state sensitive detection is to be realized by an additional 
laser beam, tuned to a cycling transition |LBMW03] . coupling state |o) resonantly to the 
P3/2 level. Since the Zeeman shift (1-lOMHz) will not be sufficient in comparison to the 
line width of the resonant transition (2p 43 MHz) to protect the population in state \b) 
against this interfering excitation, we will have to hide it in one of the F = 2 ground state 
levels prior to detection. Thus, the population of state \d) would be transferred to state 
\F = 2, mf = 0) to make room for the population of state \b) being transferred subsequently. 
To measure the population of any of the other three states we propose to use the carrier 
transition between states \a){b\ or \c){d\ instead of the side bands (and an additional carrier 
to drive |a)(c|) to allow for the transfer into state \a). Since these carrier transitions do 
not affect the motional state populations, hiding the population of state |c) and driving a 
blue-sideband transition between the population transferred into state \a) and the cleared 
state |c) would allow for deriving the motional state population. 

Taking into account the available laser intensities, all necessary Raman beams could 
be derived out of only one laser source. We would have to split the original beam and 
provide the necessary frequency offsets and switching via multi passing through acousto- 
optical modulators (AOM). The amount of necessary laser beams could be further reduced 
by using electro-optical modulators to provide red and blue sidebands at the same time. 
The additionally appearing carrier would just have to be book kept. There appear to be 
no fundamental reasons, why this experiment could not be realized with state of the art 
technology and fidelities of operation at the expense of only a slightly increased effort for the 
additionally required beams. We point out that, in order to facilitate the experiment, one 
could consider just the Dirac equation in 1+1 dimensions. The Hamiltonian to simulate in 
that case would have the form Hd = caxPx + fnc^f^y) that could be easily implemented with 
one carrier, one blue- and one red-sideband transitions. We remark that the predictions 
of Dirac equation here reviewed, like Zitterbewegung, Klein's paradox, Wigner rotations or 
mass acquisition would already take place at this level of 1+1 dimensions. 
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Entanglement of continuous 

variables 



Chapter 6 



How much entanglement can be 
generated between two atoms by 
detecting photons? 



In the world of quantum information processing it is widely accepted that, while photons 
are the ideal candidates for transmitting quantum information, this information is better 
stored and manipulated using atomic systems. The reason is that, while photons can be 
moved through long distances with little decoherence, atoms can be easily confined and can 
preserve quantum information for a long time. Consequently, an ideal design for a quantum 
network will conceivably be built upon a number of atomic or solid state devices which 
communicate through photonic quantum channels. 

There exist mainly two methods for entangling distant atoms. The first one is based 
on emission of photons by the first atom, which afterwards interact with the second atom 
generating the entanglement |(]ZKM97[ lvE(]Z97[ lGST+981 IPKOOl ICLEn2[ K^LEOSj . On the 
other hand, the second method relies on detecting the photons emitted by the two atoms 
with the subsequent entanglement generation due to interference in the measurement process 
|CCGFZ99l IBKPV99] . Some aspects of these proposals have been experimentally verified 
IJKPOll [KHR021 IMBB+Ol IBMDM041 IVWS+06l fMCJ+OGi [MMM+061 ICdF+OSl IDJD+OS" 



IBJD"'"06] . Most of the experiments with isolated atoms and light aim at entangling the 
internal state of the atom with the polarization of the photon |KHR021 [MBB+ 04[ [BMDM04[ 
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IVWS+061 lMCJ+06[ lMMM+06] . It is clear that due to the size of the Hilbert space, the 
maximum attainable entanglement is one ebit. 

In this chapter we will deal with the generation of entanglement between two atoms. 
We will focus on the second method mentioned above, in which entanglement is generated 
by measurements. To avoid the limit of one ebit, we work with continuous variables and 
seek entanglement in the motional state of the atoms. We will answer two fundamental 
questions: How much entanglement can be produced between the atoms? How can we 
achieve it? The basics of this chapter is contained in Ref. |LGRC07] . 

The first result, obtained in Section 16.11 is that by usual means — two atoms, one or 
two emitted photons, linear optics and postselection [CCGFZ991 |MMM"'"06 — , we cannot 



produce more than 1 ebit of entanglement between the atoms, even if our Hilbert space 
is larger. The second result, exposed in Section 16.21 is that we can achieve an arbitrary 
amount of entanglement using at least two emitted photons and what we call an Entangling 
Two-Photon Detector (ETPD). The ETPD is a device which combines both photons in a 
projective measurement onto a highly entangled state. Theoretically, an ETPD could be 
built using a Kerr medium and postselection, but current nonlinear materials are too inef- 
ficient for such implementation. Inspired by the KLM proposal [KLMOlj . we demonstrate 
an efficient scheme for simulating the ETPD using ancillary photons. In Section 16.31 we 
clarify our proposal with a simple example: three photons and three detectors. Our last 
result, presented in Section [67^ is that introducing N — 2 additional photons in our setup, 
together with A'^ single-photon detectors, beam-splitters and an attenuator, we can obtain 
an amount of entanglement as high as S" = log2 N ebits. Finally, at the end of the chapter 
we discuss the relevance of these results and possible implementations. 



6.1 Entanglement based on atoms and photons 

As mentioned above, current experiments with isolated atoms and light aim at entangling 
the internal state of the atom with the polarization of the photon [VWS"'"d6] . It is clear that 
in both cases, due to the size of the Hilbert space, the maximum attainable entanglement is 
one ebit. We will rather work with continuous variables and seek entanglement between the 
motional state of the atom and the photon. We have in mind the setup in Ref. |CCGFZ99j 
where two atoms are excited with a very small probability and we consider the state of the 
atoms after spontaneous emission, when both are in the ground state. The state of the 
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Figure 6.1: Schema of possible experiments for entangling two atoms, (a) Only one photon 
detected, but we do not know from which atom, (b) Two photons are detected, one from 
each atom, (c) Three photons are detected, one being supplied by the experiment (dashed 
line). Due to the setup, the probabilities of reaching each detector are balanced and the 
detectors do not distinguish between left- and right-coming photons, (d) Entangling Two- 
Photon Detector gedanken experiment. By detecting only a range of momenta we entangle 
the momenta of the atoms, pi± + p2± — 0. 



system is given by the wave function 

I*) ~ ejdpal\v&c)igiip)\-p,0) + g2ip)\0,-p)) + 

+ j dpidp2 ^i(pi)&(p2)4i42|vac)| - pi, -p2) 

+ |vac)|0,0) +0(e3). (6.1) 

Here p, pi, p2 denote the momenta of the emitted photons; cip,pi,p2 their associated creation 
operators and |vac) the vacuum state of the EM field, and e <C 1 is the excitation probability 
of the atoms. The initial momentum distribution of the emitted photons is given by Qi{p) 
for the i-th. atom. As we will sec later, we require some uncertainty in the initial momentum 
in order to generate a large amount of entanglement. Finally j — p, 0), jO, — p) and | —pi, — P2) 
denote the recoil momenta of the atoms after emitting the photons. The terms omitted in 
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Eq. (|6.1|) correspond to higher order processes where an atom emits more than one photon. 
These terms will have a very small contribution if the decay time of the atom is longer than 
the duration of the exciting pulse. 

Let us now consider a single detector placed symmetrically below the atoms |CCGFZ99] . 
as in Fig. 16.1b ,. If there is one single photon detection, this will amount to a projective 
measurement onto a single-photon state and out of the state in Eq. (|6.1|) only the term on 
the first row will survive. Since the photons coming from the atoms are undistinguishable, 
an implicit symmetrization will take place and the final state of the atoms will be of the 
form |i/;i)|0) + |0)|'(/'2)7 for some motional states ipi and ip2- Even though we work with 
continuous variables, this state can have at most 1 ebit which corresponds to {ipi\ip2) = 

(V^ilo) = {iP2\o) = o. 

We are going to show now that with two emitted photons, linear optics and two detectors, 
we cannot do better than one ebit of entanglement [See Fig. 16.1b ]. The proof generalizes the 
previous argument with a little bit more care. First of all, linear optics amounts to a linear 
transformation of the initial momentum modes, Op, to new operators, b^(j}) ■= U^apU\. 
A trivial example of this is a 50% beam splitter, which changes the photons from incident 
states a+p and a_p to (a-|_pibe*'^a_p) / \/2. Linear optics can be combined with measurements. 
Without loss of generality, all measurements will take place at the end of the process and 
they amount to a projection onto the modes a^i and for the first and second detector, 
respectively. The state after a projective measurement onto two single-photon detectors 
reads 

1*1?) = j C^Pl'^P2^l(Pl)^2(P2)(axiaa;2ft^,i(p,)f'^,2(p2))vac| -Pl,-P2>- (6.2) 

Note that the modes a^^ and detected by the first and second detector are expressed 
on an orthonormal basis different from that of the Op or b operators. We enclose this 
information, plus the initial wavefunction of the photon in the following c-numbers 

fj{xi,pj) := gj{pj)[ax,,h\^^^^^]. (6.3) 

Using these wavefunctions we define the motional states 

iV'ii) := / dpfj{x,,v)\p)- (6.4) 



The expectation value in Eq. (j6.2|) can be written in terms of the fj{xi,pj). We thus arrive 
to the following expression for the atomic state after the measurement 

|*i?)oc|Vll)|^22) + |V'2l)|V'l2). (6.5) 
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This state cannot have more than 1 ebit of entanglement, which happens when all the states 
^11 1 "012) V'21 and are orthogonal to each other. 

We must make several remarks. First of all, adding more detectors does not improve the 
outcome. Second, our proof is valid independently of the number of beam splitters, prisms, 
lenses and even polarizers we use. In particular, attenuating elements such as polarizers 
and filters can be treated as a linear operation plus a measurement and are covered by the 
previous formalism. 



6.2 Arbitrary degree of entanglement between two atoms us- 
ing postselection 

We propose now to use an ETPD to obtain an arbitrary degree of entanglement between 
the two atoms. An ETPD is by definition a device that clicks whenever two photons 
arrive simultaneously and with their momenta satisfying a certain constraint. An example 
would be a parametric up-conversion crystal, in which pairs of photons with momenta 
pi and p2 are converted with a certain probability into a new photon with momentum 
p = Pi +P2- One imposes a constraint on the initial state by post-selecting a window of final 
momenta. For example, restricting the measurement to photons with transverse momentum 
Pi. = 0) then the initial contributing momenta must be those satisfying p±i + p_|_2 = 
[Fig. 16.1b ]. In this example the ETPD ideally projects the initial two photon product state 
l^ph)~/ '^Pidp2Qi{pi)G2{P2)ci\)io)p2\\Sic) onto the probably entangled state 

ETPD\_/" , ,1^,^!^ J 



l^ph )=j dpadpbg{Pa,Pb)ala%b\"^^'^)- (6-6) 

Here g{pa,Pb) is the acceptance function of the detector or, equivalently, the constraint that 
the final detected momenta pa and pb obey. 

We claim now that with two emitted photons, linear operations and an ETPD, there is 
no limit to the attainable entanglement. To prove it we consider that after projection of the 
photon part of state in Eq. (|6.ip into |^ph"^^)=/ dpadpbg{pa,Pb)(4'aC'Pb\^^'^)^ resulting 
atomic state will take the form 

l^at^™> = / dpadpbg{Pa,Pb)\'i'{Pa,Pb)) (6.7) 

with the already entangled state 

\^{Pa,Pb)) ■■= j dpidp2[fl{Pa,Pl)f2{Pb,P2) + fl{Pb,Pl)f2{Pa,P2)]\ -Pl,-P2)- (6.8) 
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Depending on the specific shape of the functions g{pa,Pb) and fiipuPi) I = a,b, i = 1,2, 
the corresponding state may reach an unbounded degree of entanglement. For example, 
let us consider that the photons evolve freely in space without any linear optics elements, 
fiiP^Pi) = Qi{Pi)^{p — Pi)-: and assume that the detector has a very narrow acceptance 
function g{pa,Pb) = ^iPa +Pb)- The wider the initial momentum widths of the two photons, 
the larger the resulting bipartite atomic entanglement, which is not bounded from above. 
Indeed, in this ideal case the outcome will be much like the EPR pairs from the seminal 
paper Ref. [EPR35j . 

Current Kerr media are too inefficient to practically implement the ETPD introduced 
here. Motivated by this we have designed another protocol that simulates the outcome of 
an ETPD using linear optics, additional photons and postselection. As shown in the KLM 
proposal [KLMOlj . any highly entangling quantum gate can be performed this way. Some 
care is needed, though, because for our problem direct application of the KLM proposal 
results in a too large number of gates to even obtain a moderate amount of entanglement 
between both atoms. 

Our proposal starts up from the two atoms after having emitted two photons, and we 
add N — 2 additional ancillary photons. 



^ ) = / dpidp2...dpNGi{pi)G2{P2)---GN{PN)al^al^...al Ivac) ^\ -Pi,-P2). (6.9) 



The resulting state after linear operations on the photons, and A^-fold coincidence count 
on the detectors, will be, analogously to the two-photon and two-detector case [Eqs. 



where Hn denotes the set of permutations of A^ elements. This state may contain much 
more than one ebit of entanglement. In fact, an upper bound to the degree of attainable 
entanglement is 5 = log2 A^ ebits. We will show afterwards that this bound is indeed 
saturated. 

6.3 Three photons and three detectors 

As a clarifying example we consider the setup in Fig. 16.1b with three photons and three 
detectors. Photons PI and P2 come from their respective atoms, we introduce a single 
auxiliary photon, P3 and we place three detectors symmetrically to the atoms, Xi, X2, X^. 



([62])- ([63])] 




(6.10) 
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Figure 6.2: Outcome for an experiment with two atoms and three photons, as shown in 
Eq. (lUTTTD . 



The final state for the two atoms, considering that ah the three detectors are excited by 
the three photons, and fixing relative phases equal to 1 for simplicity purposes, would have 
the form 



1^- 



(3)^ 
at / 



-^(|1,2) + |2,3) + |3,1) + |1,3) + |3,2) + |2,1)), 



(6.11) 



where we denote with \i,j) the atomic state associated to detection of Pi in Xi, and P2 in 
Xj. In Fig. 16.21 we show the A^! = 6 processes that contribute coherently to the two-atom 
final entangled state. The previous state contains an entanglement of 5 = 1.25 ebits. 



6.4 Saturation of the entanglement bound log2 N ebits 

The previous example is suboptimal. The maximal amount of entanglement of S" = log2 
ebits is reachable for some of the states in Eq. (j6.10p . To prove it we consider a very 
symmetric configuration in which the detectors are located along a circle, equidistant to 
both atoms [Fig. 16.1b ]. We will assume for simplicity that the two emitted photons and 
the N — 2 ancillary ones are in s-wave states and arrive with equal probability and phase 
to every detector. In a similar fashion as in Eq. (j6.1ip . the final bipartite atomic state will 
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X, 




Figure 6.3: Quantum circuit for saturating the bound of log2 ebits as described in the 
text. 



take the form 

\^T) = Y^a,\iJ) oc - 6,,)\iJ), (6.12) 

ij ij 

where is the final bipartite atomic state after detection of photon Pi in detector Di, 
and photon P2 in detector Dj. In matrix form, the coefficients Cij are 



Cij(xNei{ei) -InxN, 



(6.13) 



where := l/\/iV(l, 1, 1)7V- Both the reduced density matrix of one atom and the 
Schmidt rank can be obtained from this matrix. The previous state can be rewritten in the 
form 



N 



a (iV- l)ei(f?i)^-^e;(e 



(6.14) 



i=2 



where {ej}, i = 2,...,N is a completion of ei to an orthonormal basis in C^. From here 
it is obvious that the density matrix has full-rank and we can with local operations obtain 
a maximally entangled state of the form, up to local phases, Cij oc InxN- To do so we 
must reduce the contribution of the term ei. As shown in Ref. |RZBB94] . a network of 
beam-splitters and phase-shifters can be used to perform a unitary operation, Ug-^, that 
maps the mode a^^ cx X^^^ a-^. to a single optical port. If, as shown in Fig. 16.31 we place on 
that port a filter F that decreases its amplitude by a factor A^ — 1, when the A^ detectors 
click simultaneously the atoms will get projected onto a maximally entangled state with 
Cij = 6*1 (ei)^ — ^i=2^ii^i)^ ■ The proof is a little bit cumbersome, as it involves studying 
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how all the photon modes in Eq. (|6.9p transform under the nonunitary operation given by 
the network in Fig. 16.31 and then ensuring that the detection of N photons does indeed give 
rise to the maximally entangled state. The state of the N photons and 2 atoms is given, 
before detection, by 

N N N 

I*) = (E^)"^"' E «]Jj'i) E (6.15) 

j = l jl = l j2 = l 

where \ ji), i = 1, 2 is the recoil state of atom i, and aj is the mode associated to detector j. 
The linear transformation Ug-^ applied to the modes before detection (see Fig. 16. 3|) . maps 
the modes aj into the modes cj, where c| := X^jLi c^. A; / 1, are a completion to 

an orthonormal set, given by c|, := X^j^i^fej'^j- "^^^ corresponding state after application 
of Ug-^ is 

TV N 

C/e-J^) = (4)^-2 ^ je^^) ^ 4|e„,)|vac), (6.16) 

711=1 n2=l 

where | e„ . ) is the recoil state of atom i transformed under Ug^ . 

The following step of our protocol is the application of an attenuator F to mode ci, 
according to 

4^yr^+vp4. (6.17) 

The meaning of Eq. (16.170 is clear: the attenuator absorbs the photon with a probability 
1 — p while leaves it unperturbed with a probability p. Postselection when detecting the N 
photons will ensure that no photon has been absorbed in the process. 
The resulting state is given by 

N N 

FUg,\^) = (4)^-2v^^-2 ^ clVnAen,) E 4K.|e„,)|vac), (6.18) 

ni=l n2=l 

where VJ^. := ^/p for rij = 1 and 1 otherwise. 

Now we apply the inverse transformation U-^^ , giving 

N N N N N 

U^^'FUg,\^) = (E«5)'^''v^'^"' E E en,n4r^n,\en,) E E a]^ | e„, ) | vac) . 

j=l ni=lji=l ra2=lj2=l 

(6.19) 

Upon detection of the photons at the N detectors, the final atomic state will be 

TV TV 

l^at) CX ^ ^ ^ ^ C'j'ij^ ^raiji ^"2^2 I ^11 ^712 ) ' (6.20) 

ni, 712=1 ii,j2=l 
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where Cj-^j^ := 1 — Sj-^^j2 was previously introduced in (|6.12|) and appears due to the fact 
that each detector just clicks for one single photon. Cj^^j^ may be recast in the form Cj-^j^ = 
Neij^eij^ — producing an atomic state 

N N 
l^at) OC ^ {Vn^Vn^N5in^5in,^ - K12 "^nma ) I Cni e„2 ) = ^ V^{N5in - l)|e„e„). (6.21) 
n\,n2=l n=l 

On the other hand, we had that Vi = y/p, and = 1 n / 1, so that our final result is that, 
fixing p = -jyi^ the final atomic state will be the maximally entangled state 

N 

l^at) OC |eiei) -^|e„e„), (6.22) 

n=2 

which contains an entanglement of log2 N ebits. This completes our proof. 

Summing up, in this chapter we have demonstrated that it is possible to achieve an 
arbitrary amount of entanglement in the motional state of two atoms by using spontaneous 
emitted photons, linear optics and projective measurements. The resulting states can be 
used to study violation of Bell's inequalities and also as a resource for quantum information 
processing. We expect that similar ideas can be used to entangle atomic clouds, replacing 
the photons with atoms, because in this case it is easy to build a two-atom detector. 

Regarding the implementation, the ideas shown here can be tested easily in current 
experiments. We would suggest using two trapped ions as target atoms. The ions should 
be either on a very weak trap, or released right before excitation. The entanglement in the 
momentum will translate into an entanglement in the position of the atoms after a short time 
of flight. In practice, with only one additional photon, 1.58 ebits can be produced, and we 
expect a value of 2 ebits to be achievable. Clearly, due to the requirement of having single 
photons, producing a larger amount of entanglement will be more difficult with current 
technology, even though, as we have shown here, there is no fundamental limit. 
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Spin entanglement loss by local 
correlation transfer to the 
momentum 

Bipartite and multipartite entanglement is considered a basic resource in most applications 
of quantum information, communication, and technology (see for instance Refs. [NCOOl 
IGMD02] ) . Entanglement is fragile, and it is well-known that in some cases interactions with 
an environment external to the entangled systems may decrease the quantum correlations, 
degrading this valuable resource |YE02l IYE03[ IYE04[ IDH04[ iNSOSl ICMPB051 ISMDZOGj . 
Momentum acts as a very special environment which every particle possesses and cannot get 
rid of. Consider for instance a bipartite system which initially is spin-entangled and with 
the momentum distributions factorized. It will decrease its spin-spin correlations provided 
any of the two particles entangles its spin with its momentum. This simple idea was studied 
in the natural framework of special relativity, where changing the reference frame induces 
Wigner rotations that entangle each spin with its momentum |PST02t ICzaOSt IAM02t IGA02t 
IPS031 [CW03l IGB A031 iLMDSOGj . However, this is just a kinematical, frame-dependent effect 
only, and not a real dynamical interaction. On the other hand, this type of reasoning is 
also related to which-path detection |SZ97] . In addition, an experiment observing photon 
polarization disentanglement by correlation transfer, in this case to the photon's position, 
was performed [SMKROT] . Can local interactions entangling spin with momentum produce 
the loss of non-local spin-spin entanglement? To our knowledge this question has not 
been explored in the literature. More remarkably, any particle owns a certain momentum 
distribution acting as an intrinsic environment, which can never be eliminated by improving 
the experimental conditions. But how does this fact affect the spin-spin correlations? This 
is the question we want to analyze |LLS06aj in this chapter. 
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In Section 17.11 we consider a bipartite system, composed by two s = ^ fermions or two 
photons, which are initially in a Bell spin state |^'~). We use a formalism that shows the 
decrease of spin entanglement whenever an interaction locally entangling spin with momen- 
tum takes place. We obtain the negativity N |VW02j in terms of an integral depending on 
the spin rotation angle conditional to the momentum. In Section [7.21 we analyze this phys- 
ical phenomenon in two specific situations: (i) Two spin-| fermions in a |^'~) Bell state, 
with Gaussian momentum distributions, that fly apart while one of them passes through a 
local magnetic field. Their spin entanglement decreases as a consequence of the transfer of 
correlations to the momentum of the latter fermion. And (ii) Two photons in a polarization 
|^~) Bell state, with Gaussian momentum distributions, which separate while one of them 
traverses an optically-active medium. This medium will entangle the polarization with the 
momentum and thus decrease the polarization entanglement. This is an unavoidable source 
of decoherence. In Section 17.31 we show that the apparent purely quantum communica- 
tion resulting from these procedures is not possible. Classical communication has to be 
exchanged for the protocol to operate. The chapter ends with our conclusions. 



7.1 Spin entanglement loss by correlation transfer to the mo- 
mentum 

We consider a maximally spin-entangled state for two s = | fermions A and B, or two 
photons A and B. The case we analyze is that in which the two particles are far apart 
already. This avoids dealing with symmetrization issues. Indeed, our state is the maximally 
entangled one for two s = ^ spins or polarizations, containing 1 ebit. 



I^p) := — [*(^)(pa)*f )(pb) - *f (paj^-f ^(Pb)], (7.1) 
where Pa and Pb are the corresponding momentum vectors of particles A and B, and 







*J(P) := X(p)li> = ( ^1^) ) . (7.2) 

with bimodal momentum distribution A4(p) := "^(^ppi + '^pp2)- We consider for the 
time being this kind of distribution for illustrative purposes. At the end of this section we 
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generalize our results to arbitrary momentum distributions of the two particles. ||) and ||) 
represent either spin vectors pointing up and down along the z-axis, in the fermionic case, 
or right-handed and left-handed circular polarizations, in the photonic case. If we trace out 
the momentum degrees of freedom in Eq. (jT.ip . we obtain the usual spin Bell state, l^*"). 

We consider a local interaction which entangles the spin of each particle with its momen- 
tum through a real unitary (orthogonal) transformation U. We choose a real transformation 
for the sake of simplicity and in order to obtain fully analytical results. The generalization 
for inclusion of complex phases is straightforward but adds nothing of relevance in this 
section. We will take it fully into account in Sec. 17. 2[ 

Each state vector in Eq. (|7.2|) transforms as 



cos6lp^ ^ 5pp^ ^ I cose^p^ j 5pp^ 




smt 



P2 



sin6lp2 \ 6pp^ 



cos ( 



Op^ ) V2 



(7.3) 



where 6p^ and 9p^ produce a spin-momentum entangled state whenever Op^ ^ 9p^. The 
effect of this local interaction is that a part of the non-local spin-spin entanglement is trans- 
ferred to the spin-momentum one, and the degree of entanglement of the spins decreases. 
To show this, we consider the state (j7.ip evolved with the transformation U, and trace out 
the momenta. 

= \ 5]./TrpJ[/(^)[^(^)(pa)]{C/(^)[*lf)(pa)]}t) 

s,s' 

®Trp,(C/(b)[M/5(pb)]{C/(^)[^t'(Pb)]}^), (7.4) 



where ss' := 5s,s> — ^s-s>- It can be appreciated in Eq. (j7.4l) that the expression is decom- 
posable in sum of tensor products of 2x2 spin blocks, each corresponding to each particle. 
We compute now the different blocks, corresponding to the four possible tensor products of 
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the states (|7.3p 



4 + 4 



ClSl + C2S2 



ClSl + C2S2 sf + S2 



-ClSl 
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where c,- 



cos 



and Sj 



of the local interaction U in the state 
equal interaction angles for the two particles, ^p'^ 
resulting bipartite spin state is 



sin(0p.). This way, it is possible to compute the effects 
p/ after tracing out the momentum. We choose 
— 9p^ , as a natural simplification. The 
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4*12/ 



where S12 := sin( 



'pi 



(7.5) 



6p^) and C12 := cos(0p;^ — Op^). The entanglement measure we will 



use is the negativity |VW02] . defined as N := max{0, — 2Amin}5 where Amin is the smallest 
eigenvalue of the partial transpose (PT) matrix of Eq. (17. 5p . It is very easily computable, 
and is found to be 



N 



cos 



(7.6) 

?P2 the entanglement remains 



^Pl ^P2 ) • 

From this expression it can be appreciated that for 6p^ 
maximal (1 ebit), and for 0p^ separating from 9p^ the entanglement decreases, until Op^^ — 
9p2 = vr/2, where it vanishes and the state becomes separable. We plot this behavior 
in Figure 17.11 showing the negativity N in Eq. (|7.6|) as a function of 9p^ and 9p^ . Every 
local interaction producing spin-momentum entanglement will in general diminish the initial 
maximal spin-spin entanglement of the two particles, thus degrading this resource. This 
result is valid either for s = i fermions or photons, as they both have a two-dimensional 



internal Hilbert space. The generalization of Eq. (j7.6|) to a uniform distribution with n 
different momenta is straightforward, and the spectrum of the corresponding PT matrix is 



1 1 

2' 2'^ 



1 1 " 



(7.7) 
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Figure 7.1: Negativity N in Eq. (j7.6p as a function of ^p^^ and Op^. 



with a resulting negativity 



2 



COS 



n 



We take now the continuous hmit, for an arbitrary momentum distribution ^(p) for each 
particle. We suppose for the sake of simplicity l^/''-'^^ (p)| = I^^^^P)!) although the spatial 
distributions do not overlap, as the two particles are far away. Accordingly, will be 



N 



1-2 / / dViV'(p)nv'(p')rcos 



'm2 



(7.9) 



Notice that this expression involves an integration over the momentum variables p and p', 
associated to the same particle (not to each of them). We point out that, according to (|7.9p . 
in the case where momentum does not entangle with spin (i.e. whenever 6p is a constant), 
then = 1 and thus the spin-spin entanglement remains maximal. Otherwise, the spin- 
spin entanglement would decrease due to the transfer of correlations to the spin-momentum 
part. 

The loss of spin entanglement under a spin-momentum entangling transformation can 
take place in a variety of possible situations. Wigner rotations that appear under relativistic 
change of reference frame entangle each spin with its momentum producing loss of spin- 
spin entanglement |PSTn2[ IClzanSl IAMn21 \GA()2\ IPSM \CWO^\ \GBAO^\ ILMDSnOj . This is 
just a kinematical-relativistic effect, not due to a dynamical interaction. In the rest of the 
chapter we focus on two relevant examples of these interactions, taking fully into account 



86 



Spin entanglement loss by local correlation transfer to the momentum 



the complex phases: a local homogeneous magnetic field, for the two-fermion case, and a 
local optically-active medium, for the two-photon case. 

7.2 Applications 

7.2.1 Two fermions and a local magnetic field 

In this section we analyze a bipartite system, composed by two s = ^ neutral fermions A 
and B, which are initially far apart and in a Bell spin state |^'~), with factorized Gaussian 
momentum distributions. We consider that one of them traverses a region where a finite, 
homogeneous magnetic field exists. As a result, it will transfer part of its spin correlations 
to the momentum. 

The initial spin-entangled state for the two fermions A and B is 

l^p) := -^[M/;^)(pa)^f )(pb) - (Pa)*f ^(Pb)], (7.10) 

where pa and pb are the corresponding momentum vectors of particles A and B, and 

%(p) ■■= Q{pm = (^^^^^y 

^i(p) := a(p)li) = ^H^^ ^ , (7.11) 

with Gaussian momentum distribution ^(p) := 7r~^/^(7~^/^ exp[— (p — po)^/2(7^]. In Eqs. 
(jT.lip we are not indicating explicitly the particle index. In the center of mass frame, 
Pq*^^ = — Po^'*, and we consider that the two particles are flying apart from each other. ||) 
and II) represent spin vectors pointing up and down along the z-axis, respectively. If we 
trace out momentum degrees of freedom in Eq. (|7.10p . we obtain the usual spin Bell state, 

Suppose a local interaction which entangles the spin of fermion A with its momentum 
through a unitary transformation U. In this case we choose a magnetic field Bq which is 

(a) 

constant on a bounded region D of length L, along the direction of Pq , vanishes outside 
P, and extends infinitely with a constant value along the other two orthogonal directions, 
as shown in Figure 17. 2i We take Bq along the direction orthogonal to Pq'^^ , so it is diver- 
genceless, V • Bq = 0, and we quantize the spin along Bq so that o" • Bq = sBq, with s 
the corresponding spin component. Due to the rotational invariance of the spin singlet. 
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Figure 7.2: Sketch of the two-fermion case explained in the text. Fermion A traverses a 



constant magnetic field Bq located in region V with a width L along the direction of Pq 



(a) 



this choice is completely general. In momentum space, the problem reduces to one dimen- 
sion, the one associated to the direction of Pq . We will denote from now on p to the 
corresponding momentum coordinate. 

The system Hamiltonian can be written as 

H = ^ + j{Bo-S)e{x)e{L-x), (7.12) 

where 7 is the magnetic moment of our neutral particle. Accordingly, we get 

p = i[H,p] = (-7(Bo-S)[5(x) -,5(L-x)],0,0), (7.13) 

S = i[H, S] = -7 (Bo A S) 9{x) 9{L - x), 

S = i[H,S] = -7 [(Bo-S)Bo-Bo^S]0(x)0(L-x) 

- i7 (Bo A S) [p, {6{x) - d{L - x))]) . 

From the first of the above equations we obtain £{py2m) = -7 (Bo • S) [S{x) - S{L - x)] 
like using matching conditions at x = 0, L. The second and third equations give the spin 
evolution. By inspection we see that i) the spin remains parallel to Bq if it was initially 
so and, ii) the spin is constant in this case S = S = 0. Hence, in spite of choosing a 
case where the spin is conserved, its entanglement with the momentum decreases the spin 
correlations with the idle fermion. 

The effect of the magnetic field on particle A can be seen in its state. Behind the region 
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P, the resulting state vector as transformed from the one in Eq. (|7.11|) is 

^t(p)^f/[^T(p)]='^T(p) 
^i(p)^C/[*i(p)]=Ti(p) 




(7.14) 



where T^{p) (Ti{p)) is the transmission coefficient associated to the mesa (well) potential 
induced by Bq, for initial spin up (down). It is given by 



Tsip) :-- 



2pps cos{psL) - i{p^ + pI) s'm{psL) ' 



(7.15) 
for spin 



where Ps{p,Bq) := \J p^ — 2s'm'jBQ, as given by (j7.13p . Bq := |Bo| and s = ^(- 
t (I). As expected for Bq = 0, T^ip) = Ti{p) = 1. The initial state is preserved so no 
spin-momentum correlations are generated. In general, for Bq / 0, T|(p) ^ '^[(p), produc- 
ing spin- momentum entanglement. We are considering here just transmission through the 
region T>, without taking into account the reflection of the wave packets. We suppose all 
the measurements will take place beyond D so we may normalize the final transmitted state 
to 1. Finally, the net effect of this local interaction is the reshuffling of spin-momentum 
correlations in the state of the active fermion. Accordingly, the degree of spin-spin entan- 
glement between both particles decreases. As was done in Eq. (j7.4p . we evolve the state 
()7.10p with the transformation U and trace out the momenta 



Tr 



Pa,Pb 



i 5] ..'TrpjC/[^(^)(pa)]{C/[M/J)(pa)]}t) 



®TVp,(^L'i(Pb){*L°i,(Pb)r), 



(b) 



(7.16) 



where ss' :- 
states, 



normalized spin-spin state is 



The traces corresponding to particle B give just the initial spin 
because U is just the identity for B. The resulting, properly 



(7.17) 
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where 



:= J d^p\gip)\' 



\h{p)\' + \Tap)\ 
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(7.18) 
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Figure 7.3: Negativity N in Eq. (f7l9|) as a function of 7S0 for m = 100, = 10, L = 3, 
and cr^^) = 1, 2, and 3. The higher curves corresponds to the thinner u's. All quantities are 
measured with respect to a global arbitrary energy scale. 

The negativity for this state is found to be 

iV = 2|%|>0. (7.19) 

This expression for N is rather illuminating and its behavior easy to understand. For the 
initial state ()7.10p = 1 (1 initial ebit), and, as long as the magnetic field is increased, 
7|(p) and T^{p) become more different, making the term I-^j smaller, and diminishing A^. 
On the other hand, the wider a^^^ the more destructive interference between ^^(p) and T^^p) 
will occur, reducing N. We plot in Figure ES] the negativity as a function of 7-605 with 
Bo := |Bo|, for m = 100, p^^^ := |pf^| = 10, L = 3, and cj(^) = 1, 2, and 3. All quantities 
are measured with respect to a global arbitrary energy scale. The entanglement goes to zero 
with increasing Bq, and the wider o"^^\ the lesser the entanglement. A similar behaviour 
arises from the cumulative effect of the barrier; the larger L, the smaller the entanglement. 
We show in Figure [73] this behavior, plotting as a function of L for m = 100, p'^^ = 10, 
7S0 = 0.2, and <t('^) = 2. 

7.2.2 Two photons and an optically-active medium 

In this section we analyze a bipartite system, composed by two photons A and B, which are 
far apart in a polarization Bell state l^'") with factorized Gaussian momentum distributions. 
We consider that the photon A traverses a region where a finite, optically-active medium, 
exists. As a result, part of its spin correlations will be transferred to the momentum. 
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Figure 7.4: Negativity N in Eq. (fH^ as a function of L for m = 100, p''^^ = 10, 7^0 = 0.2, 
and 0"^'^) = 2. All quantities are measured with respect to a global arbitrary energy scale. 

Basically the mathematical formalism used for the two-fermion case is also valid here, 
with I and | indicating right-hand and left-hand circular polarizations, which we will denote 
by R and L indices. The transmission coefficient in the WKB approximation, at lowest order, 
takes now the form of a complex phase, depending on the polarization 



Ts{w) := exp[rans(u')L], s = R,L, 
where the refractive indices are 

nR,L{w) := \/l + Xii ±Xi2, 
and xui X12 are two of the matrix elements of the susceptibility 

X ■= -iXi2 Xu 
V X33 



(7.20) 



(7.21) 



(7.22) 



X is produced, for example, by an isotropic dielectric placed in a magnetic field Bq directed 
along z, which is also the direction of photon propagation. L is the dielectric length that 
the photon traverses, xii and xi2 are 



Xuiw) :-- 



Xi2iw) 



Me" 
meo 

meo 



w 



WWc 



(7.23) 
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where the cyclotron frequency Wc '■= e|Bo|/m, e is the electron charge, m its mass, wq the 
resonance frequency of the optically-active medium, M the number of electrons per unit 
volume, and eo the vacuum electric permittivity. 



mission coefficients. However, the approximation that considers these factors coincides ex- 
actly with the lowest order one when taking into account just linear terms in Bq. We will 
consider the realistic case in which Wc is small as compared to the photon average energy. 
Thus we will work from the very beginning just with the transmission coefficients ()7.20p . 
The negativity A^, obtained for this case in analogy with the two-fermion case, is 



where pq is the average momentum of photon A, a <^ po its momentum width, and B := 



We plot in Figure 1731 the negativity N in Eq. (I7.24p as a function of BL for po = 10, 
a = 2, and wq = 10. All quantities are measured with respect to a global arbitrary energy 
scale. The entanglement decreases as the magnetic field B or the length L of the dielectric 
increase. We plot also in Figure [7l6] the negativity as a function of a for pQ = 10, BL = 4, 
and Wq = 10. Surprisingly, and opposite to the two-fermion case, the entanglement increases 
as the momentum width a of the photon is larger. This effect comes from the fact that 
for larger widths, centered at wq, the contribution from the region around the resonance 
frequency wq, in which the effect of the medium is appreciable, becomes smaller. On the 
other hand, in the limit of negligible width, the spin could not get entangled with the 
momentum so in this limit the spin-spin entanglement would remain maximal. We observe 
then that there is a region of intermediate widths a in which the spin-spin entanglement 
becomes minimal. Finally, we plot in Figure 17.71 the negativity as a function of wq 
for Pq = 10, BL = 2, and a = 0.5,1,2. The larger a corresponds in this figure to the 
higher curves. These graphics show that the entanglement decreases mainly for resonance 
frequencies wq around the average momentum pQ. It also shows the surprising behavior 
mentioned above: For wider a, the entanglement is larger, and the interval of wq for which 
the entanglement decreases is wider, as expected according to the previous analysis. 



The next order correction would include factors ^tir^l in the denominators of the trans- 




,2 



(7.24) 



J\fe^\Bo\/{mho). 



7.3 Is purely quantum communication feasible? 



A cautious reader may immediately object that, in principle, our preceding analysis seems 
to suggest the feasibility of communication through a purely quantum channel, i.e. without 
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Figure 7.6: Negativity in Eq. ()7.24p as a function of a for po = 10, BL = 4, and wo = 10. 
All quantities are measured with respect to a global arbitrary energy scale. 
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Figure 7.7: Negativity in Eq. (j7.24p as a function of wq for po = 10, BL = 2, and 
a = 0.5, 1,2. The larger a corresponds to the higher curves. All quantities are measured 
with respect to a global arbitrary energy scale. 

classical communication, as in all quantum-informational protocols. Let us illustrate this 
point with the following proposal based in the previous two-fermion case. As usual, Alice 
and Bob will be the corresponding observers of each fermion. The joint spin-spin state 
is given by equation (j7.17p . which immediately drives one to the subsequent reduced spin 
state for fermion B: 



But the quantities Iss' depend on the magnetic field Bq {cf. equations (j7.18p and (|7.15p ). 
which can be controlled by Alice. This allows them to agree on the following procedure. They 
agree on communicating with a binary alphabet with classical bits and 1. If Alice were to 
communicate 0, she would adjust Bq so that the reduced spin state for Bob is, for example. 



They prepare a statistically significative amount of pairs of fermions under such condi- 
tions. Then Bob, when measuring the spin upon his fermion, will typically obtain spin up 
in the 75% of the measurements and spin down in the remaining 25%. He thus deduces 
that Alice is sending the bit 0. On the contrary, if Alice were to communicate the bit 1, she 




(7.25) 




(7.26) 
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would adjust Bq so that the reduced spin state for Bob is, for example, 




(7.27) 



They also prepare a statistically significative amount of pairs and Bob perfoms his 
measurements. He will detect 75% of them in the spin down state and 25% in the spin 
up state. He thus deduces that Alice is sending the bit 1. Notice that this information 
transmission is carried out without the assistance of classical communication. 

The flaw stems from the disregarding of the wave packet reflection in Alice's site. This 
can be seen in two complementary ways. On one hand, to perform a genuine information 
transmission in which Bob's fermion is actually carrying information encoded by Alice, he 
must be able to discern between those fermions whose pairs have been reflected in Alice's 
barrier potential, so that he can securely discard them (they are not carrying information 
at all). And this is only possible if Alice classically communicates this information to Bob. 

On the other hand, a detailed calculation taking into account the reflection coefficients, 
hence without the normalization appearing in the denominator of ()7.18p . shows that Bob's 
reduced spin state will be given by: 



where TZs{p) denotes the corresponding reflection coefficient for spin s and momentum 
p. The calculation reveals that Bob gains no information whatsoever from Alice's deci- 
sions, unless she classically informs Bob about them. Mathematically this can be expressed 
through the unitary character of the process. If no classical information is exchanged, 
the evolution is locally unitary (^f^") (g) vJ/C') [/(a)[\^(a)j ^ [/('') [^(f*)]) and thus cannot 
change the entanglement shared by both parties (the entanglement class is invariant un- 
der locally unitary evolution). Consequently no information through the purely quantum 
channel can be obtained. On the contrary, if Alice classically sends information to Bob, 
she is actually selecting a subset of her incoming fermions, i.e. she is projecting her state 
(\I/{«) (g) ^(b) (g, vl/W]^ where P^"-''^ is an orthogonal projectoi0), which is a 

^More generally, it can also be a POVM, depending on whether the information provided by Alice is 
complete or not [Per95) . 



PB 




(7.28) 
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nonunitary operator which changes the entanglement class. This fact allows them to ex- 
ploit the initial quantum correlations between their fermions to establish a communication 
protocol. This example shows once more the impossibility of using quantum correlations, 
i.e. entanglement to exchange information without the aid of classical communication. 

In summary, we showed the spin entanglement loss by transfer of correlations to the 
momentum of one of the particles, through a local spin-momentum entangling interaction. 
This phenomenon, already analyzed for a non-interacting particular case in the context of 
Wigner rotations of special relativity, may produce decoherence of Bell spin states. The 
momentum of each particle is a very simple reservoir and indeed it is one that cannot be 
eliminated by improving the experimental conditions, due to Heisenbcrg's principle. We 
show that an s = ^ fermion (photon), which initially belongs to a Bell spin state, may 
lose its spin correlations due to this physical phenomenon when traversing a local magnetic 
field (optically-active medium). These specific media entangle each component of the spin 
state of the particle with its momentum, like in a Stern-Gerlach device. This could have 
implications for quantum communication and information processing devices. 



Chapter 8 



Schmidt decomposition with 
discrete sets of orthonormal 
functions 

While multipartite entanglement and mixed states entanglement are resulting very difficult 
to study theoretically, entanglement of pure states of systems with just two components is 
much better known and very interesting for apphcations |Eke91[IBW92llBBC+93[lBBM+98[ 
lBPM+97j . The Schmidt decomposition [SchOGl [EK95] is a valuable tool for analyzing these 
states entanglement. This procedure (not to be confused with the Gram-Schmidt orthog- 
onalization procedure, which is a different one) relies on the singular value decomposition 
of a square matrix [NCOOj . It allows to express an arbitrary pure bipartite state as 'sum of 
diagonal biorthogonal products'. This way the entanglement of this state is easily evaluated 
and ready to use. In Appendix El we explain in detail the Schmidt decomposition procedure, 
both for finite dimensional systems and for those described by continuous variables. 

In this chapter we obtain an alternative way |LL05aj to compute the Schmidt decom- 
position for systems described with continuous variables, different to the usual one. It has 
the advantage of preserving the analytical dependence of the states, and it is well moti- 
vated due to the research community concern in entanglement. Here we consider the case 
of continuous variables entanglement. For us, these variables may be {a + a^,i(a^ — a)} 
which commute as phase space variables do. We also refer to continuous variable entan- 
glement in systems described by momentum and/or energy observables. Precisely, the 
entanglement of continuous variables stems from the original EPR article |EPR35j . The 
systems with continuous variables have been studied thoroughly both theoretically and ex- 
perimentally, |Vai94[ iFSB+gsl lLB99l iGieOTl IGECP031 IAB051 iRvLOS] and references therein. 
However, the treatment of the entanglement of systems with continuous variables, from a 
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practical point of view, is far from straightforward. Until now, obtaining the Schmidt de- 
composition in the continuous case required solving the corresponding integral equations 
[PBPOOl ILWEOOl ICLE021 ICLE03llLE04] . They had to be discretized, losing the continuous 
dependence of the initial state. Here we propose a method [LLOSaj to perform the Schmidt 
decomposition for this case, to the accuracy desired, keeping the continuous character of 
the variables. This method consists of two steps: 

1. We decompose the bipartite system wave function, f{p, q), by using two denumerable 
and complete sets of orthonormal functions, {On\p)}, {On\q)}, of L^, in the form: 



The purpose of this step is to transform the continuous problem into a discrete one 
(a necessary step for the numerical computation), while preserving the continuous 
dependence of f{p, q). 

2. Then we apply the (finite dimensional) Schmidt procedure to ()8.ip in order to write 
the wave function f{p, q) as diagonal sum of biorthogonal terms: 



The orthonormal functions 'i/'n {p) , (?) -the modes- will be some particular linear com- 
binations of On\p), On\q), respectively. Notice that we are using the Schmidt procedure 
for discrete systems to obtain the decomposition for the continuous case. This is much more 
tractable, as it implies diagonalizing matrices instead of solving integral equations. 

The rationale for this procedure is the expectation that only a few On will suffice: 
A handful of appropriate orthonormal functions will approximate f{p,q) to the desired 
accuracy. We finish by pointing out some properties of this method, namely 

• We obtain complete analytic characterization of the modes il^n\p), '>pn\q) to the 
desired precision. Our method surpasses the standard numerical procedures in that 
keeps the continuous features present in f{p,q). 

• We remark the portability of the attained modes ipn\p), ipn\q) that are ready for 
later uses. 

• For the physical systems analyzed, we found that the number of On functions required 
is small. For example, in the biphoton case analyzed in Chapter[9l the accuracy (error) 




(8.1) 



m.n 



n 
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was of around 2% with 26 x 26 Cmn matrices. Considering the case of two electrons 
which interact via QED, studied in detail in Chapter O the error was of 0.7% with 
12 X 12 matrices. 

We consider a bipartite quantum system formed by two subsystems and 82- Some 
examples are two photons entangled by parametric down-conversion, a photon emitted by 
an excited atom and as a result entangled with it or two charged particles which interact 
electrically. This system is described by the pure vector state 



where a|^^(p), a^2)('i) creation operators of a particle associated to the subsystems 

5*1 and 5*2. p and q are continuous variables associated to 5*1 and 5*2 respectively, which can 
represent momenta, energies, frequencies, or the like. In general, the analysis is made in 
an ad hoc kinematical situation in which p and q turn out to be one-dimensional variables, 
p £ (oi,6i), q £ (a2)&2)- In this chapter we assume this is the case. In addition, there can 
be discrete variables (like the spin) to be treated with the Schmidt method, that we do not 
include here to avoid unwieldy notation. 
Our method works as follows: 

We consider two denumerable and complete sets of orthonormal functions {On\p)}, 

(2) 

{On (q)} n = 0,1,..., 00, each one associated to each particular subsystem Sa (a = 1,2). 
These functions obey 



1. Our first step is to expand the wave function f{p,q) as a linear combination of the 
translating the continuous problem into a discretized one. Thus we work with 
the discrete coefficients of the linear combination, though the continuous character of 
the state is preserved in the k dependence of the oi"^ functions. The expansion reads: 




(8.2) 




(8.3) 



(8.4) 



n 



00 




(8.5) 



m,n=0 



where the coefficients Cmn are given by 




(8.6) 
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2. Our second step is to apply the Schmidt decomposition (see Appendix to the 
discretized pure bipartite state ()8.5p . as is usually done for finite dimension Hilbert 
spaces (diagonalizing matrices, instead of solving integral equations). In order to do 
this, it is necessary to truncate the expansion (jS.Sp . something that is possible to 
a certain accuracy due to the fact that J dpdq\f(j),q)\'^ < oo {f{p,q) is in principle 
normalizable) , and the expansion is in orthonormal functions, so the coefficients Cmn 
go to with increasing m,n (see below). 

We truncate the series (|8.5p at m = mo, n = uq, with mo < no, without loss of 
generality. The Schmidt procedure leads to 

mo 

/(p,'?)^Ev^^J'^(^')^f^(9)' (8.7) 

1=0 

where 

mp 

4\P) ^ Y.V,^Og\p), (8.8) 

m=0 

mo no 

4'\q) ^ Y.Y.V*^CmnO^^\q) (8.9) 

^■^^ m=On=0 

i = 0, ...,mo. 

Here the matrix V is the (transposed) eigenvectors matrix of Mij = M*^ = Xlnlo ^inC^n- 

mo 

^ MimVjm = XjVji, (8.10) 
m=0 

and {Aj}j=o,...,mo are the eigenvalues of M. 
There are two sources of error in this procedure: 

• Truncation error: This is the largest source of error in our method. Inescapably, 
the series (18. 5p must end at some finite m, n when attempting to obtain some 
specific result. This step is possible to a certain accuracy because the function 
/(p, q) is square-integrable and we are expanding it into orthonormal functions, so 
I]m=o S^o \Cmn? < OO and thus Cmn ^ wheu m, n ^ oo. 

The particular choice of the orthonormal functions O^") will affect how fast the Cmn 
go to zero. Hence, the election of these functions for a particular physical problem 
will be a delicate task. To reach the same accuracy with different sets {O^"-*} it will 
be necessary in general to consider a different pair of cut-offs {mo, no} for each of the 
sets. 



8.1 Localization point of view 
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• Numerical error: This is a better controlled source. It includes the error in calculating 
the coefficients Cmn via (j8.6p and the one produced when diagonalizing the matrix 
M = CCl 

The suitable quantity to control the convergence for a particular f{p,q) and a specific 
set {O^"^} is the well known (square) distance d]^^^^^ between the function f{p,q) and the 
Schmidt decomposition obtained with cut-offs {mQ,no} (mean square error): 

,1 _ c dpi': diifip^i) - T2U VK.i^^} {p)i^^:^\<i)? 

This expression gives the truncation error. It will go to zero with increasing cut-offs ac- 
cording to the specific {O^"^} chosen. 

Another easily computable, less precise way of controlling the convergence is given by 
the fact that (with no cut-offs) X],^=o ~ 1 1/(^5 9)11^ ^'^d thus 

— , - "1=0 ^"^ /Q 19^ 

is other measure of the truncation error, where here Am is calculated with cut-offs {?n.Oi ^o}- 
Would we compute the exactly, then d} = d?. In practice this can not be done because 
our Xn are the eigenvalues of the mo x mg matrix Mij, that depend slightly on mo, riQ. Both 
distances behave in a very similar way, as we show in FIG. 19.11 and FIG. 19.21 though d"^ is 
more easily computable than d^. 

The choice of the two sets of orthonormal functions for a particular physical problem, 
{0^'^^}a=i,2 can be approached from two different points of view, according to the feature 
one desires to emphasize: Localizability properties or convergence improving. 



8.1 Localization point of view 

The choice of the orthonormal functions in a particular problem can be done according to 
the specific intervals in which the variables p, q take values for that case. Typical examples of 
discrete sets of orthonormal functions are the orthogonal polynomials, defined in a variety 
of intervals. For example, a possible choice to describe one dimensional momenta p E 
(—00,00), are the Hermite polynomials, On\p) ~ Hn{p)- The equivalence sign indicates 
here that the polynomial must be accompanied by the square root of the weight function 
in order to be correctly orthonormalized, and normalization factors must be included. If, 
on the other hand, the variable of interest in a specific problem is bounded from below. 
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like the energy of a free massless particle p € (0, oo) , then the election could be Laguerre 
polynomials, On\p) ~ Ln{p). 

The criterion for choosing the orthonormal functions O^"^ according to the intervals 
in which p, q are defined has a fundamental character. For example, the localizability in 
configuration space of the Fourier transforms of the modes (|8.8|) . (|8.9|) . depends critically 
on the intervals in which these modes are defined |BB98l ILam02j . Only if we choose the 
functions 0^"^ to be defined exactly in the same intervals as the amplitude f{p, q), may the 
Fourier transforms of the modes have the right localization properties. In spite of that, this 
point of view may not be the most suitable one, as it may give slower convergence than the 
point of view presented below. 

8.2 Convergence point of view 

In this case, the choice is approached with the goal of improving the convergence. The 
are chosen here according to the functional form of f{p,q)- The closer the lowest 
modes are to / the lesser the number of them necessary to obtain the required accuracy. 
We are looking for O^"^ that maximize J^^^ dpOm* {p) dqOn^* {q)f{p, q) for low m, n. In 
some cases this practical point of view will be more useful than the fundamental one: the 
convergence will be faster. For example, suppose the amplitude for a particular problem is 
of the form f{p,q) = g{p, q)e~P^ ^^e~'^^ ^'^ , with g{p,q) a slowly varying function of p, q. In 
this particular case it is reasonable to choose the functions O*-"^ as Hermite polynomials, 
because their weight functions are indeed gaussians. This leads to Or?\k) oc Hn{k)e~^^ /'^ . 



Chapter 9 



Entanglement in Parametric 
Down- Conversion 



In this chapter we consider a realistic case of biphotons already studied in the literature 
[LWEOOl ILE04j : two photons entangled in frequency through parametric down-conversion. 
We apply our method to this physical system in order to obtain the Schmidt decomposition 
and the structure of modes without losing the analytic character within the target accuracy. 

The system under study is a biphoton state generated by parametric down-conversion 
(PDC) of an ultrashort pump pulse with type-II phase matching. The amplitude in this 
particular case takes the form [LWEOOj 



where u^o, We G (0, oo) are the frequencies associated to the ordinary and extraordinary 
fields respectively, k'^ and k'^ are the inverse of group velocities at the frequency a), k is the 
inverse group velocity at the pump frequency 2(D, L is the PDC crystal length and a is the 
width of the initial pulse. Typical values for these parameters are [k — k'^)L = 0.213 ps, 
{k - k'^)L = 0.061 ps, uj = 2700 ps~^ L = 0.8 mm and cj = 35 ps'^ 
We perform now the following change of variables: 



f{uJo,iOe) 



= exp[— (Wo + We — 2tj)^/o"^] 

X sinc{L[(a;o - uj){k'^ -k) + {uJe - ^){k'^ - fc)]/2} 



(9.1) 



UJn — UJ 



Lp = {k- k'^)La, 



(9.2) 



P = 



a 



Lq = (k — k'^)La, 



(9.3) 



a 



and thus obtain 



/(p,g)=e-(^+^)' 



smc[{Lpp + Lgq)/2]. 



(9.4) 
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Figure 9.1: dmo.no ^ function of the cut-offs {mo, no}. 

We liave applied our method to the function (j9.4p (once normaUzed) according to Chapter 
[51 in the following way: 

We choose as orthonormal functions Hermite polynomials, because their weights are 
gaussians and a Gaussian appears in (j9.4p . These polynomials were used in |UBW03j for 
PDC in some particular cases which are exactly solvable. The orthonormal sets we chose, 
looking for maximizing the Cmn (|8.6p for the lowest m, n, were 

OW(fc) = (^2"n!)-i/2^„(A;)e-^'/2 ^ = 1,2. (9.5) 

This choice of polynomials is suitable for the convergence approach (section 18. 2p . taking 
into account that (D 3> cr and thus the interval of definition of /(wo, ^e) can be restricted to 
a region centered in Co of width ~ a in Uq, Wg. 

We have considered cut-offs mo = no taking values {5 — 25} and followed the steps of 
Chapter [8l We have computed the eigenvalues A„ of the Schmidt decomposition (j8.7p for 
each pair {mo, no}. We have also computed the modes (18. 8p and (18. 9p . 

In FIG. 19.11 we plot the distance (j8.1ip as a function of mo = no, to show how 

fast the convergence is. With mo = no = 25 the truncation error is of 2%. We also plot in 
FIG. E2] the distance d^^ ^^ ^KTIh . which serves as another measure of the convergence, as 
a function of mo = no. We obtained 25 = 2%. 

Regarding now the most precise case considered, mo = no = 25, we plot in FIG. 19.31 the 
eigenvalues A„ for different values of n, observing good agreement with the results existing 
in the literature |LWE00j . For this case we also plot in FIG. [93] the modes ((SSl) and ([5:9]) 
for i = 0, 1, 2, 3, confirming the validity of the method when comparing with [LWEOOj . 



5 10 15 20 25 

Fi gure 9.2: c/^^^^p as a function of the cut-offs {mo,?^o}- 




2 4 6 3 



Figure 9.3: Eigenvalues A„ versus index n. 
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Figure 9.4: Modes ^p^n'{p), ^pn'{q) as a function of j3 = and g = forn = 0,1,2,3. 
The modes are given explicitly by: 

25 

^l^Hk) = e-''/'^{V^2^n\)-'/'A<;^lHn{k) (9.6) 

n=0 

m = 0,...,25 a = l,2, 

where the values of the coefficients Ami are obtained through (j8.8|) and (|8.9|) . The actual 
properties of the modes ()9.6I) depend on these values. In fact, the parity and number of 
nodes is determined by them, taking into account that Hn is a polynomial of degree n, 
parity (— )" and having n nodes. 

A good approximation to the V'o^^ (p) obtained with our procedure is 

^W(p) = e-^''/2 (0.81395 - 0.1476V + 0.00821/). (9.7) 

This expression has a deviation (square distance) of 10~^ from the whole mode obtained 
including terms until p^^, which is the greatest power appearing for rriQ = uq = 25. On the 
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'mo.np 



mo = no 


P = 0.5 


p= 1.0 


/? = 2.0 


25 


0.13 


0.020 


0.037 


20 


0.19 


0.024 


0.041 


15 


0.27 


0.032 


0.050 


10 


0.38 


0.062 


0.064 



Table 9.1: for p = 0.5, 1.0, 2.0 and mo = no = 25, 20, 15, 10. 

other hand, 44 - 45 25 = 0-213 > 10"^. From ([HZD it can be seen that in this mode the 
even components are greater than the odd ones (these are neghgible), so it is an even state, 
as shown in FIG. 19.41 

Another example is the approximation to tJj\^ {q)'- 

^pPiq) = e"'?'/2(2.91088g - 3.54070g^ + 1.29062g^ 
-0.204029"^ + 0.01598g^ - 0.00063g^^ + O.OOOOlg^^). (9.8) 

This has a deviation (square distance) of 10~^ from the whole mode obtained including 
terms until q"^^. On the other hand, ^3 — 45 25 ~ 0.020 ^ 10~^. More terms are needed 
in (j9.8p . because they go to zero more slowly with increasing powers of q. Here the most 
important components are the odd ones (the even ones are negligible), leading to an odd 
parity state, as shown in FIG. 19.41 

To show how the convergence of the method depends on the specific family pairs of 
orthonormal functions {On\p)}, {On\q)} chosen, we consider the cases of Hermite or- 
thonormal functions depending on a parameter /3 related to the width of the Gaussian, 
fixed for each family pair: 

0(")(A.) = ^^^/7„(/?A:)e-(WV2 ^ = 1^2. (9.9) 
V/7r2"n! 

We applied our method to the amplitude (j9.4p with these sets of orthonormal functions, 
for P = 0.5, 1.0, 2.0, and cut-offs mo = no = 25, 20, 15, 10. We show in table 19.11 the values 
of „g for these specific parameters. 

Clearly, the convergence is better for the case P = 1.0, which we used in the preceding 
calculations. In case we chose another type of orthonormal function for (j9.4p (Laguerre, 
Legendre,...), the convergence would have been much worse because of the specific shape of 
that amplitude. 



Chapter 10 



Momentum entanglement in 
unstable systems 



10.1 Time evolution of bipartite entanglement 

Writing about time in scattering theory is rather tricky, because the in and out states are 
asymptotic states which come from and leave to the spatial infinite, i. e., very large distances 
compared to the interaction distance. Thus the time between emission and detection of the 
particles is very large compared to the interaction period. Interesting situations arise when 
the time lapse between preparation and detection is finite. Some examples are entangled 
systems which fly apart from each other but still lie in the same region, or the effect of the 
interaction between particles during a finite time interval. In these cases, it is necessary to 
use a finite time evolution formalism (see Chapter . This implies that the time integral 
of the action is restricted to a finite interval [— Associated to this finiteness of the 
integration interval for the time coordinate is the non-conservation of energy at short times. 
The point is that, instead of a Dirac delta 5{AE), a function of the kind 



is obtained when performing the integral in t, where AE is the difference between the final 
and initial energies. The interpretation of this is that the finite time evolution allows non- 
conservation of energy for short elapsed times compared with the interaction time interval. 
In the lavge-t limit, the sine function (110. ip converges to the 6 and thus the usual scattering 
expression is recovered. The situation changes when we initially have an unstable system, 
like an excited atom or a resonance. In this case, the system has a certain width T wherever 
inelastic decay channels are open (see below). For these systems, there is an additional 
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term —Ft in one of tlie exponentials of the time-dependent amplitude, so this exponential 
goes to with increasing t. The effect is that, at large times, the remaining amplitude of 
decay from the initial state \i) to the final state |/) is just a Lorentzian curve |CLE02| with 
a global time-dependent phase 

lim(/|[/(t,-t)|.)oc^J|i^. (10.2) 

In this chapter we study |LL05b] the V dependence of the time evolution of momentum 
entanglement for unstable systems. In the F = case, as we have indicated, a function 
(jlO.ip appears in the transition amplitude between the initial and final states, for a finite 
time lapse. The entanglement in energies increases as this function goes to (5(A£^) with 
increasing time, because the 5-function is the (non-normalizable) state with highest (infinite) 
entanglement (see Chapter II ip. On the other hand, if T 7^ 0, the entanglement in energies of 
the decay amplitude increases with t until it reaches a maximum. This is the entanglement 
of (jl0.2p . which remains constant in the subsequent time evolution. Notice that the smaller 
the value of T, the larger this constant value is and the later it is reached, as we will explicitly 
show in this chapter. 

The elastic (non decay) amplitude has been extensively studied in the literature [Kha57t 
ICSM77j . These analysis show the presence of power corrections at long times. Here we 
consider the decay amplitude, obtaining power corrections in t that introduce entanglement, 
as they are not factorizable. This is of great interest as these corrections may be sizeable 
in some relevant cases |JMS+05| . 

We consider a non-elementary unstable system A in an excited state |e) with mass me 
which decays into a ground state \g) with mass vfig emitting a particle 7 with mass (see 
FIG. 110. ip . This system could be an excited atom which emits a photon and gets down 
to the ground state, an unstable nucleus which radiates, or the like. In the atom example 
me — TUg is the atom mass, me — rUg = ujq, where loq denotes the energy difference of the 
internal atom levels |e) and \g), and m^ = 0. For simplicity we consider that the emitted 
particle is a scalar. Our aim is to study the evolution with t of the entanglement in the 
final state momenta for different values of F. The Hamiltonian we consider H = Hq + Hj is 
customized for analyzing the process e gj. Hq is the free Hamiltonian of system A and 
particle 7. Hj is the interaction Hamiltonian. 

Ho = {Te+u;o)\e){e\+Tg\g){g\+ J d^p^T^al_ap^, (10.3) 

Hi = A(|57)(e| + |e)(<?7l)- (10.4) 

Te, Tg are the kinetic energies of system A in states |e) and \g). is the kinetic energy of 
particle 7. cuo is the energy difference between levels e and g, including the mass differences 
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1 j/ 















Figure 10.1: Decay process of an unstable system which emits a particle 7. 

rUe — rUg — m^. A is the coupling constant between the states |e) |0) {A in state |e) and 
no particles 7) and \g) (8> I7) {A in state \g) and 1 particle 7), that are the only ones that 
are taken into account in the present approximation. 

We denote in the following Ee = Te + loq, Eg = Tg and w-y = T-y. |pe, 0) is the state of 
system A with momentum Pe and internal state |e), with no particles 7. |pg, p^) is the state 
of system A with momentum pg and internal state \g) and one particle 7 with momentum 
P7 

The amplitude for the decay e — > 57 during the time interval [—t, t] can be given 
|CTDKG92j as 

Sfilt = exp[i{Eg + uj.f + Ee)t]{pg,p^\Uit,-t)\pe,0), (10.5) 

where the time evolution operator U can be obtained integrating the resolvent operator 
G{z) := l/{z - H) along C+ (see FIG. [102]), i.e. 

(p3,p^|C/(t,-t)|Pe,0) = / dzexp(-2izf)(p3,p^|G(z)|pe,0). (10.6) 

2vrz Jc+ 

For the case we are interested in, where the 7 particle in the final state behaves as a 
mere spectator, it is useful to split the space of states by means of the projectors 



P = |e)(e| + |^^)(g|, 
Q = I -P. 



(10.7) 
(10.8) 
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Q includes states of the form \g){g\ |7)(7|- We consider only transitions between states 
|e,0) and \g,-i) induced by Hi. Thus QHQ = QHqQ and PHjP = 0. 
After some algebra it follows 

QG(z)P = QHiP -— , (10.9) 

^ ^ ^ z-QHQ^ z- PHqP - PR{z)P' ^ ' 

where 

R{z) = Hi + HiQ I-—QH1. (10.10) 

z - QHQ 

Near the real axis 

(e,0|i?(£;±ie)|e,0) = l^e{E) ^ iT{E) /2, (10.11) 

where Ae(-E) and r(£^) are respectively the energy shift and the width of the excited state 
|e), given by 

Ae(i?) = p(e,0|g/Q ^_^^^^ Qff/|e,0), (10.12) 
T{E) = 27r{e,0\HiQ6{E -QHoQ)QHi\e,0). (10.13) 

Due to (jlO.lip G{z) needs two Riemann sheets because there is a discontinuity in R for 
E > rUg + m^, where T{E) does not vanish. 
The resolvent matrix element is then 

(p,.p,|G(i;±i.)|p.,0) = ^E-E,-E,±U)lEL^-AAE)±.HE)r <1"'"' 

where h{' = (p,,, p^|iJ/|pe, 0). 

With all this, integrating along the contour of FIG. 110. 2|, by Cauchy theorem, (jlO.Sp 
transforms into 

i{Ee-Eg-LU^)t _ -i{Ee-Eg-w^-iV)t 

= E,-E,-.,-^T|2 + 

where we consider F as a constant parameter and we neglect the energy shift Ae(-E'). The 
integral along R vanishes when the contour is taken to infinity. Icut is the integral along 
path L which lies on the first and second Riemann sheets. It gives the power-law corrections 
in t for large times. It was discussed in [Kha57t [CSM77j . For instance, for F/£'e <C 1, it is 

= -^r—^-^^^^{a(/>(a,2t)-a(/<(6,2t)- 6-1/2/4 

X (6-a)[2\/2^-e2^*7r(46t- l)erfc(\/26t)]}, (10.16) 
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where a = —i{Eg + uj^), b = —iEf, and 

0(x,t) := -^e'^^erfcf^/^). (10.17) 
2 ^Jx 

The dependence of (jl0.16p on the final energies is different from that of the first term on 
the rhs of (|10.15p . and thus the entanglement in final momenta is also different. Eberly et 
al. |CLE02t ICLE03j analyzed the dominant part of this term showing the leading role of 
the Lorentzian in the features of the final state entanglement. Here we want to connect 
these results, valid for unstable systems, to the case where the initial state is more and 
more stable. The contribution from (jl0.16p gets lesser as T gets thinner. Thence we firstly 
neglect it, analyzing the entanglement associated to the amplitude F{t, T) := Sfi\t — Hplcut- 
We will justify this approximation later for the particular cases analyzed. The behavior of 




Figure 10.2: Integration contour in the z complex plane for evaluating (110. 6p with residues. 



F{t, r) for critical values of t and T is 



lim F(t,T) 
r^o 



hm F{t,T) 
r^oo 

limF(t,r) 
lim F{t,r) 



-2iH 



-iH 



fi Sm\{Ee - Eg- LO^)t\ 

Ee — Eg — LO^ ' 

fiex.^[i{Ee - Eg- u:^)t\ 



r/2 

/j, exp[i(£'e - Eg- uj^)t\ 
Ee — Eg — oj^ — iV /2 



-H 



(10.18) 

(10.19) 
(10.20) 
(10.21) 



(jlO.lOp and (|10.20p correspond to amplitudes with no entanglement in the final energies Eg 
and uj^y: they are separable amplitudes. (jlO.lSp contains a sm{AEt) / AE function leading 
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asymptotically to d{AE), already discussed |LL05aj (see Chapter [TT|) . Finally, (jlU.2ip 
produces the Lorentzian entanglement in Eg and uj^ multiplied by a separable phase. The 
entanglement of (jl0.2ip does not evolve with time, because the separable phase can be 
factorized. Summarizing: For F = the entanglement increases monotonically with time 
towards the maximum entanglement possible, associated to the <5-function. For F 7^ 0, the 
entanglement grows with time until it reaches a maximum, which is the corresponding value 
for the Lorentzian (jl0.2ip |CLE02j . The greater the value of t for fixed F, the greater the 
entanglement. The wider the value of F for fixed t, the lesser the entanglement. The fact 
that the asymptotic attainable entanglement decreases for increasing F is rather surprising 
at first sight. A wider width F is associated to a stronger interaction, that would in principle 
generate more entanglement. We interpret this result in the following way. A wider width 
F has associated a shorter mean life so that the particles reach faster the asymptotic limit. 
Accordingly, they interact during a shorter period and generate less entanglement. This 
is related to the fact that the exponential term that goes to zero, exp(— Ft), decreases 
exponentially faster for a linear growth in F. 



10.2 A specific example 



To explicitly illustrate the reasonings of section 110.11 with a specific case, we consider an 
initial state |e) for system A with Gaussian distribution in momentum pe, centered in 



(Pe|e) oc e-P'/"' 



(10.22) 



Using (fT0:il) we have Ilj — X6^'^^(jpg + p^y — Pe)- Then, we obtain the amplitude f{Pg, Pj} 
for the decay of the state |e) into a state \pg, p^) in a time 2t computing the relevant matrix 
element (|1U.15I) . 

For the case where A (in state \g)) and 7 are detected in opposite directions from the 
initial position of system A (in state |e)) we get 



/(p,g)oce-(P-'')'/^' 



iAi 



-i(A-ir)i 



A - iF/2 



cut 



(10.23) 



where 



Pg P-f 

p = — ; q = —, 

rup nip 



A = ^/{p- qY + 1 - \ jp^ + ml - J + 



(10.24) 
(10.25) 
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Above we have used dimensionless parameters o obtained by multiplying or dividing by mg, 
for instance t = met. 

We obtained the Schmidt decomposition of the bipartite amplitude (|10.23p . according 
to Chapter [HI for rhg = 0.7 and rh^ = 0.1. We observed the contribution from Icut was 
always lesser than 10"'^ times of the total amplitude (jl0.23p (for the considered cases). 

In FIG. we plot the Schmidt number K = (E^=o |GRE94j versus i for f =0, 

0.005, 0.015, and 0.03. K is a measure of the entanglement of a pure bipartite state, that 
gives the number of effective terms in its Schmidt decomposition. K > 1, K = 1 for 
separable states, and, the larger is K, the larger the entanglement. The graphic was made 
by computing the Schmidt decomposition for five values of t and four values of T (with an 
error lesser than (P = 1%, see Chapter [SJ in all the points but two, and lesser than d? = 5% 
in these two points, which were harder to compute; however, we remark that the attained 
precision is good for our present purposes, which consist of showing a general behavior of 
X as a function of T and t) and afterwards a polinomial interpolation for obtaining the 
continuous figure. As we have shown in this chapter, the entanglement grows with t in the 
beginning, and keeps growing for F = while for F 7^ it reaches a maximum. For each 
t, the wider the value of F, the lesser the entanglement. This is a surprising result because 
F is wider for stronger interactions, which presumably would create more entanglement. 
However, the stronger the interaction, the faster the asymptotic limit is reached. As a 
result, the system saturates before at a lower degree of entanglement, as we show here. 

In FIG. [Ton we plot the Schmidt number K versus f for i =50, 100, 150, and 200. 
The graphic was made by using the computed points of the previous figure (see above) and 
afterwards a polinomial interpolation for obtaining the continuous figure. The entanglement 
decreases with F. For each F, the longer the time t, the greater the entanglement. 



10.3 Remarks 



For the sake of completeness we point out that the term (jl0.16p may be important when F 
is not so small, and there is a time in which this term becomes larger than the exp(— Ft) 
term. This value of t is 



to 



ilog 



8V27r 



Ee ^ {Eg + Ey) 



{Ee — Eg 



E^y + f 2/4 



(10.26) 



The value of F for which the term (jl0.16p becomes larger than the other terms (exp(— Ft) 
term and Lorentzian term) for the case here considered of rhg = 0.7, = 0.1 and for 
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l.i 



1.0 




100 

i 



300 



Figure 10.3: Schmidt number K versus t for T 
curves correspond to the thinner F's. 



0, 0.005, 0.015 and 0.03. The higher 



t = 50 is around T = 1, i.e., of me order, which does not seem to be a reahstic value in 



principle. However, some recent results [JMS"'"05] point to the possibility of decay in S-wave 
in such a way that the width is larger than the difference in energies between excited and 
ground states. In such the contribution from Icut would be relevant. 

In summary, we analyzed the generation of momentum entanglement in the decay of un- 
stable systems described by a decay width F. We verify that, as expected, the entanglement 
grows with time until reaching an asymptotic maximum. On the other hand, unexpectedly, 
the wider the decay width F, the lesser the asymptotic attainable entanglement. We explain 
this apparently surprising result in terms of the fact that for wider width the mean life is 
shorter, so that the system evolves faster (during a shorter period) and can reach lesser 
entanglement than with longer mean lives. 
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Figure 10.4: Schmidt number K versus T for t = 50, 100, 150 and 200. The higher curves 
correspond to the longer t's. 



Chapter 11 



Maximum entanglement: The 
Dirac delta 

One interesting (non-normalizable) amplitude when studying continuous variable entan- 
glement is the Dirac delta. In this chapter we give |LL05aj a concise description of its 
entanglement content based on the method we proposed in Chapter [8] for computing the 
Schmidt decomposition of an amplitude with continuous variables. 

In this case we have (see Chapter [8]) f{p,q) = S{p — q) and we take the same interval 
(a, b) for p and q to apply our method, for analyzing its entanglement. We consider complete 
sets of orthonormal functions satisfying On\k) = On^*{k). A particular case is when they 
are real functions, as for example the typical orthogonal polynomials (Legendre, Hermite, 
Laguerre, Chebyshev,...) are. We must take into account that the Dirac delta is not 
a function but a distribution, and indeed is outside L?'. However, we can calculate the 
Cmn and study how much entanglement does this state have. We obtain straightforwardly 
Cmn = ^mn- This gives 

oo 

5{p-q) = Y,0^:^>{p)0^:^\q), (11.1) 

n=0 

which is just the resolution of the identity as given in (j8.4p . The Schmidt decomposition 
of the Dirac delta is not unique, because all the weights are equal to one (they are 
degenerate). In fact, the decomposition can be done with any complete, denumerable set 
of orthonormal functions, in the form ()8.4p . This expression can be seen as an infinite 
entanglement case, in the sense explained below. The fact that all the weights are equal to 
one, makes sense only because we are considering a distribution, not an LP' state. The sum 
of the squares of the weights, which must be equal to the square of the norm of the function 
f{p,q), diverges because the Dirac delta is not square-integrable. 

A possible measure of the entanglement of a state f{p, q) in its Schmidt decomposition 
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(|A.14p is given by the von Neumann entropy [NCOOj 



oo 




(11.2) 



Tliis is usually called the entropy of entanglement. 

The state of closer to (jll.ip is the case of an entangled state with N diagonal terms 
with equal A„, when goes to infinity. To be correctly normalized it verifies = 1/N, 
n = 0,...,N - 1 and 



This is the maximum entanglement case. This provides an estimate of the entropy of the 
Dirac delta (were it in L^). Another well-known way of obtaining the amount of entan- 
glement of the Dirac delta is considering it as the limit of infinite squeezing of a bipartite 
Gaussian state [GieOlj . However, to our knowledge, the Schmidt decomposition of the Dirac 
delta has never been obtained before. We believe ours is an elegant result, that may clarify 
the entanglement structure of this distribution. 




N-l 



(11.3) 



ra=0 



Part III 



Multipartite entanglement 



Chapter 12 

Sequential quantum cloning 



12.1 Quantum cloning sequentially implemented 

Multipartite entangled states stand up as the most versatile and powerful tool to perform 
information-processing protocols in quantum information science |BDOO| . They arise as an 
invaluable resource in tasks such as quantum computation |DE981 IRBOlj , quantum state 
teleportation [BEAOO] . quantum communication [ABH"'"Ol] and dense coding [BW92j . As 
a result, the controllable generation of these states becomes a crucial issue in the quest for 
quantum-informational proposals. However, the generation of multipartite entangled states 
through single global unitary operations is, in general, an extremely difficult experimental 
task. In this sense, the sequential generation studied by Schon et al. jSSV"'"05| ISHW"'"06] . 
where at each step one qubit is allowed to interact with an ancilla, appears as the most 
promising avenue. The essence of this sequential scheme is the successive interaction of 
each qubit initialized in the standard state |0) with an ancilla of a suitable dimension D to 
generate the desired multiqubit state. In the last step, the qubit-ancilla interaction is chosen 
so as to decouple the final multiqubit entangled state from the auxiliary D-dimensional 
system, yielding [SSV"'"05) 

l^)= E {^F\Vl:yVi^\^i)\h---in), (12.1) 
ii---j„=0,l 

where the Vj^^ are D— dimensional matrices arising from the isometries Vj^j : TIa ® {|0)} 
'Ha^'Hb^, with Ha = <C^ and Hb^, = being the Hilbert spaces for the ancilla and 
the kth qubit, respectively, and where \(pj) and \(Pf) denote the initial and final states 
of the ancilla, respectively. The state (|12.ip is, indeed, a Matrix-Product State (MPS) 
(c/. [EckOSl IPGVWC06| and multiple references therein), already present in spin chains 
|AKLT87] , classical simulations of quantum entangled systems |Vid03j and density-matrix 
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renormalization group techniques |VPC04j . Moreover, it was proven that any multiqubit 
MPS can be sequentially generated using the recipe of Ref. [SSV+05| . Notice that in this 
formalism, the mutual qubit-ancilla interaction in each step k completely determines the 
matrices Vj^*, ik = 0, 1, whereas we enjoy some freedom to build such an interaction from 
a known V^^ . This freedom stems from the very fact that in the proposed scheme only the 
initial state |0) for each qubit will be relevant. 

Here we consider [DLL^OG] the possibility of implementing quantum cloning based on a 
sequential protocol with the help of an ancillary system. This problem is certainly far from 
being an application of Ref. [SSV"'"05] . given that the initial and final states are unknown. In 
this sense, any proposed strategy will be closer to the open problem of which global unitary 
operations (certainly not all of them) can be implemented through a sequential procedure. 
Despite the fundamental no-cloning theorem |WZ821 IDie82j . stating the impossibility to 
exactly clone an unknown quantum state, there exist several cloning techniques with a 
given optimal fidelity [SIGAOSj (see Appendix [B]) . These procedures differ either from the 
initial set of states to be cloned or from symmetry considerations. In general, an optimality 
condition of the cloning procedure is obtained via the maximization of the fidelity between 
the original qubit and each final clone state. We will show how to perform sequentially 
both the universal symmetric |BH961 IGM97j and the economical phase-covariant symmetric 
quantum cloning [DM031 IBDM05] from one qubit to M clones. In the first case, a global 
unitary evolution transforms any input state {ip) in a set of M clones whose individual 
reduced states pout carry maximal fidelity with respect to IV'): i^i,M = ^f^j^- This cloning 
procedure is fully described by the evolution 

1^)0 15)^ Y^a,\{M-j)i^,ji^^)s^\{M-j-l)r,jr^)s, (12.2) 

j=Q 



where \B) denotes the initial blank state, aj = y m{m+i) ~ j)(t^:j4^'^)s denotes 

the normalized completely symmetric state with (M — j) qubits in state (p and j qubits in 
state (p'^ . As a relevant feature it must be noticed that the presence of M — 1 additional 
so-called anticlones is necessary in order to perform this cloning procedure with the optimal 
fidelity. The anticlone state ip* refers to the fact that they transform under rotations 
as the complex conjugate representation. For concreteness sake we have chosen \tp*) = 
cos^/2|l) + e"**^ sin 0/2 1 0) in coincidence with the seminal paper by Buzek and Hillery 
|BH96) . where \'ip) = cos0/2|O) + e*'^sin0/2|l). In the second case, motivated by quantum 
cryptoanalysis, the goal is to clone only those states belonging to the equatorial plane of 
the Bloch sphere, i.e. those such that 6 = tx 12. Furthermore, we have only focused upon the 
cases where no anticlones are needed (hence the term economical). Under this assumption, 



12.1 Quantum cloning sequentially implemented 



125 



imposing the purity of the joint state, the number of clones M must be odd [DM03]. The 
cloning evolution is now given by 

1 



\{k + l)0,kl)s + e'^\kO,{k + l)l)s , (12.3) 



where k = (M — l)/2 and where we have followed the same convention as above. In order 
to employ the sequential ancilla-qubit device as a quantum cloning machine we will firstly 
elucidate the minimal dimension required for the ancilla. The basic idea is to express the 
final states (|12.2p and (112. 3p in its MPS form, following Vidal's recipe |Vid03j (see Appendix 
ED: 

I^)=E( E r[i]L\A[iur[2]j.^^„^A[2],,r[3]^3^„3...r[n]LvJ ^ (12.4) 

ii...in \ai...a„-i j 

We identify the matrices Vj^^ by matching indices in expressions (jl2.ip and (jl2.4p . The 
indices Uj run from 1 to X) where x = ™cix-p{x-p}, Xv denoting the rank of the reduced 
density matrix p-p for the bipartite partition V of the composite system |Vid03) . To clone an 
arbitrary input qubit state |^) = a|0)+/?|l), we exploit linearity and determine the minimal 
dimension Di of the ancilla to perform the cloning for the state |0) and then similarly for 
the state |1). Then we combine both results in a single ancilla to obtain its desired minimal 
dimension D. Let us focus upon the symmetric universal cloning of |0). To determine 
the minimal dimension Dq of the ancilla we need to compute which can be undertaken 
without the exact MPS expression for the state. We need some previous results: 

Proposition 1. Let |^) and \(^) he multipartite states of the same system related through 
an invertible local operator F-p = Fa Fb for the partition V = A\B: 

\^) = Fp\^). (12.5) 

Then xv{4') = Xvi'4')- 

Proof. Recalling that the rank of p-p coincides with the rank of the coefficient matrix Cp of 
the corresponding state for that partition V, the application of the invertible local operator 
Fp amounts to changing the local basis of each part A and B in which the coefficient 
matrix is expressed. Since the rank is invariant under local changes of basis, we will have 
r (C-p((/>)) = r (C-p(-0)) for the bipartite partition V. Hence xvi'P) = Xv{'^)- D 
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Proposition 2. Let C^|j\,f_i be the coefficient matrix of the state p2.2p for the partition 
M\M-1. Then 

^{Cm\m-i) = M. (12.6) 

Proof. In virtue of the preceding proposition and using the invertible local operator Sm ® 
•Sm-i, where Sk denotes the normalized symmetrizing operator for K qubits, we only need 
to compute the rank of the coefficient matrix of 

M-l 

aj\{M - j)0,jl) (g) \iM - j - 1)1, jO), (12.7) 

i=o 

where the states are no longer completely symmetrized. Given the orthonormality of the 
involved states and the number of different components (M), it is clear that there are only 
M different rows, whereas the rest are all null, i.e. i'(CAfiM-i) = D 
Finally we prove the following 

Proposition 3. Let Cjt|2M-i-fc be the coefficient matrix of state (jl2.2p for the partition 
k\2M -k-1, where A; = 1, 2, . . . , 2M - 2. Then 

^{Ck\2M^i-k) < ^{Cm\m-i) (12.8) 

Proof. The key point is to realize that the matrices C^2M-i-k obtained from Cm\m-i 
by appropriately adjoining rows or columns to make them longer. From the preceding proof 
it is clear that there are only M different rows in C^|jvf_i, the rest being all null, thus the 
reordering procedure to build the other matrices cannot increase the former rank. Hence 
the stated result. □ 

With all these propositions it elementarily follows that x = that the minimal 

dimension Dq to clone the |0) state is Dq = M, namely the number of clones to produce. 
Repeating the same argument for the initial state |1) we also conclude that the minimal 
dimension of the ancilla to clone the |1) state is Di = M, as expected. Now we must 
combine both results to find D for an arbitrary unknown state {ip) = a\0) + /3|1). It is a 
wrong guessing to think that it should also he D = M and, consequently, a different scheme 
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must be given. The MPS expression of (|12.2|) for the original state |0) determines the 
/^-dimensional matrices V^^i^y whereas the corresponding MPS expression for the original 
state |1) determines ^/[^p 



\GMMiO)) 
\GMm{1)) 



E (^?Voi;i---^oyi]io)^ 



III ■■■ir. 



ii.. .1,1=0,1 



E {4''\y!i:^^■■y!(^\o)D\^l...^n). 



ii...i„=0,l 



Here, and \tp^p) can be calculated explicitly and will play an important role below. 

We propose now to double the dimension of the ancilla, C'^ C"^ C^, in order to 
implement a deterministic protocol of sequential quantum cloning. 



Protocol 1. Let V^j^^ = |0)(0| "l^J*^, + V'^^j^y Then 



o[k] 



i. Encode the unknown state in the initial ancilla state \ipi) 



|0)d. 



ii. Allow each qubit k to interact with the ancilla according to the preceding 2D -dimensional 
matrices Vj^'j . 

Hi. Perform a generalized Hadamard transformation upon the ancilla: 



i=[|0)®|4°)) + |l)®|(^«)^ 



;i2.10a) 



Note that the choice — > (based on pedagogical reasons) could be changed, 
equivalently, to . In this way, Eq. ^2.10a\) would not display entangled 

states but simple linear superpositions. 

iv. Perform a measurement upon the ancilla in the local basis {|0) \'f^p^), |1) W^p)}- 

V. If the result is |0) Iv^f^) (which happens with probability 1/2), the qubits are already 



in the desired state; if the result is |1) (probability 1/2), perform a local 



IT 



phase gate upon each qubit, then they will end up in the desired state. 
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Proof. After the first two steps, the joint state of the ancilla and the qubits is a y\0) (8) \Vp 
\GMm{0)) + f3 \GMm{1)), where originally = a\0) + /3\1). After the 

Hadamard rotation (step iii), this state becomes 

(|0) [q|GMm(0)) +/?|GMm(1))] + 
+i= [a|GMM(0)) -/3|GMm(1))] 

The remaining steps follow immediately from this expresion and from linearity |GM97j . 

□ 

Notice that despite the measurement process in step (iv), the final desired state is 
obtained with probability 1. In summary, the minimal dimension D of the ancilla for 
cloning M qubits is -D = 2 x M, i.e. it grows linearly with the number of clones even although 
the dimension of their space grows exponentially (2^''^). The reader can straightforwardly 
convince himself that were we to clone a d-dimensional system, the minimal dimension for 
the ancilla would he D = dx M, with an obvious generalization of the preceding argument. 
For the symmetric phase-covariant cloning, the same arguments can be reproduced provided 
we realize that the first term of the r.h.s. of Eq. (jl2.3p can adopt a similar form to (jl2.2p : 

k 

\{k + 1)0, kl)s = J2^j\{k + l-j)0,jl)s ^\{k- j)l,jO)s, (12.11) 

j=0 

where jj ^ for all j, and similarly for the second term. Thus for symmetric phase-covariant 
cloning the minimal dimension for the ancilla is D = 2 x (k + 1) = 2 x ^^^^ = M + 1. 
We see that the dimension of the ancilla D also grows linearly with the number of clones, 
although it is now lesser than above. This is a direct consequence of reducing the set of 
possible original states to clone. 



12.2 A specific illustrative case: 1^3 

In this section we explicitly compute the isometries for the universal symmetric quantum 
cloning from 1 qubit to 3 clones. We will show how the dimension of the ancilla is indeed 
2x3 = 6. We will be mainly using the formalism for computing the MPS form of an 
arbitrary multipartite state by Vidal |Vid03] . The final state after the cloning procedure 
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for the 11)0, according to Gisin and Massar [GM97] is 



1^. ittniih^ ITTi) + inT) + liTT)in) + liT) , ITii) + liti) + liiT) 
1^) = "olTTDlii) + oil 7= j= h 02 7= ITT) 



V2 



V3 



(12.12) 



We now compute the Schmidt decomposition for the entangled state ()12.12p in the 
partition 1|2345: 



^IT) 



«o ^ aiim±mm±m + - ^n) 



T 



V2 



T 



[2...5] 

T 



+V^li) 



«i .III) + +^1^(111) + |||))in) 



(12.13) 



[2 51 [2 51 

where " and "' are the Schmidt vectors of the subsystem corresponding to qubits 2 



to 5, and A''^, A''^ are the Schmidt coefficients, given by N-^ := Oq + + := ^ + 



According to the protocol developed in |Vid03| . we have, following its notation, (see Eq. 
(fTTI]) above) 



[rWT] = (i,o); [rW^] = (o,i) 



(12.14) 
(12.15) 



We proceed now iteratively, computing the Schmidt decomposition for state ()12.12p in the 



^Along this section we will be using mainly the notation jt) instead of |0) and ||) instead of |1) for the 
sake of clarity, although sometimes and in the following section we will switch from one notation to the 
other. 
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partition 12|345. It easily follows 



I*) = v^ITT) 



TT 



V2 



+ 



[3...5] 
TT 



Ti+iT- 



-IT D + liT) 

V2 



«i .|^^in) + liT) 



02 



Ti+TT 



:|iTT) 



[3...5] 

TT+TT 



+v^|ii) ^g=|TTT) 



(12.16) 



Ti 



[3...5] 
TT 



[3 51 [3 51 [3 51 

where , and are the Schmidt vectors of the subsystem corresponding 



to qubits 3 to 5, and A^f^, 2N^^i^ and A''^ are the Schmidt coefficients, given by N^-^ :- 
«0 + 2iV|T+TT - V + ^TT ■= f • 



According to [Vid03| . we have in this case 



[Ag] = iy^, V^N^i+ih yWi) (12-17) 



(12.18) 
(12.19) 




We proceed now with the partition 123|45 of the state |^'). Eq. (jl2.12p is indeed 
the Schmidt decomposition for this partition, where the Schmidt vectors of the subsystem 
associated to qubits 4 and 5 are <^f^f := ^f^f^^ := liMU, and := |tT). 

With the notation (...) we denote here the combination maximally symmetrized. 

Following |Vid03j . we have 
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[ALI] = (ao,ai,a2) 



L 0:20:3] 



/N. 



TT 





1 



\ 





Ti+iT 







(12.20) 
(12.21) 








\ 







(12.22) 



Finally we compute the Schmidt decomposition of |\1') in the remaining partition 1234|5, 
which turns to be 



«olTTT)li) , aiKTU))! T) 



+ 



li> 



+ ViV4:l,T 



«i|(m))|i) , a2|(TU)>lT> 



+ 



4:1,T 



IT), 



(12.23) 

(12.24) 



where the Schmidt vectors associated to the subsystem of qubit 5 are -=11) 

^4^1 1 •= IT), aiid the Schmidt coefficients are N^-i^^ := + iV4:i^-f := ^ + Q!|. 
In this case we have 



[aL1] = (V^,V^) 

/ \ 

1 

1 



as 04 J 







V 

/ 1 



L 03 04 J 



V 



Finally we have the remaining matrices [Faj ] and [F, 



V^4:l,i 




[5]Ti 



\ 

1 



V2JV4:1,T 





mi 







(12.25) 
(12.26) 



(12.27) 



(12.28) 
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The last step left is to match the X's and T's with the V's (isometries) of Eq. (jl2.ip . 
An easy comparison between Eqs. (jl2.ip and (jl2.4p yields 



0[1]-I«J 



0[2]Jy 



0[3]J«J 
















1 

V2 








V2 j 









\ 









3 





1 

V2 / 

































V 








1 




^ 



0[2]J«J 



7V- 



TX+J.T 
Aft 








0[3]J»J 














02 











n+J-T 



0[4]Jy 



0[5]Jy 




0[4]J«J 



0[5]J«J 




(12.29) 



The interested reader may straightforwardly check the isometry condition + 
-^3x3, A; = 1, 5 that the matrices verify, as expected. Some remarks are to be 



pointed out: 



The dimension of the ancilla Hilbert space is in this case Z?o = 3 (the matrices V 
explicitly calculated have dimension 3x3), in concordance with the results of the 
previous section: the dimension equals the number of clones M for the universal 
symmetric cloning (3, in this particular example). The global dimension for this case, 
D, is, according to our protocol stated in previous section, Z? = 2 x Dq = 6. 

Notice the change of order between the matrices V and Vidal's matrices F and A: we 
have expressed the V matrices in order to be applied from right to left, sequentially, 
upon the initial state := (1,0,0)"^. The matrices 1 to 3 create the three clones in 
a highly entangled state with the ancilla, and the remaining matrices 4 and 5 uncouple 
the clones from the ancilla by creating the two anticlones. The final state of the ancilla 
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after successive application of the five matrices ^0(1] results \ipF) '■= (1,0,0), as the 
interested reader may readily check. 

• The previous matrices are the corresponding to the sequential generation of state 
(|12.12p . which is the cloning of state ||). In order to clone an arbitrary superposition 
a|t) + , according to our protocol, the sequential generation of the cloning of state 
II) is needed. This state is obtained from state (jl2.12p by doing the substitution 1. 
Accordingly, the matrices V^i^j are obtained from the ones (|12.29|) by V^i^j = ^o[A;]' 
where i := i © 1. The global unitary operations are 6x6 matrices built from Vgjfc] and 
^i[fc] by diagonal blocks, according to our protocol [TJ 



12.3 Isometries in the general case 1 ^ M 

For the symmetric universal cloning 1 ^ M we now give in detail the 2D— dimensional 
matrices V^^ driving us to a concrete sequential scheme. The procedure for the calculation 
in this case follows the same lines of the previous example, although here the method relies 
in a highly symbolic notation based on maximally symmetrized states. 

To explicitly obtain the isometries, we firstly compute the Schmidt decomposition (see 
Appendix for the multipartite state |(M — j) | j l)s, where here \k ] k' l)s denotes the 
maximally-symmetrized state with k qubits | and k' qubits |. We will be interested in the 
Schmidt decomposition for the 1|2...M bipartition, that gives 



m-j) T j i)s = \i^^mM-j - 1) T J i)s + Jj^mM-j) t u - 1) ds. 



(12.30) 



The coefficients y ^^^^ and y -j^ follow directly from straightforward combinatorics. 

Expression ()12.30p will be the basis for our computation of the isometries. 

We consider now the general state with M clones and M — 1 anticlones for the initial 
qubit t |GM97j 

M-l 

IT) ^ IV') = 2Z -i) T,J i)s ®\{M-j-i) i,j T)5, (12.31) 
i=o 

where the first vector in the explicitly indicated tensor product corresponds to the clones 
and the second one to the anticlones, and Qj — • /^iM-^L 



J ~ y M(M+i) • 
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The Schmidt decomposition of (|12.3ip for the V\2...2M — 1 bipartition, is, according to 
(fT230|) 



H) = IT) 



ao\{M - 1) ])s\{M - 1) i)s + ai\l^^\{M - 2) T 1 i>5|(M - 2) j 1 T)5+ 



.. + aM-i\lj^\{M - 1) i)s\{M - 1) T)5 



a,^l — \{M-l) T)5|(M-2)ilT)s 



+a2^^m - 2) T 1 Dsm - 3) i 2 T)5 + ... + aM^i^J T (M - 2) 1)5 
x|(M-l) T)5 



(12.32) 



The (normahzed) Schmidt vectors of the previous decomposition associated to the subspace 
of qubits 2...2M — 1 are directly read from it and are given by 



'2...2M~1] 



ao\{M - 1) ])s\{M - 1) l)s + aixl^^m - 2) j 1 i)s 



x|(M - 2) i 1 ])s + ... + aM-i\lj^\{M - 1) i)s|(M - 1) ])s 



(12.33) 



'2...2M-1] 
I 



- 1) T)5|(M - 2) i 1 T)s + «2a/^|(M - 2) T 1 i)s 



x|(M - 3) i 2 T)5 + ... + aM-i\l^^\l T (M - 2) l)s\{M - 1) T>5 



, (12.34) 



where A^^ ■= «o + " + «2(^ - 2)/M + and iV| := af/M + 2a^/M + .... 

According to the protocol developed in |Vid03j . we have, following its notation, (see Eq. 
(fTTIj) above) 



[AS] = (^/iW>^/^) 

[rS^] = (i,o); [rWi] = (o,i). 



(12.35) 
(12.36) 



We now proceed iteratively, just as in the 1^3 case shown above. We compute 
the Schmidt decomposition for every bipartition, 12. ..k\k + 1...2M — 1, k = 2,...,2M — 2, 
extract the Schmidt coefficients, \^''\ and the Schmidt vectors, <i)[fc+i...2A/-i]^ g^^^^ applying 
Vidal's recipe, we obtain the rl*']*'' matrices. We made the calculation along this lines. 
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M 
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by splitting it in 3 steps: the partitions with k = 1,...,M — 1, that create the clones. 
The partition with k = M, that attains the multiqubit state of the clones with highest 
entanglement with the ancilla of all the sequential protocol, and the remaining partitions 
with fc = M + 1, 2M — 2, that create the anticlones. The need for this splitting of the 
computation is due to the different role of clones and anticlones: they are not mutually 
symmetric. The result of our computation is 



1 < n < M - 1 



N, 



(n-l)T 



N, 



N, 



(n-l)Tli 



(n-2)TU 



■p[n]J,^[n] ^ 



n = M 








N, 



("-i)TU 



(n-l)T 





N, 



N, 



(?i-2)T2i 



(n-2)TU 







y/N^M-l)} 









N, 



(M-2)TU 



\ 




lT(n-l)i r, 
^(n-l)i ^ / 

/ nxn 



+1 







N, 



nxn+1 








^J^{M-k-V)^ki 



CHM-\ 



h/M_ 



/ MxM 
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: 



M-fe)T(fc-l)i 



V 

< i < M - 2 



/ 



p[M+l+i]T_i^[M+l+i] 



1 



M-l-j 







■v/^lT(M-2)i 





\ 




MxM 



M-l~j 



V 

/ 1 





p[Af+l+iU^[Af+l+i] 
















1 / 
\ 



M~jxM-j-l 



M-j-k 
M-l-j 









V 







1 



M-l-j 





M-jxM-j-1 

where the coefficients -/V/ji^/j are the normalization coefficients of the corresponding Schmidt 
vectors as outhned for the 1^3 case in this chapter, and previously in this section. They 
are given by 



1 



M-l 

(M^ 

k=j ^i+j' 



(12.37) 



The corresponding isometrics Vj^,] for the universal symmetric cloning 1 ^ M of qubits, 
built upon the T matrices and A vectors, according to our protocol, are given in Table [T2TTJ 
The matrices Vj^* coincide also with the ones for the symmetric phase-covariant cloning 



just by doing the substitutions M 



^^2^ and CKj 



7j := 



' fc+l \fk\ 
Kk+l-i)\j) 
(2k + l\ 
\ fc+l / 
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A; = 


k = l 


roiiiJ,,= 


1 


5ijC{2-i,i-l) l<i,j<2 
^ Sij otherwise 


( 5i,3_jC(2-i,i-l) l<i,j<2 
1 -ij 5ij otherwise 


r , ~\ 


\ 


C r C(n + l-i,i-l) i ^ • • ^ 
Sij c(n-i,i-l) ^<hJ<n 

Sij otherwise 


1 


75 i = l,,- = n + l 

Si,j+i cfn-jj-l) 2<i<n + l;l<i<n 






ij . ^ l<i,j <M 

C(M-i,i-l)^(,^l) 


r <5,,,+i ^ — == i<i,j<M 




1 


l<i<M-n;2<j<M-n + l 
i = M-n + l;l<j<M 
ij otherwise 


1 


■ ^^^V'^M+n"' l<i,J<M-n 

i = M- n + l;l<j<M 
-^Sij otherwise 



Table 12.1: Explicit form of the isometries for universal symmetric quantum cloning 1 — > M. 



Here C{i,j) := = J^^^ri |a^|2L.^^ (^e use the convention Q = 

if g > p) and 1 < n < M — 1. Complementarily, we also have V^^^^^j = Vj^^p where by i we 
indicate i := The reader can readily check that the minimal dimension for the ancilla 

is 2 X M. When sequentially applying these matrices to the initial state \(pi) of the ancilla, 
one can check, as expected, that if we were to stop at the Mth step, the M clones would 
have already been produced with the desired properties, although in a highly entangled 
state with the ancilla. To arrive at a final uncoupled state, the remaining M — 1 anticlones 
must be operated upon by the ancilla. 

In conclusion, following the idea that some (certainly not all) unitary operations upon 
a multiqubit state can be implemented sequentially by successive interactions between each 
qubit and an ancilla, we have shown how to reproduce both the symmetric universal and 
symmetric phase-covariant cloning operations. For the universal cloning we have proved 
that the minimal dimension for the ancilla must be -D = 2M, where M denotes the number 
of clones, thus showing a linear dependence. The original state must be encoded in a 
i?— dimensional state. For the phase-covariant case, the dimension D of the ancilla can be 
reduced to D = M+1, even lower. In both cases, the ancilla ends up uncoupled to the qubits. 
Along similar lines, this sequential cloning protocol can be adapted to other proposals, such 
as asymmetric universal quantum cloning or other state-dependent cloning. This procedure 
can have notable experimental interest, since it provides a systematic method to furnish 
any multiqubit state using only sequential two-body (qubit-ancilla) operations. 



Chapter 13 



Inductive classification of 
multipartite entanglement under 
stochastic local operations and 
classical communication 

A comprehensive understanding of entanglement is still lacking, mainly because it is a 
highly counterintuitive feature of quantum systems (non-separability |Bel87j ) and because 
its analysis can be undertaken under different, although complementary, standpoints. As 
prominent examples the subjects of deciding in full generality whether a given state carries 
entanglement or not and how much entanglement the system should be attributed to are 
vivid open questions (cf. e.g. jBru02j and references therein). This state of affairs is critical 
in multipartite systems, where most applications find their desired utility. 

Among others, part of the efforts are being dedicated to classify under diversely moti- 
vated criteria the types of entanglement which a multipartite system can show. It is in this 
sense desirable, independently of these criteria, to have classification methods valid for any 
number N of entangled systems. One of these most celebrated criteria to carry out such 
a classification was provided in [DVCOOj . In physical terms Diir et al. defined an entan- 
glement class as the set of pure states which can be interrelated through stochastic local 
operations and classical communications (SLOCC hereafter) or equivalently, as those pure 
states which can carry out the same quantum-informational tasks with non-null possibly 
different probabilities. They also proved the mathematical counterpart of this characteri- 
zation: two states ^ and ^ of a given system belong to the same entanglement class if, and 
only if, there exist invertible local operators (ILO's hereafter; that is, nonsingular matrices), 
which we agree on denoting as such that ^ = F^^^ • • • F^^\'^). Moreover, they 
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provided the first classification under this criterion of tripartite multiqubit entanglement, 
giving birth to the two well-known genuine entanglement triqubit classes named as GHZ 
and W classes. Later on, exploiting some accidental facts in group theory, Verstraete et al. 
|VDMVd2] gave rise to the classification of 4-qubit states. 

Regretfully none of the previous works allowed one to succeed in obtaining a generaliz- 
able method. In the second case, the exploitation of a singular fact such as the isomorphism 
SU{2) SU{2) ~ 50(4) is clearly useless in a general setting; in the first case, the use of 
quantitative entanglement measures specifically designed for three qubits, as the 3— tangle 
[CKWOd] , to discern among different entanglement classes discourages one to follow up the 
same trend, since we would have to be able to build more generic entanglement measures, 
per se a formidable task. However, Verstraete et al. |VDM03| succeeded in this approach 
by introducing the so-called normal forms, namely those pure states such that all reduced 
local operators are proportional to the identity matrix. These authors also provided a sys- 
tematic, mostly numerical, constructive procedure to find the ILO's bringing an arbitrary 
pure state to a normal form. Furthermore, the use of these normal forms allowed them to 
introduce entanglement measures (entanglement monotones [VidOOj . indeed), which offered 
the possibility to quantify the amount of entanglement in the original state. For complete- 
ness' sake let us recall that classification under SLOCC is coarser than that using only local 
unitaries, that is in which every F^^'^ is unitary. Nevertheless relevant results in this realm 
can be found in the literature |GRB98l lAAC+OOi ICHSOOllGXFWOGj . 

Here we offer |LLSS06] an alternative and complementary approach to the classification 
under SLOCC based on an analysis of the singular value decomposition (SVD) of the co- 
efficient matrix of the pure state in an arbitrary product basis. The coefficient matrix is 
chosen according to the partition 1|2 . . . with the subsequent goal of establishing a recur- 
sive procedure allowing one to elucidating the entanglement classes under SLOCC provided 
such a classification is known with one less qubit. The key feature in this scheme is the 
structure of the right singular subspace, i.e. of the subspace generated by the right singular 
vectors of the coefficient matrix, set up according to the entanglement classes which its 
generators belong to. As a secondary long-term goal, the approach seeks possible connec- 
tions to the matrix product state (MPS) formalism (c/. jEckOSj [PCVWCOGj and multiple 
references therein), which is becoming increasingly ubiquitous in different fields such as 
spin chains |AKLT87j , classical simulations of quantum entangled systems |Vid03) , density- 
matrix renormalization group techniques |VPC04] and sequential generation of entangled 
multiqubit states [SSV+05l[SHW+06| . 

This chapter uses a more mathematical notation. The canonical orthonormal basis in 
will be denoted by {ej}j=i^,..^7v (correspondingly in physics the kets |j — 1)). Normalization 
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is not relevant in elucidating the entanglement class which a state belongs to. Thus we will 
deal with unnormalized vectors and non-unit-determinant ILO's. In the SVD of an arbitrary 
matrix (cf. appendix), V and W will denote the left and right unitary matrix, whereas T, 
will stand for the diagonal possibly rectangular matrix with the singular values as entries. 
In the multiqubit cases, we will agree on denoting by small Greek letters (j),ip,... vectors 
belonging to C^, whereas capital Greek letters "f, <I>, ... will denote a generic entangled state 



The chapter is organized as follows. In section 113.11 the entanglement of two qubits is 
revisited with a reformulation of the Schmidt decomposition criterion in terms of the singular 
subspaces. In section [13.21 the extension to the three-qubit case is developed in detail and 
the principles of the generalization to multipartite and arbitrary-dimension systems are 
discussed in section flS. 31 Some concluding remarks close the chapter. 

13.1 Bipartite entanglement 

13.1.1 The Schmidt decomposition criterion revisited 

The determination of entanglement of pure states of bipartite systems in any dimensions, 
in general, and in two dimensions (qubits), in particular, was solved long ago with the aid 
of the well-known Schmidt decomposition |Sch06l |EK95| (see Appendix [A]), by which any 
bipartite state can be written as a biorthogonal combination 



where Ai > A2 > • • • > for all n and Ni denotes the dimension of subsystem i. If A„ = 
except for only one index Ai 7^ 0, then the state is a product state; on the contrary, if An 7^ 
for two or more indices, then the state is an entangled state. Furthermore, A„ coincides with 
the common eigenvalues of both reduced density operators. Thus, to practically determine 
the entangled or separable character of a given pure state all we must do is to compute the 
spectrum of pi or p2 or equivalently to analyze the dimensionality of their ranges. This is 
the backbone in the study of 3— partite entanglement carried out in [DVCOOj . 

Followingly in order to pave the way for a generalization to multipartite systems, we 
will reformulate the Schmidt decomposition criterion for bipartite systems focusing upon 
the subspace generated by the singular vectors. We need the next 



in C2. 



min(Afi,Ar2) 





n=l 



Inductive classification of multipartite entanglement under stochastic local operations 
142 and classical communication 



Definition 5. We will denote by 23 (resp. SUJ the subspace generated by the left (resp. 
right) singular vectors, i.e. 23 = spanjfi, . . . , v^} (resp. W = span{u)i, . . . , w^} )■ 
We can now state the fohowing 

Theorem 1. Let ^ G C™ C" and C(^') denote the matrix of coefficients of ^ in an 
arbitrary common product basis. Then ^ is a product state if and only if dim W = 1 (or 
alternatively dim 23 = Ij. 

Proof. Let {ej}j=i^...^m and {/j}j=i,...,n denote bases in C"^ and C", respectively. Then any 
vector ^ G C™ (8) C" can be written as 

m n 

= Y^Y,^i:i('i^fj' (13.2) 

i=i j=i 

where Qj are the complex coeficients of ^, which we arrange as: 



^ Cll ■■■ Cin^ 



(13.3) 



yCjni . . . Cmn J 

The matrix C(^') = C always admits a SVD, given by C = ySVl^^, where V and W are 
unitary matrices and S is a diagonal matrix with entries ak (the singular values, indeed). 
Thus 

min(m,n) 

Cij = ^ VikakW*^. (13.4) 
k=l 

Inserting (jl3.4p into (|13.2p and identifying new bases {ej}i=i^2 and {fj}j^i 2 we arrive 
at the well known Schmidt decomposition 

min{m,7i) 

^= Yl '^kek^fk- (13.5) 
fc=i 

The number of non-null singular values coincides with the rank of S, which in turn 
coincides with the dimensions of 23 and 2IJ. 

□ 
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From the proof we can deduce a practical method to recognize where a bipartite system 
is entangled or not: 

Corollary 1. Let ^ € C"* (X" C" denote the state of a bipartite quantum system and C{^) 
its coefficient matrix in an arbitrary product basis. Then is a product state if, and only 
if i{C{^)) = 1. 

13.1.2 Classification of two-qubit entanglement under SLOCC 

We only need one further tool to find the classification of bipartite entanglement under 
SLOCC, which is established as follows: 

Proposition 4. Let ^f,^ S (^C^ denote two two-qubit states related by SLOCC, i.e. 

^ = (13.6) 

where F^'^'^ and -Pt^' are non-singular operators upon C^. Then their corresponding coefficient 
matrices C, C in an arbitrary product basis are related through 

C = (FW^y)S(F[2]tvF)t. (13.7) 

Proof. Just substitute "if = Ylij=i fj ™ (|13.6p and identify indices. □ 

The key idea in our analysis is to recognize the effect of the ILO's upon the singular 
vectors. If Vj (resp. Wj) is a left (resp. right) singular vector for the matrix coefficient 
C, then F^^^'^{vj) (resp. F^'^^'^ (wj)) is a left (resp. right) "singular vector" for the new 
matrix coefficient C. In order to ease the notation, we will agree hereafter on relating ^ 
and ^ through ^ = F^^]^ F[^]"f(^'), which allows us to drop the transpose and Hermitian 
conjugatioio in future considerations. 

The case of two qubits is elementary, since there is no much space to discuss. The bases 
in which the coefficient matrix will be expressed are the canonical orthonormal basis {ei, 62} 

^Notice that they cannot rigorously be singular vectors, since the ILO's are not in general unitary, thus 
they do not preserve the orthogonality of {vj} and {wj}. We will understand these "singular vectors" in a 
loose sense, in which they substitute the original singular vectors in the SVD of the coefficient matrix. 

^The transpose and Hermitian conjugation are referred to the chosen product basis in which C is 
constructed. 
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in C^. Only two options are present: either dim 221 = 1 or dim2IJ = 2. In the first case, 
after choosing Ft^l such that 

= — ei, (13.8a) 

F^'^\wi) = ei, (13.8b) 

the new coefficient matrix will turn into C = ( q o ) ' which corresponds to the product state 
^ = ei (g) ei . We will agree on stating that belongs to the entanglement class denoted by 
00. 

In the second case, where cti > (72 > 0, after choosing F^^^ and F^l such that 

FW(^vi) = —ei, F^^\v2) = —e2, (13.9a) 

C7i (72 

fP1(u;i) = ei, F^'^\w2) = 62, (13.9b) 

the new coefficient matrix will be C = ( q ? ) , which corresponds to the entangled state 
^ = ei (g) ei + 62 (8) 62. Now we say that * belongs to the class 

The reader can readily check by simple inspection how in the first case the canonical 
matrix C has rank one, whereas in the second it has rank 2, as expected. In summary, only 
two classes are possible, namely 00 and 



13.2 Tripartite entanglement 

The classification of tripartite pure states is performed along the same lines, namely choosing 
the ILO's so that the final coefficient matrix reduces to a canonical one. In order to 
find such canonical matrices, we must be exhaustive in the considerations of all possibilities 
when discussing about 03 and W. 

The analysis of tripartite entanglement can be undertaken upon three possible coefficient 
matrices, arising from the three different ways to group the indices, that is, since = 
Ylii 12 13=1 2^n»2i3^n ® %5 where as before {e^} denotes the canonical orthonormal 

basis in C^, we have 

CW^Ci|23 = '^'^ '^''V (13-lOa) 

\C211 C212 C221 C222/ 



c(2)^c2|i3 = , (13.10b) 

\Cl21 Ci22 C221 C222/ 



C(3)^C73|i2 = "''M. (13.10c) 

\Cll2 Ci22 C212 6222/ 
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There is no loss of generality in choosing one of them, since the analysis will be exhaus- 
tive. Hereafter we will choose C = C^^h Notice that now the left singular vectors of C 
belong to whereas the right singular vectors are in C^^C^. Also, we immediately realize 
that only two possibles options arise, namely dimSH = 1 or dim SB = 2, since there are at 
most two positive singular values. The recursivity appears when classifying the different 
structures which the subspace 2IJ can show. The classification of these subspaces is per- 
formed according to the entanglement classes which their generators belong to. In order to 
do that we need the following result, which was firstly proved in the context of entanglement 
theory in |STV98j . We offer an alternative proof in order to illustrate our methods. 

Proposition 5. Any two-dimensional subspace in contains at least one product 
vector. 

Proof. Let y be a two-dimensional subspace of C^. With no loss of generality two 
entangled vectors can be chosen as generators of V with coefficient matrices given by Ci = I 
and C2 being an arbitrary rank-2 matrix in the product canonical basis. Then it is always 
possible to find non-null complex numbers a and (3 such that al + (3C2 has rank onqf|. □ 



In other words, this proposition shows that span{\I'i, ^'2} always equals either span{(/)i(g) 
'4'i,(p2 'S' 4^2} or span{(/) f^ip,'^}, where implicit are the assumptions that different indices 
denote linear independence and in the last case only one product unit vector can be found. 
Thus, with the same convention, the right singular subspace SU can show six different 
structures, namely span{(/) (?) ip}, spanl^f}, span{(/) •tpijCp ® 1P2}, span{0i ip,4>2 ip}, 
span{(/)i 'SiijJi,(p2 (8) V'2} and span{(/)(8)V) We can now state our result, already contained 
in [DVCOO| with different criteria: 

Theorem 2. Let G O O be the pure state of a tripartite system. Then ^ can 
be reduced through SLOCC to one of the following six states, which corresponds to the six 
possible entanglement classes, according to the following table: 



Proof. We discuss depending on 2IT: 



^Notice that —a/f) must be chosen to be an eigenvalue of C2. 
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Clciss 


Canonical vector 


Canonical matrix 


an 


000 


ei ei (g) ei 




/'I o^ 
^0000^ 




span{^ V'} 




ei (8) ei ei + ei 62 62 




n i\ 
^0000, 




spanjvl/} 


02*^3 


ei (g) ei ei + 62 ei (g) 62 




?i o^ 
^0100, 




(/)0C2 


03*^2 


ei ei ei + 62 62 ei 




(I o^ 
,0010^ 




€^0 V 


GHZ 


61 ei ei + 62 62 62 




A o\ 
^0001, 




span{0i 02 V'2} 


W 


61 61 62 + 61 62 61 + 62 61 61 




?o 1 1 o^ 
^1000^ 




span{0i V'l, ^"1 



TaftZe Genuine entanglement classes for three qubits 

1. 2U = span{^(g)^}. In this case, wi = (j)®-^;. Choose the ILO's F^, A; = 1, 2, 3 so that 

= — ei, (13.11a) 

0"! 

f[2](</,) = ei, (13.11b) 
f[31(^) = ei. (13.11c) 



Then the new coefficient matrix will be 



/l o\ 









/ 


c = 










1° 


•J 


V 























(71 





(13.12) 



7 



(13.13) 



which corresponds to the state ei(8)ei(8)ei, and where the dots • indicates the irrelevant 
character of that entry. 
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2. W = span{^}. In this case tt^i = 0i ® + (/'2 ^S) ip2- Choose the ILO's so that 



FW(^i) = -ei, 
F[3l(V;i) = ei, 



F[21((/,2) =62, 
i^'^k^2) =62. 



(13.14a) 
(13.14b) 
(13.14c) 



Then the new coefHcient matrix will be 













( 











/l l\ 



7 



which corresponds to the state ei ei (8) ei + ei (8) 62 (8) 62. 



(13.15) 



(13.16) 



3. 2U = ^ (8 = span{0 tpijcf) V'2}. In this case wi = fj-ucf) ipi + lii2<P <8 V'2 and 
W2 = <8 V'l + M22^ <8 V'2) where the matrix [/Xjj] has rank 2, since wi and W2 are 
linear independent (orthonormal, indeed). Choose the ILO's so that 



Ff^(t;i) = ^ei, Ffi(^2) = ^e2, (13.17a) 

= [i.W(e,) i.W(e2)] = (13.17b) 

F[i1=fWfW, (13.17c) 

F[2](<^) = ei, (13.17d) 

fPI(V'i) = ei, F[^1(V'2) = 62. (13.17e) 
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Then the new coefficient matrix wiU be 



C 



1 * 

Mil 


/^12 1 


-1 


0] 


/ 


YM21 






i 


V 



CTi 
(72 



/x*2 0\ 

11*21 /^22 



V • 



• 7 



'^1000^ 



(13.18) 



1 0^ 

which corresponds to the state ei (8> ei (8) ei + 62 (8) ei (8> 62- 

4. 2U = (g) V' = span{0i V'}- In this case w\ = //n^i (8) + A*i2^2 ^ V' and 

'W2 = /U2101 iX" V' + /^'2202 "S" V') where the matrix [fiij] has rank 2, since wi and W2 are 
linear independent (orthonormal, indeed). Choose the ILO's so that 

FW = [FW(ei)FW(e2)] = K-r\ 

F[2l(0i) = ei, FP](</>2) = e2, 
F[3](^) = ei. 



(13.19a) 
(13.19b) 
(13.19c) 
(13.19d) 
(13.19e) 



Then the new coefficient matrix will be 



C 



II* II* 
Mil M12 

yM21 M22 y 



'^1 0^ 




Mil 
M21 M22 



V • 



7 



(13.20) 



yO 1 0^ 

which corresponds to the state ei (8 ei (8 ei + e2 <8) 62 (8 ei . 
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5. 2H = span{(/)i (g) ^i,(/)2 <^ ^2}- In this case wi = /iii0i V"! + Mi202 ® V'2 and 
■"^2 = A'21^1 ® V'l + A'2202 ® V'2) where the matrix [/x^j] has rank 2, since wi and ■u;2 are 
linear independent (orthonormal, indeed). Choose the ILO's so that 



F|'k^i) = ^ei, F|^J(^2) = ^e2, (13.21a) 

= [i.W(e,) Fl'\e,)] = [^^*^]-\ (13.21b) 

FW=fWfW, (13.21c) 

F[2l(0i) = ei, FP](,/>2) = e2, (13.21d) 

F[3](^i) = ei, F[31(V2) = 62. (13.21e) 



Then the new coefficient matrix will be 



C 





^fil^ /it2^ 
nil 1J.I2 



'^1 0^ 



1 



(13.22) 



which corresponds to the state ei (8) ei ei + 62 (8) 62 ® 62. 



6. 2B = span{0i (g V'l)^}- In this notation, implicit is the assumption that only one 
product unit vector can be found in 2IJ. In this case ^ can be chosen so that 
= 01 <8) V'2 + 02 <8) V'l ■ Thus the singular vectors can always be expressed as wi = 
Mil (01 <8) V'2 + 02 <8) V'l) + M1201 <8) ■01 and W2 = At21 (01 ^02 + 02 <8) V'l) + At2201 <8) ■01) 
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where the matrix [fJ-ij] has rank 2, since wi and 11)2 are hnear independent (orthonor- 
mal, indeed). Choose the ILO's so that 

i^fk«i) = ^ei, fW(^2) = J-e2, (13.23a) 

fW = [FW(ei) FW(e2)] = [^^*,]-\ (13.23b) 

Fm=FPFP, (13.23c) 

FPI {cj)i) = ei (</.2) = 62, (13.23d) 

F[3] (^^) = jrl3] (^2) = 62. (13.23e) 



Then the new coefficient matrix wih be 



c = 






-1 


1" "]■( 






^22/ 




1° ^) \ 




^ fJ'U Mil A^n 0^ 

/^22 /^21 ^^21 

V • • • 



V 



0110 
1000^ 

which corresponds to the state ei CjD 61 (8) 62 + ei (8) 62 CjD ei + 62 (8) ei (8) 61. 



7 



(13.24) 



Since there is no more options for the subspace 2IJ we have aheady considered all possible 
alternatives. □ 

In conclusion, we have found that there are six classes of entanglement, named after 
|DVCOO| as 000, Oj^^^j^, GHZ and W . The theorem also indicates how to practically clas- 
sify a given state ^: compute the SVD of its coefficient matrix and elucidate the structure 
of s^an{wi,W2}- We include a further proposition comprising the practical implementation 
of this result. We need to introduce the following definition. 

Definition 6. Let Wj = 61 (8 Wji + 62 Wj2 S (8 he an arbitrary vector. We associate 



a two-dimensional matrix Wj to wj by defining 



Wj = [Wji Wj2]. 



(13.25) 
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This definition will be mainly applied to the right singular vectors of the coefficient 
matrix C. As usual, the singular values of C will be denoted by a^, in nonincreasing order, 
and (t{A) denotes the spectrum of a matrix A. Our proposal to implement the preceding 
result is 

Theorem 3. Let ^' denote the pure state of a tripartite system and (7^*^ its coefficient 
matrix according to the partitions i\jk (cf. (jl3.10ap - (|13.10cp ). Then 

1. IfT{C^'^) = 1 for all i = 1,2,3, belongs to the 000 class. 

2. IfT{C^^^) = 1 and t{C^^'i) = 2 for k = 2,2,, ^ belongs to the Oi^-^g class. 

3. IfT{C^^'^) = 1 and r(C(^)) = 2 for k = ^ belongs to the 02^-^3 class. 

4. //r(C(3)) = 1 and r{C^^^) = 2 fork = 1,2,^ belongs to the ^^^"[^ class. 

5. Ifv{C^'^) = 2 for alli = 1,2,3 andv{Wi) = i{W2) = 1 , ^' belongs to the GHZ class. 

6. IfT{C^^) = 2 for all i = 1,2,3, r(VFi) = 2, r{W2) = 1 and (t{W{^W2) is non- 
degenerate, ^ belongs to the GHZ class. 

7. Ifi{G^^) = 2 for all i = 1,2,3, r(VFi) = 2, t{W2) = 1 and a{W{^W2) is degenerate, 
^ belongs to the W class. 

8. Ifi{G^^) = 2 for all i = 1,2,3, v{Wi) = 2, r{W2) = 2 and (t{W{^W2) is non- 
degenerate, ^ belongs to the GHZ class. 

9. IfT{G^^) = 2 for all i = 1,2,3, t:{Wi) = 2, v(W2) = 2 and a(W^^W2) is degenerate, 
^ belongs to the W class. 

Proof. The first four cases are elementary, since it is a matter of detection of the vector 
which factorizes. The final five cases correspond to true tripartite entangled states. If 
^(Wk) = 1 for k = 1,2, it is clear that there exist two product vectors belonging to 2IJ, thus 
^' belongs to the GHZ class. If r(PFi) = 2 and t(W2) = 1 we need to check whether an 
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ILO applied upon the first qubit can reduce the rank of the transformed Wi. As it can be 
deduced from the preceding proofs, an ILO upon the first qubit amounts to constructing a 
hnear combination between the two right singular vectors, whicli is equivalent to find new 
matrices Wj = F^^Wi + F2^^W2, with j = 1,2. If r{Wi) = 2, then by multiplying this 
expression to the left by W^^ , we have 

F[]h + W{^W2. (13.26) 

It is immediate to realize that it is possible to reduce the rank of Wi to 1 and to choose 
F,7 such that fI^] is nonsingular provided the spectrum of is non-degenerate, in 

which case ^ belongs to the GHZ class. If the spectrum is degenerate, thus both eigenvalues 
being null, no further reduction is possible and ^ belongs to the W class. 
Finally if r(I^x) = ^(^^2) = 2, reasoning along similar lines if both eigenvalues of W^'^W2 
are equal, only one rank can be reduced keeping the nonsingularity of F^^l and ^' belongs 
again to the W class, whereas if the eigenvalues are different, both ranks can be reduced to 
1 keeping the nonsingularity of F^^l and ^ belongs to the GHZ class. □ 

As a final remark let us indicate how close, despite the apparent differences in the 
approach, our analysis runs to that performed in [DVCOO] : the ranges of the reduced density 
operators are indeed generated by the corresponding singular vectors, and the study of these 
ranges drove them and has driven us to the same final result. The change of method is 
motivated by the attempt to find a generalizable criterion not using entanglement measures 
specifically built upon the number of qubits of the system, such as the 3— tangle [CKWOO] . 
With this approach it is not necessary to consider at any stage the reduced density matrices 
and entanglement measures upon them. 

13.3 Generalizations {N > 4) 

The generalization of the preceding approach to pure states of arbitrary multipartite systems 
is two-folded. On the one hand, the generalization to multiqubit states can be implemented 
inductively: 
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Theorem 4. // the entanglement classes under SLOCC are known for N qubits, the cor- 
responding entanglement classes for A'" + 1 qubits are also known. 

Proof. We proceed by induction. We have proved in preceding sections that this statement 
is true for N = 2 and have explicitly found the entanglement classes for = 3. For a 
given (N + l)-qubit system, write the coefficient matrix Ci|2...Ar-|-i = C. Because of the 
induction hypothesis one knows in advance the classification of the right singular subspaces 
of C according to 2IT = spanj^j} if dimSU = 1 and 2)3 = spanj^j, ^j} if dim!2B = 2, where 
each and belong to one (possibly the same) of the entanglement classes of qubits. 
Choose the ILO's i^'^^ (8) • • • (8) ^[^+-"^1 so that the two first columns of W (the transformed 
right singular vectors) are expressed as linear combinations of the canonical vectors of the 
entanglement classes corresponding to the structure of 2U and choose the ILO F^^'^ so that 
VTiW'^ drops out as many non-null entries as possible (typically F^^l will be the inverse of 
a rank— 2 submatrix of W^). The result is the canonical matrix for an entanglement class 
of iV + 1 qubits. □ 

There is an important remark in the preceding inductive construction, already stated in 
[DVCOO| and explicitly shown in |VDMV02] : there will be a continuous range of states with 
a similar right singular subspace but with no ILO's connecting them. Let us illustrate this 
peculiar fact with an explicit example. When considering 4— partite entanglement, there will 
exist 45 a priori structures of the right singular subspace of the coefficient matrix, arising 
from 6 possible one-dimensional right singular subspaces 2n = span{^}, where ^ belongs 
to one of the six entanglement classes of = 3, times four possible sites for the fourth 
added qubit, plus 21 possible bidimensional right singular subspaces 2)3 = spanj^'i, ^'2}; 
corresponding to the (^^2^^) ways to choose the classes for A^ = 3 which ^'i and ^'2 belong 
to. An example will be 2)3 = spanjOOO, with the already convention that only one 

product vector and no Oj^'jfc belongs to 2)3, i.e. 2)3 = spanjc^i v^i ■01 j '?^2 <8 ¥'2 "SD V'2 + 1^2 ® 
(^2 <Xi V'2}) where the vectors with "are pairwise linearly independent. In order to only have 
one product vector and the rest being GHZ vectors, we must have (up to permutations) 
2)3 = span{(/) (d(p(d'4'',4"^^'^'4' + 4"^'f'^ "0}) with -0' / ip, -0. 
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Recalling that 



W2 



fin4> (8) ® V'' + /^12 ® ® "0 + ® </5 ® ■0) > 
fi2l4' ® ^ ® + 1^-22 {<t> ® ® i' + (f) ® ^ ® i^) ■, 



;i3.27a) 



(13.27b) 



where the matrix 

i?[2]^i7[3]^^[4] g^pjj ^^^^ 



, M21 M22 ) '^ill be non-singular, it is immediate to find ILO's 



i?[2] ^[3] _p[4] (-y^^-) ^ ^^^g^ (g) e2 V + Aii2 (ei O ei (g) ei + 62 <8) 62 O 62)(,13.28a) 

i?[2] (g ^[3] (g _p[4] (y^^) ^ e2 V' + /i22 (ei 61 O 61 + 62 62 (g) 62Xl3.28b) 

which corresponds to a coefficient matrix given by 

Mn^2 



C = VT. 



/^22 M2lV'l /^21^2 ^22 



v • • 



(13.29) 



8x8 



where the coefficients ipi corresponds to the coordinates of the transformed %[)' in the canon- 
ical basis. Choosing F^^l so that 

= (13.30) 



we arrive at 



C 



V'l ^^2 0^ 



^1 1^ 
which corresponds to the canonical vector 



(13.31) 



6l ® 6l 62 V'* + ^2 ® ei 



61 (g) 61 -I- 62 (gl 62 62 (8) 62 = (V"* / Cl, 62) 

= |001V*) + 11000) + 11111) (1^*) / |0),|1)) (13.32) 



Thus, different ip will yield different entanglement classes under non-singular local operators 
(g) • • • (gi F^^'^ . Notice that this vector belongs neither to the GHZ^ class nor to the W4 
class nor to the $4 class (containing the cluster state of four qubits -see below). It is a 
peculiar feature that two infinitesimally close states could belong to distinct entanglement 
classes, so a deeper elucidation of this point is on due and will be carried out also elsewhere. 
For the time being, we will agree on attributing all states reducible to ()13.32p by ILO's 
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FI^I eg) • • • (gii^'W, independently of the particular vector -0, the same entanglement properties 
under SLOCC and analogously for arbitrary A^-partite multiqubit systems. 

This allows us to find an upper bound for the number of genuine (A^ + l)-partite entan- 
glement classes. Firstly, notice that e.g. the right singular subspace SU = spanjOOO, 000} 
in the 4— partite case actually contains structures with different properties, namel}!^ 21? = 
(g) (and permutations), W = span{(/) (^i ipi, (/>(8> 992 V2} (and permutations) 
and 2IJ = span{<j)i (8) 931 (X" ^1 , 02 (X" 9^2 ® V'2} • All of them drives us to at least one factor qubit 
in the final canonical state, except one, that is, there will correspond one right singular 
subspace structure spanj^fi, ^'2} to each genuine (A^ + 1)— entanglement class. 

This is rigorously proved in the following 

Proposition 6. Let Wj\f be the right singular subspace of the coefficient matrix in an 
arbitrary product basis of an N -qubit pure state. If Wn is supported in a product space 
Wn = ijj ® Wm-i, then the state belongs to a product class 02^, where ^ denotes a class 
of (N — 1) -partite entanglement. 

Proof. Under the above assumption, wj = il)®Wj., j = 1, 2, with ■0 £ and Wj G C^^-^^^^ 
We can always find an ILO F'^' such that 

Wj eiCgiiij, (13.33) 

where also wj G C^(^~^\ hence Wn = En (g) Wn-i, where En denotes the Weyl matrix 
Ell = Since we can always write Sjv = En (8) T,i\i^i, the coefficient matrix can 

always be written as 

On = V^nW^ = V {Ell ® ^N-i) {Ell ® Wn-i)^ = En ® {vT^n^iWI^.^ ■ (13.34) 
The remaining ILO's Fill and FI^I j > 2, can always be chosen so that 

Cn = En^CN-i, (13.35) 

where C^-i denotes a canonical matrix of an (N — l)-partite entanglement class. This 

proves that the second qubit factorizes, as the reader may check. □ 
''As usual, different indices denote linear independence. 
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With appropiate permutations, this result apphes to any qubit. 

If we denote by M{N) the number of A^— partite entanglement classes, there will be at 
most 



genuine entanglement classes for + 1 qubits. Besides, the number of degenerate {N + 
1)— entanglement classes will be at most (A^ + 1) x M{N) (corresponding to the + 1 
possible factor positions which the (A^ + l)th qubit can occupy), thus 

Corollary 2. Let M{N) denote the number of N —partite entanglement classes under 



The equality will be in general unattainable, since, as in the case of tripartite entangle- 
ment, only a few distinct true entanglement classes exist, coming out from the only actually 
different structures which the right singular subspace can adopt (only two in the case of 
tripartite systems; cf. proposition [5|). 

Another benefit of the present approach arises when deciding whether two states belong 
to the same entanglement class or not. As an example, let us include a one- line proof 
that the 4-qubit GHZ state \GHZi) = (|0000) + |1111)) and the cluster state {cp^) = 
i (|0000) + |0011) + |1100) - |1111)) [BROlj do not belong to the same class [WZOOj . Their 
respective right singular subspaces are WGHZ4, = spanjei (8) ei (gi ei, 62 18) 62 (SD 62} and 
^(j>4 = spanjei (8> 62 ^~}, where ^'^ denote two-qubit Bell states. It is immediate to 
conclude that they are different, since none ej ej ej belong to OT^^ (write the coefficient 
matrix of a generic vector in SD^^ in terms of two coordinates a and /3 and check that it is 
impossible to choose the latter so that the matrix corresponds to ej (g ej (g ej). These states 
belong to the respective so-callecil GHZ4 and ^>4 classes, characterized by the above right 
singular subspaces. 

On the other hand, to find a wider generalization one can focus upon arbitrary dimen- 
sional entangled systems. The leit motiv is still the same, with the important exception that 
the dimension of the right singular subspace can grow up to the dimension of the Hilbert 
space of the first subsystem. Thus the analysis of the possible structures which 2IJ may 
adopt is now much more complex. 

We include as an illustrative immediate example the analysis of all entanglement classes 
under SLOCC of any bipartite (A'^i x A^2)-dimensional system: there exist min(A'^i, A^2) 
^The first one is named by a natural extension of the tripartite case; tlie second, after its representative 




(13.36) 



SLOCC. Then 



M{N + 1) < -M{N) [M{N) + 2N + 2,]. 



(13.37) 



\4>i)- 
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entanglement classes, which can be denoted as 00 = ^~[, ^2 , ^'3 , • • • , ^^jin^jVi N2)' "^^ose 
canonical states will elementarily be Yli=i ^i® ^i-, for each class They correspond to 
canonical matrices given by Yl\=i -^ni ^^^^ can state the following 

Theorem 5. Let ^ e C^i (g) C-^^ ^/jg 

pure state of a bipartite quantum system with 
coefficient matrix in an arbitrary product basis denoted by C{'^). Then belongs to the ^'^ 
class, k = 1,2, . . . , min(A^i, N2), if, and only if, r(C(*)) = k = dimQJ = dimSH. 

Proof Let 9J = span{(f)k}k=i,- ,n<min{Ni,N2) and 2IJ = span{ipk}k=i,- ,n<min{Ni,N2)- Choose 
and fI^I so that 

F^'KM = —ek, (13.38) 

F^^\ipk) = ek. (13.39) 

Then the coefficient matrix (in blocks) will turn out to be 

- I In 0jV2-n 
( / = 

^Ni-n OiVi-n.ATa-n j 

For more general cases, the difference stems solely in the higher computational complex- 
ity. 

We have developed a recursive inductive criterion to classify entanglement under SLOCC 
in multipartite systems in pure states which allows one to find the entanglement classes for 
N + 1 qubits provided this classification is known for TV qubits. The method rests on the 
analysis of the right singular subspace of their coefficient matrix, which is chosen according 
to the partition 1\2. . . N, hence a 2 x 2^~^ rectangular matrix. Then one must elucidate 
the classification of the one- and two-dimensional right singular subspaces according to the 
entanglement classes which their generators belong to. As a consequence, this construction 
reveals a systematic way to detect the entanglement class of a given state without resorting 
to quantitative measures of entanglement. In arbitrary-dimensional generalizations, the 
same scheme must be followed with the exception that the dimension of the right singular 
subspaces is higher and their structure now depends on several generators. 

For AT > 4 it has been showed that within each right singular subspace structure, there 
could exist a continuous infinity of states not connected through invertible local operators. 



(13.40) 
□ 
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Additionally, up to this continuous degree of freedom within each right singular subspace 
structure, we have found an upper bound for the number of classes on + 1 qubits in terms 
of the number of classes of N qubits. 

13.4 The N = 4 classes 

Here we give, in Table 113.21 and without explicit calculation, our result for the = 4 
entanglement classes under SLOCC according to our method. The specific details of the 
calculation may be found in |LLSS07a] . 

Agreeing to consider each structure of 2U as a single entanglement class, we have found 
18 degenerate and 16 genuine classes (totally 34 classes), where permutation is explicitly 
included in the counting. Taking into account the permutation among the qubits, 8 genuine 
classes can be considered, recopiled in table 113.21 As expected, in most of the classes 
a continuous range of strictly non-equivalent states is contained, although with similar 
structure. 

This result allows us to predict that there will be at most 765 entanglement classes (per- 
mutation included) for 5-partite systems, 595 at most genuine and 170 at most degenerate 
(cf. |LLSS06j ). This classification stands up as a formidable task. 

As a final remark, let us conjecture that a possible connection with the MPS formalism 
is probable to exist. In this formalism (c/. [EckOSl IPGVWC06j and multiple references 
therein) any pure state is written as 

^ tr(4^1...4-])e,,®---®e,^, 

ji...ijv 

SO that adjoining a further {N + l)-th qubit amounts to adjoining a further aJ^^^^ matrix 
in the trace giving the coefficients. In the analysis carried out above, this last added qubit is 
equivalent to increase the dimension of the right singular subspace dimWN — > dimSHAr+i = 
2 X dimWiy- Our conjecture is that the structure of 2)3 at should be read from the properties 
of the N matrices A^^'^\ so that the succesion of structures of Wn should run parallel to 
that of the matrices A^l^\ . . . , 4^^' . 
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Class (an) 


Canonical States 


Name 


Notation 


span{000, 000} 


|0000) + 


GHZ4 


2Hooo,ooo 


span{000,Ofc^'} 


loooo) + |1100) + |1111) 
0000) + 1101) + 1110) 


- 


SSJooo.Ofc* 


span{000, GHZ} 


|o^v?V) + |iooo) + |1111) 




WoOO,GHZ 


span{000, W} 


|1000) + |0100) + |0010) + |0001) 


W4 


2nooo,w^ 


span{Ofc^',Ofc*} 


|0000) + |1100) + AilOOll) + A2IIIII) 
0000) + 1100) + Ai 0001) + Ai 0010) 

+A2I1101) +A2I1110) 


$4 




span{Oi*, Oj^'l 


|0(^00) + \0(f)l^) + |1000) + |1101) 
0(/)0V') + |O01O) + 1000) + 1101) 






span{Ofe*, GHZ} 


|0(/)^') + |1000) + |1111) 






span{GHZ, W} 


|0001) + |0010) + |0100) + + ll^V') 




^GHZ,W 



Table 13.2: Genuine entanglement classes for four qubits 



Chapter 14 

Conclusions 



In this Thesis I have obtained a series of results related to quantum entanglement. The 
interest of entanglement is manifest, both from a fundamental point of view in quantum 
mechanics, or in order to process and transmit information with quantum systems more effi- 
ciently than with classical ones. Its study is thus very relevant, and here I have contributed 
with some developments that I believe may help in the understanding of this mysterious 
physical property of quantum systems. In order to do this, I have mainly followed three 
research lines related to three different aspects of entanglement: 

• Entanglement and Relativistic Quantum Theory. 

— Dynamics of momentum entanglement in lowest-order QED |LLS06b] . 

We have analyzed the momentum entanglement generation among two electrons 
which interact in QED by exchanging virtual photons. We show that surprisingly, 
S matrix theory produces pathological results in this case: the entanglement 
in M0ller scattering would be divergent for incident particles with well-defined 
momentum. In order to manage with these divergences, that would be physical 
(entanglement is a measurable magnitude, with a physical meaning) , we made the 
calculation for electrons with Gaussian momentum distributions which interact 
for a finite time. The divergences disappear, but, remarkably, the attainable 
entanglement would not be bounded from above. 

— Generation of spin entanglement via spin-independent scattering |LL06j . 

Here we have considered the spin entanglement among two or more identical 
particles, generated in spin-independent scattering. We show how the spatial 
degrees of freedom act as ancillas creating entanglement between the spins to a 
degree that will depend in general on the specific scattering geometry considered. 
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This is genuine entanglement among identical particles as the correlations are 
larger than merely those related to antisymmetrization. We analize specifically 
the bipartite and tripartite case, showing also the degree of violation of Bell's 
inequality as a function of the scattering angle. This phenomenon is unrelated 
to the symmetrization postulate but does not appear for unlike particles. 

Relativity of distillability [LMDSOGj . 

In this work we have studied the Lorentz invariance of usual magnitudes in 
quantum information, like the degree of entanglement or the entanglement dis- 
tillability. We introduce the concepts of relativistic weak and strong isoentangled 
and isodistillable states that will help to clarify the role of Special Relativity 
in the quantum information theory. One of the most astonishing results in this 
work is the fact that the very separability or distillability concepts do not have 
a Lorentz-invariant meaning. This means that a state which is entangled (dis- 
tillable) for one observer may be separable (nondistillable) for another one that 
propagates with a finite v < c speed with respect the first one. This is an 
all-versus-nothing result, in opposition to previous results on relativistic quan- 
tum information, which showed that a certain entanglement measure was not 
relativistically- invariant (but always remained larger than zero). 

— Dirac equation and relativistic effects in a single trapped ion |LLSS07b] . 

Here we have introduced a method for simulating Dirac equation, a quantum- 
relativistic wave equation for massive, spin-^ particles, in a single trapped ion. 
The four-component Dirac bispinor is represented by four metastable, internal, 
ionic states, which, together with the motional degrees of freedom, could be con- 
trolled and measured. We show that paradigmatic effects of relativistic quantum 
mechanics unaccesible to experimental verification in real fermions, like Zitterbe- 
wegung, Klein's paradox, Wigner rotations, and spontaneous symmetry breaking 
produced by a Higgs boson, could be studied. 

• Continuous- variable entanglement. 

— How much entanglement can be generated between two atoms by de- 
tecting photons? |LGRC07| . 

We have proved that in experiments with two atoms an arbitrary degree of en- 
tanglement between them may be reached, by only using linear optics and post- 
selection on the light they emit, when taking into account additional photons as 
ancillas. This is in contrast to all current experimental proposals for entangling 



two atoms, that were only able to obtain one ebit. 

— Spin entanglement loss by local correlation transfer to the momentum 
|LLS06a] . 

We have shown the decrease of the initial spin-spin entanglement among two 
s = ^ fermions or two photons, due to local correlation transfer from the spin to 
the momentum degree of freedom of one of the two particles. We explicitly show 
how this phenomenon works in the case where one of the two fermions (photons) 
traverses a local homogeneous magnetic field (optically active medium), losing 
its spin correlations with the other particle. 

— Schmidt decomposition with complete sets of orthonormal functions 
|LL05aj . 

We have developed a mathematical method for computing analytic approxima- 
tions of the Schmidt modes of a bipartite amplitude with continuous variables. 
In the existing literature various authors compute the Schmidt decomposition 
in the continuous case by discretizing the corresponding integral equations. We 
maintain the analytical character of the amplitude by using complete sets of or- 
thonormal functions. We give criteria for the convergence control and analyze 
the efficiency of the method comparing it with previous results in the literature 
related to entanglement of biphotons via parametric down-conversion. 

— Momentum entanglement in unstable systems |LL05bj . 

We have analyzed the dynamical generation of momentum entanglement in the 
decay of unstable non-elementary systems described by a decay width T. We 
study the degree of entanglement as a function of time and as a function of T. 
We verify that, as expected, the entanglement grows with time until reaching 
an asymptotic maximum, while, the wider the decay width F, the lesser the 
asymptotic attainable entanglement. This is a surprising result, because a wider 
width is associated to a stronger interaction that would presumably create more 
entanglement. However, we explain this result as a consequence of the fact that 
for wider width the mean life is shorter, so that the system evolves faster (during 
a shorter period) and can reach lesser entanglement than with longer mean lives. 

Multipartite entanglement. 

— Sequential quantum cloning [DLL"*" 06] . 

Not every unitary operation upon a set of qubits may be implemented sequen- 
cially through successive interactions between each qubit and an ancilla. Here 
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we have analyzed the operations associated to the quantum cloning sequentially 
implemented. We show that surprisingly the resources (Hilbert space dimension 
D) of the ancilla grow just linearly with the number of clones M to obtain. 
Specifically, for universal symmetric quantum cloning we obtain D = 2M and 
for symmetric phase covariant quantum cloning, D = M + \. Moreover, we ob- 
tain for both cases the isometries for the qubit-ancilla interaction in each step 
of the sequential procedure. This proposal is easily generalizable to every quan- 
tum cloning protocol, and is very relevant from the experimental point of view: 
three-body interactions are very difficult to implement in the laboratory, so it is 
fundamental to reduce the protocols to sequential operations, which are mainly 
two-body interactions. 

— Inductive classification of multipartite entanglement under SLOCC 

Here we have proposed an inductive procedure to classify A^-partite entanglement 
under stochastic local operations and classical communication (SLOCC) when 
the classification for — 1 qubits is supposed to be known. The method relies in 
the analysis of the coefficients matrix of the state in an arbitrary product basis. 
We illustrate this method in detail with the well-known bi- and tripartite cases, 
obtaining as a by-product a systematic criterion to establish the entanglement 
class of a pure state without using entanglement measures, in opposition to what 
has been done up to now. The general case is proved by induction, allowing us 
to obtain un upper bound for the number of entanglement classes of A'^-partite 
entanglement in terms of the number of classes for — 1 qubits. We also give 
our complete classification of the A^ = 4 case. 
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The Schmidt decomposition 



In this appendix we briefly review the Schmidt decomposition procedure [SchOGl lEKQSj 
IPBPOOj to express an arbitrary bipartite pure state as 'sum of biorthonormal products'. 

A.l Finite-dimension Hilbert space 

We begin with an arbitrary pure, normahzed state, pertaining to a finite dimension Hilbert 
space, associated to a bipartite system of subsystems Si and 5*2 



m=0 71=0 

where ttiq < uq with no loss of generality, and {|m)}, {|n)} are two orthonormal bases 
associated to Si and ^2 respectively. 

Expression (jA.ip does not show the degree of entanglement of the state \ip), because 
one cannot tell in principle whether this state is a product state (in which case the degree 
of entanglement would be zero) or not (in which case the state would be entangled). To 
address this question, and to quantify the degree of entanglement ('how much'), it is useful 
to express the state in its Schmidt decomposition. To do this, firstly we construct the 
density matrix of IV') 



and now we take traces over the largest subspace (between the subspaces related to Si and 
S2), in this case, the subspace associated to subsystem ^2 {mo < tiq by hypothesis). 



mo no 




(A.l) 




(A.2) 




(A.3) 



/'=0 



mm' 



,'=0 n"=0 
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Now we diagonalize the matrix Mmm' = ^m'm = Sn''=o ^i'nn"C^,^t,, obtaining the (trans- 
posed) matrix V of eigenvectors and the eigenvalues {Am}: 



mo 



m'=0 



Thus the eigenvectors \ipm) of M, resuh 



|V'^^> = Yl ^rnm'lm'), (A.5) 



m'=0 



and we can write (fX2|) . expressing |m) as a function of iV'm ) via , in the form 

mo no mo mo 

E E C'mnC;;,^, Yl OmlV'^.^') E ^™"'m'(V'i;l|®|n)(n'|. (A.6) 
mm'=Onn'=Q m"=0 m"'=0 

Now we define the (unnormalized) states 

. mo no 

E Y^mm'C^'nln), (A.7) 

m'=On=0 

and thus it fohows 

mo mo 

m=0 m'=0 

Expression (jA.Sh represents the density matrix of the pure state 

mo 



m=0 

f2) (2) 

The last step is to correctly normalize the states \Tpm getting the orthonormal states \iprn ) 

|V^^2))^ 1 I^W). (A.IO) 

V ^m 

Finally, the Schmidt decomposition we have obtained is 

mo 
m=0 

where {|V'm^)}, {|v^m^)} are orthonormal bases by construction and thus ()A.lip is a decom- 
position of {ip) in diagonal biorthonormal terms. 
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A quick inspection of (jA.lip tells us whether the state It/') is entangled or not. If the 
decomposition has just one term, the state is a product: not entangled. If the decomposition 
has more than one term, the state is entangled. The situation is different for some specific 
cases that we do not review in this appendix, like for a pair of identical particles jSCK^Olj 
IESBL02"1 IGM04] (see Chapter [2]) or whenever superselection rules are present |SVC04a[ 
ISVC04bj . 

A suitable magnitude to quantify the degree of entanglement of the state is the 
entropy of entanglement 

mo 

5 = - ^ A^loggA^. (A.12) 



m=0 



For a product state, 5 = 0. For an entangled state, S" > 0, and the more entangled is 
the greater is S. 



A. 2 Infinite-dimension Hilbert space 

In this section we review the Schmidt decomposition procedure for continuous-variable 
states. 

Now we consider a continuous pure bipartite state, pertaining to an infinite dimension 
Hilbert space, associated to a bipartite system, of the form 

|^)= J dpdqf{p,q)al^{p)a\^^{q)\0,0) (A.13) 
(\\f{p,q)f^ J dpdq\fip,q)\^ <ocy 

where aj^-^ (p) , a^^^-^ (q) are the creation operators of a particle associated to the subsystems 

51 and 5*2 which form the system, p and q are continuous variables associated to 5i and 

52 respectively, which can represent one dimensional momenta, energies, light frequencies, 
or the like. In general, p G (oi, 61), q £ (02, 62). 

The amplitude f{p,q) can then be expressed, via the Schmidt decomposition, as 'sum 
of biorthonormal products' in the form 



f{p,q)=Y,VKi^'aHp)^i?KQ), 
n 



(A.14) 
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where ipn , ipn and An are solutions of the integral eigenstate equations 

r dp'K^^\p,p')^^^\p') = \n4^\p)^ (A.15) 

r dq'K(^\q,q')^lji^\q') = A„V'i')(g), (A.16) 

and K(^\ 

are given by 

kW(p^p') ^ r dqf{p,q)np',q), (A.17) 

Ja2 

K(''\q,q') ^ r dpf{p,q)r{p,q'). (A.18) 
J ai 

A. 3 The singular value decomposition 

We include the relevant properties of the SVD of an arbitrary matrix and suggest the 
interested reader to consult e.g. [HJOl] for a comprehensive analysis of this decomposition 
with the corresponding proofs. The set of m x n complex matrices will be denoted as usual 
by -A4m,n(C) = ■Mm,n and the group of unitary matrices of dimension k will be denoted by 
U{k). The main result can be stated as 

Theorem 6. (Singular Value Decomposition) Let Q G M-m,n- Then Q can always be 
decomposed as 

Q = VT.W\ (A.19) 

where V £ U{m), W £ U{n) and S G M^, n is a diagonal matrix with non-negative entries, 
i.e. Sjj = aidij, with i = 1, . . . ,m, j = 1, . . . ,n and fX/t > for all k. 

The columns of V and W and the positive entries of S receive a special name: 

Definition 7. The columns of V = [vi V2 ... Vm] (resp. W = [wi W2 ... Wn]) are the left 
(resp. right) singular vectors of Q. The positive entries ofT, are the singular values of Q. 

Notice that with this definition any m x n dimensional matrix will have m left singular 
vectors and n right singular vectors; since the relevant singular vectors will be those asso- 
ciated to non-null singular values, we agree, as usual, on referring as singular vectors only 
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to the latter, i.e. to those Vk and Wk for which cjfc > 0. Another common convention is the 
decreasing order of the singular values in the diagonal of S: cji > ct2 > • • • > 0. 

The singular vectors are highly nonunique or equivalently there always exist another 
unitary matrices V and W such that Q = VTiW\ where these new unitary matrices depend 
of the former V and W and the multiplicities of each singular value jHJ91| . However this 
fact has not been exploited in Chapter [13l 

One of the main consequences of the SVD is that the rank of a given matrix Q coincides 
with the rank of S, i.e. with the number of positive singular values, which, in turn, coincides 
with the dimension of the subspace generated by the left (or right) singular vectors. This 
is the basis to the analysis of entanglement of a pure state upon its coefficient matrix in a 
product basis performed in Chapter [T31 
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In this appendix we briefly review the no-cloning theorem of quantum mechanics and two 
types of approximate quantum cloning (to a certain fidelity) : Universal symmetric quantum 
cloning, and phase-covariant quantum cloning. For a thorough review of the field, see 
[SI(^An5] . 

B.l No-cloning theorem 

Quantum mechanics forbids to exactly copy quantum states while leaving unperturbed the 
original state |WZ82t IDie82| . This is a fundamental property lying at the very core of 
quantum mechanics, and proofs can be found based in unitarity and also in linearity. 

Theorem 7. (No-cloning theorem) No quantum operation exists that can duplicate with 
fidelity 1 an arbitrary quantum state. 

Linearity-based proof: 

Proof. We proceed with a reductio ad absurdum by considering the 1 — > 2 case (two clones). 
We begin by considering that perfect cloning can be realized by a unitary operation such 
that 




(B.l) 



where \R) is the blank state into which the clone would be produced, and the internal 
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state of the cloning machine. In particular, for two orthogonal states |0) and |1), it is verified 

|0) \R) \M) ^ |0) |0) |A^(0)), (B.2) 

|1) «) \R) \M) ^ |1) |1) (B.3) 
But then, due to linearity, it is verified 

(|0) + |1)) \R) \M) ^ m\M{0)) + (B.4) 
which is incompatible with 

(|0) + |1))(|0) + \1))\M{0 + 1)) = (|00) + |01) + |10) + |11))|A^(0 + 1)). (B.5) 

This implies that Eq. (jB.ip may hold for the states of an orthonormal basis, but not for 
their superpositions. □ 
Unitarity-based proof: 

Proof. We proceed with a reductio ad absurdum by considering the 1 — > 2 case (two clones) . 
We begin by considering that perfect cloning can be realized by a unitary operation such 
that 

mR)\M) ^ um\R)\M)) = m^Mm, (b.6) 

for a certain state |^'). But for a different input state |$) of the cloning system, it should 
be verified 

mR)\M) ^ um\R)\M)) = \<^m\Mm- (b.t) 

Taking the inner product of Eqs. ()B.6p and ()B.7p . and imposing unitarity, we have 

(^1$) = (B.8) 

But Eq. (jB.Sp can only hold (for non-orthogonal states |^) and \^)) in case (^|<&)(A^(^)|A^(<I>)) 
1, which is impossible given that |^') and \^) are different states. Eq. ([R8l) may only hold 
in case |^') and |$) are orthogonal, but that would only allow to clone, as in the previous 
proof, the states of an orthonormal basis, and not their superpositions. □ 
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B.2 Optimal approximate cloning 

Since the seminal work by Buzek and Hillery |BH96j , which obtained the optimal Symmetric 
Universal Quantum Cloning Machine (SUQCM) for the 1 — > 2 cloning of qubits, a lot of 
research has been done in order to obtain the optimal unitary operations that clone an 
arbitrary quantum state maximizing the fidelity (for a thorough review, see [SIGA05] ). 
Here we review the optimal SUQCM |BH961 IGM97| for 1 — > M cloning of qubits, and 
the economical phase-covariant symmetric quantum cloning |DM03l IBDM05] . which has 
applications in cryptography. These are mainly the two cases considered in our proposal of 
sequential quantum cloning exposed in Chapter [121 

Definition 8. A Quantum Cloning Machine (QCM, unitary operation) is called universal 
if it copies all the states with equal fidelity. On the other hand, it is called symmetric if at 
the output all the clones have the same fidelity. 

B.2.1 Symmetric universal quantum cloning machine: 1 ^ M case for 
qubits 

The general formula for the 1 — > M cloning of qubits was provided by Gisin and Massar 
|GM97] . The unitary operation associated to the cloning machine is 

M-l 

\i;)(^\B) ^\GMMm ■■= Y.aj\{M - j)^l^,j^l^^)s(^\{M - j - l)r ,jr^)s, (B.9) 

j=0 

where \B) denotes the initial blank state, aj = ^ m{m+^ \ {M — j)ip, jip-^) s denotes the 
normalized completely symmetric state with (M — j) qubits in state (p and j qubits in state 
0"*". It must be noticed that the presence of M— 1 additional so-called anticlones is necessary 
in order to perform this cloning procedure with the optimal fidelity. The anticlone state ip* 
refers to the fact that it transforms under rotations as the complex conjugate representation. 

B.2. 2 Economical phase-covariant symmetric quantum cloning: 1 M 
case for qubits 

This is a case of state-dependent cloning. Here the motivation is to clone at best an arbitrary 
state of a subspace of the whole Hilbert space of one qubit. In fact, the phase-covariant 
QCM is defined as the QCM that copy at best states of the equator (x — y) of the Bloch 
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sphere, i.e., those states of the form = l/\/2(|0)+e*'^|l)). We have only focused upon the 
cases where no anticlones are needed (hence the term economical). Under this assumption, 
imposing the purity of the joint state, the number of clones M must be odd [DM03]. The 
optimal machine in this case was obtained by D'Ariano and Macchiavello [DM03] . and is 
associated to the operation 



\B) 



1 

72 L 



|(A; + l)0,A;l)5 + e*^|A;0,(A; + l)l)5 



(B.IO) 



where k = {M - l)/2. 



Appendix C 

Matrix-Product States 



In this appendix we review the protocol |Vid03| for expressing a multipartite pure state in 
its matrix-product form (MPF, cf. |Eck051 IPGVWCOGj and multiple references therein), 
already present in spin chains |AKLT87] , classical simulations of quantum entangled systems 
|Vid03j and density-matrix renormalization group techniques |VPC04j . 

We begin with a multiqubit state pertaining to a Hilbert space Tif^, expressed in the 
computational basis, 

1 1 

= 12"'Y1 Cii...ijk- ■ -in)- (C.l) 

«1=0 in=0 

Our aim is to obtain its MPF, i.e., to express it in the form 

i^)=e( E nl]c^AMa,m%aA'^um%a,■■■nn]^-^\h (c.2) 

where {r[l], . . . , r[n]} are two- and three-tensors and {A[l], . . . , X[n — 1]} are vectors. The 
indices ik and Ok take values in {0, 1} and {!,..., x}i respectively, being x = max^ XA- Here 
XA is the rank of the reduced density matrix of the partition A : B of the multipartite 
state (fUn) . 

In order to obtain Eq. ()C.2p we first compute the Schmidt decomposition (SD, see 
Appendix [X]) of state (jC.ip according to the partition 1 : 2 • • • n 

m = E^w-ii'^S)i'^*L'r"') = E nmAMaAh)\^^^r^), (c.s) 

where in rhs of last equality we have expanded each Schmidt vector |^>a|) = 
in terms of the computational basis for the qubit 1. 
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The next step is to expand eacli Schmidt vector \^ai" ) in the computational basis for 
qubit 2 

lC-"') = El^2)|r|;t')- (C.4) 

Next we express each vector {rj^^l'^^) in terms of the at most x Schmidt vectors {|*J'a2 "^)}a2=i 
of the bipartition 12 : 3 ... n and the associated Schmidt coefficients A[2]q2 

Now we insert Eq. (|C.5p in Eq. (|C.4p and the resulting expression in Eq. (|C.3p , and obtain 

IV'>= E ^[l]j,^^A[l],,^[2]j,^^„2A[2]„,|il.2)|<^>&•"')• (C-6) 

Proceeding iteratively in this way, by making the n — 1 SD associated to the successive 
bipartitions of state (|C.ip . we arrive straightforwardly to Eq. (jC.2p . This is basically the 
protocol developed in [Vid03| . 
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